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Abstract. Converting a pioneering idea of V. V. Sokolov et ai (60] to a geometric 
object, we introduce a well-defined notion of linear matrix operators in total deriva- 
tives, whose images in the Lie algebras 9 of evolutionary vector fields on jet spaces are 
closed with respect to the commutation, 




The images generate involutive distributions on infinite jet bundles J°°(7r) over fibre 
bundles tt, and the operators induce Lie algebra structures [ , ]a on f2 = dom A/ ker A. 
We postulate that the operators are classified by the vector and covector transforma- 
tion laws for their domains dom A, which are related to sections of 7r by Miura substi- 
tutions. The gauge of dom A may be independent from transformations of the images 
im A. This is a generalization of the classical theory with dom i ~ imi C sym £ for 
recursion operators for integrable systems £, and with dom A ~ imA 3 cosym£ for 
Hamiltonian operators (here the Miura substitutions are the identity mappings in both 
cases, and the gauge transformations arc uniquely correlated). In particular, recursion 
operators A £ Endk sym £ that satisfy (0) are solutions of the classical Yang-Baxter 
equation [Aipi, A1P2] = A(Upi, <P2\a) for the Lie algebra sym£ 9 tpi, (p 2 . 

If, for r linear differential operators Aj : fl — > g with a common domain f2, pairwise 
commutators of their images hit the sum of images again, we endow the spaces A = 
k • Ai of the operators with a bilinear bracket 



that satisfies the Jacobi identity. A class of such algebras is given by operators that 
generate Noether symmetries of hyperbolic Euler-Lagrange systems of Liouville type; 
we calculate explicitly the operators Ai and the structural constants c^- . Thus we give 
an exhaustive description of higher symmetries for all 2D Toda chains associated with 
semi-simple Lie algebras of rank r, completing the results of 43, 48, 53, 561 [60] . 

The bracket [ , ] of operators is anti-symmetric for the domains f2 composed by 
vectors and is symmetric for f2 with covector transformations. In the latter case, we 
reveal a flat non-Cartan affine connection in the triples — > g such that symmetric 
bi-differential Christoffel symbols are encoded by the structural constants c*j- , and 
such that completely integrable commutative hierarchies 21 C imA are the geodesies. 

We demonstrate that the notion of Lie algebroids over infinite jet bundles J°°(Tr) 
does not repeat the construction over finite-dimensional base manifolds, when condi- 
tion Q is fulfilled by the anchors. To correlate them with Lie algebra homomorphisms 
A: £1 g, we introduce formal differential complexes over Lie algebras (f2, [ , }a)- 
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Introduction. Let M n be a smooth n-dimensional real manifold and ir : E n+m — > M n 
be a fibre bundle. The space of linear differential operators over R on the infinite jet 
bundle J°°(7r) — ^ M n contains operators A in total derivatives whose images in the 
(7°°(J°°(7r))-module r(7r^(vr)) = r(?r) <g>c~(M™) C 00 ^ 00 ^)) of evolutionary derivations 
are closed with respect to the commutation, 

[im A, im A] C im A. (1) 

We say that such operators A are Frobenius in view of the fact that their images 
constitute involutive distributions on J°°(7r). Each Frobenius operator transfers the 
Lie algebra structure of evolutionary vector fields to the bracket [,}a on the quo- 
tient of its domain by the kernel and determines an IR-homomorphism of Lie algebras. 
The Hamiltonian operators that correspond to the Poisson bi-vectors with vanishing 
Schouten bracket are an example, among many others; e.g., the 2D Toda systems 
^Toda = {u l xy = exp(K l jui)} associated with the Cartan matrices K determine a class of 
Frobenius operators that take values in the spaces sym £roda of their (Noether) higher 
infinitesimal symmetries. 

The main problem we address is the coordinate-independent definition of operators 
that satisfy (PQ), and we perform the analysis of their algebraic and geometric properties. 
We conclude that, up to a special case, there are two principal classes of Frobenius op- 
erators, which generalize the recursions sym £ — > sym £ and the Hamiltonian structures 
cosym£ — > sym£ for integrable systems £, respectively. The essential distinction of 
generic Frobenius operators A from the recursion or Hamiltonian operators is that the 
coordinate transformations in the domains of A can be not correlated with the trans- 
formations in their images. In this sense, the two arising classes of Frobenius operators 
correspond to vector and covector-type transformation laws for their domains f^ 1 (£„.), 
which are defined by ir and another fibre bundle £ over M n . 

For example, if £ = 7r and hence fi 1 (7r 7r ) = r (71^(71")), then Frobenius endomorphisms 
R G EndjR r(7r^ C) (7r)) of generating sections <p of evolutionary vector fields are not only 
recursions for the Lie algebra {q{tt), [, ]) of derivations. Indeed, they induce nontrivial 
deformations [, ]r of the standard structure [, ] on Y^n^in)) via the classical Yang- 
Baxter equation 

[Rsp 1 ,Rtp 2 ] =R([tpi,<P2U), fi,f2 e r(7r^(7r)). (2) 

Frobenius operators of second kind, which generalize the Hamiltonian structures, are 
assinged in this paper to hyperbolic Euler-Lagrange systems of Liouville type. To 
do that, we give an exhaustive description of higher symmetries for these integrable 
systems. In particular, our explicit formula (11051) for the operators is valid for all 
2D Toda chains associated with root systems of complex semi-simple Lie algebras of 
rank r. These operators determine commutative hierarchies of KdV-type evolution 
equations which are Noether symmetries of the Euler-Lagrange systems. The operators 
specify factorizations of higher Hamiltonian structures for the hierarchies. 

We introduce a Lie-type algebra structure on spaces of linear differential opera- 
tors in total derivatives. We consider the linear spaces A = ©[ =1 k • A4 of operators 
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A4: fi 1 (^ 7r ) — > r(7r^ (7r)) with a common domain ft 1 ^) such that the sums of images 
of Ai are closed under the commutation of vector fields, 



[E, im ^E/ m ^-] ^ J2 k imAk - ( 3 ) 



By setting [A iy Aj](p,q) := [A^p), Aj(q)\ for p,q G fi 1 (^ 7r ) and 1 < i,j < r, we endow 
the linear spaces A with a Lie-type bracket 



that satisfies the Jacobi identity. The bi-differential structural constants • : il 1 (£„•) x 
r2 1 (^ 7r ) — > fi 1 (^ 7r ) are skew- symmetric for vector pre-images of evolutionary fields and 
are symmetric for covectors that constitute the domain f2 1 (^ 7r ) of A+. In the latter case, 
we calculate the constants explicitly. This algebraic problem and our approach are 
different from Sato's formalism [20] of pseudo differential operators because, e.g., the 
compositions of Frobenius operators fi 1 (^ 7r ) — > T^n^n)) are not defined even in the 
Hamiltonian case, and hence they do not constitute multiplicative associative algebras 
with unit. Therefore the Lie-type algebra structure relies heavily on the intrinsic ge- 
ometry of the jets J°°(7r) as the base manifolds for all fibre bundles. The bracket we 
propose, in particular, for Frobenius recursions R is different from the Richardson- 
Nijenhuis bracket [, ] that determines their Lie super-algebra structure [301 §5.3] by 
using the composition of operators, although we use similar geometric techniques. 

In this paper, we pass from the (wide) category of fibre bundles over M n to a more 
narrow category of bundles and horizontal modules over the infinite jets J°°(7r). Then 
we investigate which objects survive and which there appear anew. We conclude that 
the operators that obey (CEJ) are closely related to Lie algebroids over smooth finite- 
dimensional manifolds M n . We show that Lie algebroids over J°°(7r), with the anchors 
given by Frobenius operators, can not be defined as straightforward generalizations of 
the finite-dimensional case, but this is achieved with a much more refine construction. 
To this end, for each Frobenius operator A, we generate a formal differential com- 
plex over Lie algebra (fi 1 (£„.), [, ]a) and suggest to determine the Lie algebroids using 
representations of the differential in this complex through homological vector fields. 

On the other hand, we discover that Frobenius operators determine a generalization 
of the affine geometry such that bi-differential Christoffel symbols T^- are encoded by 
the structural constants of the algebras A. We prove that Frobenius operators A 
determine flat connections in the triples T^n^-n)), A) consisting of two Lie 

algebras and a morphism, and we interpret commutative hierarchies in the images of 
Frobenius operators as the geodesies. This is different from the classical realization of 
integrable systems as geodesic flows on infinite-dimensional Lie groups [2J. 

Our main result is the following. We formulate a well-defined notion of linear differen- 
tial operators whose images are closed w.r.t. the commutation, see (CD), and we associate 
infinitely many such operators to Euler-Lagrange hyperbolic systems of Liouville type. 
Here, in particular, we solve an old problem [321 S3], [56] in geometry of Leznov-Saveliev's 
2D Toda chains related to semi-simple complex Lie algebras: We describe their infin- 
itesimal symmetries by an explicit formula and calculate all commutation relations in 
these Lie algebras. 




(4) 
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The paper is organized as follows. First in section [T] we outline our basic concept. 
Then in section [2] we summarize important properties of the Poisson structures for 
completely integrable evolutionary systems. We introduce the notion of nondegenerate 
operators in total derivatives and formulate a conjecture that simplifies the search for 
Poisson structures of PDE. 

In section |3] we define the Frobenius operators, establish the chain rule for them, 
describe the differential Frobenius complex, and construct flat connections in the triples 
(fi 1 (^ 7r ), 5(71"), A) of two Lie algebras and a morphism. Also, we describe an inductive 
method for the reconstruction of the Sokolov brackets on the domains of nondegenerate 
Frobenius operators. We illustrate it using the dispersionless 3-component Boussinesq 
system fl5*U|) . see [TH1 Elj, that admits a family of Frobenius operators floT|) and yields 
the Frobenius recursion flHTT) . 

In section H] we propose the definitions of the linear compatibility and the strong 
compatibility of Frobenius operators. We endow the linear spaces of both linear and 
strong compatible Frobenius operators with the Lie-type algebra structure. We show 
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that the Magri schemes provide commutative examples of such algebras. The structural 
constants encode the bi-differential Christoffel symbols, which are transformed under 
reparametrizations by a direct analogue of the rules for the classical connection 1-forms. 

In section [5] we assign the algebras of Frobenius operators A to Liouville-type Euler- 
Lagrange systems (in particular, to the 2D Toda chains associated with the root systems 
of semi-simple Lie algebras pUl H2J [56]). To this end, we describe the generators ip = 
A(-) of their higher symmetry algebras and calculate all commutation rules. As a by- 
product, we find the Hamiltonian structures for KdV-type systems related to these 
hyperbolic Darboux-integrable equations. 

Finally, we discuss open problems in the theory of Frobenius structures. 



Let us fix some notation; the language of jet bundles is contained, e.g., in [51 HH US 
I5T] . In the sequel, everything is real and C°°-smooth. By J°°(vr) we denote the infinite 
jet space over a fibre bundle n: E n+m — > M n , we set vr^: J°°(7r) — > M n , and denote by 
[u] the differential dependence on u and its derivatives. Put J-(ir) = C°° (J°° (n)) , which 
is understood as the inductive limit of filtered algebras [5] and hence each function from 
J~"{tt) depends on finitely many coordinates on J°°(7r). The tt- vertical evolutionary 
derivations 8. v = ^2 a D a (ip) ■ d/du a are described by the sections cp G T (n^n)) = 
r(7r) <g>c°°(M n ) C°° (J°° (7r)) of the induced fibre bundle tt^tt). The shorthand notation 
for this jF(7r)-module is x(7r) = r(7r^ (7r)). For all ip such that £ V ,(V') makes sense, 
the linearizations £^ are defined by £^\<p) = £tp(ip), where <p G x(ir). Denote by 
A n (7r) the highest 7r-horizontal forms on J°°(7r). For any ^ r (7r)-module f) = r(7r^(£)) = 
r(^)(S>(7oo( M n)C 00 (J 00 (7r)) of sections of an induced bundle over M n , we use the notation 
£„. = 7r^ (^) and denote by (j = Hom^) (f, A n (7r)) the dual module. Examples of £ will 
be given later. 

Frobenius operators map f — * x(7r), where the module f is one of the following 
restrictions of 1} C r(£ 7r ) onto the image of a Miura differential substitution w = 
w[u] : J°°(7r) -> r(f): We have that either 

f = X (0L : joo( 7r )^ r ( ? ) or f — MOL; J°c(7r)-»r(5)- 

In particular, £ = vr for recursion operators R G End]Rx(7r), and f = x(ir) for Hamil- 
tonian operators; here we set w = id: r(7r) — > T(£) in both cases. 

Denote by q(tt) = (x(7r), [, ]) the Lie algebra of evolutionary vector fields £^ with 
the standard bracket of sections. Let (fl 1 }a) denote the Lie algebra that is 

isomorphic to f/ ker A as a vector space and is endowed with the Lie bracket [, ]a by a 
Frobenius operator A: f — > q(tt). Definitions are discussed in detail in sections [2] and [31 

1. Basic concept 

The problem of construction and classification of operators in total derivatives that 
satisfy (pQ) was suggested first in [60J. The operators A were regarded there as non- 
skew-adjoint generalizations of the Hamiltonian operators, whose images are closed 
with respect to the standard Lie algebra structure on the modules of evolutionary 
fields. The operators appeared in [60] in local coordinates in the context of the scalar 
Liouville-type equations £l = { u% xy = F(u,u x ,u y ; x,y),l < i < m}, whose infinite 
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groups of conservation laws [f(x, [w]) dx] + [f(y, [w]) dy\ are differentially generated 
by finitely many densities 

Wi, . . . , w r G ker D y \ , W\ , . . . , w f G ker D x | . (5) 

Let us recall the main motivating example that begins the theory of Frobenius operators. 

Example 1 (The Liouville equation). Consider the scalar Liouville equation 

^Liou = {U xy = exp(2W)}. (6) 

The differential generators w, w of its conservation laws [rj] = [f(x, [w]) dx] + [f(y, [w])dy] 
are 

w = U X — U xx and w = U y — U yy (7) 

such that D y {w) = and D x {w) = by virtue (=) of ^Liou and its differential conse- 
quences. The operators 

□ = U x + \D X and □ = U y + \D y (8) 
factor higher and Noether's symmetries 

„ = □(««,[«])), ^ = n(^M) ; t=a(to,w», 9c = a(?%*M) 

of the Euler-Lagrange equation (jSJ) for any smooth 0, and 7i, TC. Note that the 
operator □ = \D~ l o is obtained using the adjoint linearization of w, and 

similarly for □. 

Each of the images of (|HD is closed w.r.t. the commutation such that 

[□0), □(<?)] = □(£n( p )(g)-£n( ff )(p)+{{p, ?}}□), here {{p,g}}n = D x (p)-q-p-D x (q), 

and same for □; the evolutionary derivations £(.) are given in (1281) on p. [13] in local 
coordinates. The symmetry algebra sym^Lkm — imD + imQ is the sum of images 
of ©, and the two summands commute between each other, [imd, imO] = on ^Liou- 
Therefore, 

[im □ + im □, im □ + im □] C im □ + im □. (9) 

The Frobenius operator □ factors higher symmetries of the potential modified KdV 
equation 

£pmKdV = {U t = -\U XXX + Ul = D{w)}, (10) 
whose commutative hierarchy is composed by Noether's symmetries <pc ^ im(Do5/(W) 
of the Liouville equation (J6j). The operator □ factors the second Hamiltonian structure 
B 2 = □ o A 1 o □* for £ pmK dv, here A 1 = D x l . 

The generator w of conservation laws for Suou provides the Miura substitution (J7J) 
from £ pm Kdv to the Korteweg-de Vries equation 

£xdv = {w t = ~\w xxx + 3ww x }. (11) 

The second Hamiltonian structure for £i<dv is factored to the product A 2 = □* o B\ o □, 
where B\ = D x is the first Hamiltonian structure for the modified KdV (see dia- 
gram (11131) on p. [321 . The domain of the Frobenius operator □ contains sections of the 
cotangent bundle cosym£ K dv to the KdV equation. The bracket {{ , }}□ on the domain 
of □ is equal to the bracket {{ , }}^ 2 induced on the domain of the operator A 2 (which 
is Hamiltonian and hence its image is closed under commutation) for £i<dv- 
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Example [T] is reproduced for all m-component 2D Toda chains 



associated with semi-simple complex Lie algebras [121H3] . For example, all the 2D Toda 
chains (112j) with the matrix K symmetrizable by a vector a (that 
no summation over i) admit the conserved density w l = (ku x ,u x ) /2 — (a,u xx ). At the 
same time, the chains (fT2"j) may admit other conserved densities. In the fundamental 
paper [SB], A. B. Shabat et al. proved the existence of maximal (r = f = m) sets (J3J 
of conserved densities for ^Toda if and only if the matrix K is the Cartan matrix of 
a root system for a semi-simple Lie algebra of rank r. Further, Ref. [55] contains an 
explicit procedure that yields special systems of jet coordinates for ( {121) ; with respect to 
them, all coefficients of the characteristic equation D y (w) = on ^Toda become linear, 
whence its first integrals w 1 , . . ., w r are obtained. The differential orders of w,w grow 
as r grows, and the formulas are big already for the Lie algebra G 2 , see [HH EI]. We 
claim that the orders of w l with respect to Shabat 's momenta m,,- (see below) coincide 
with the gradations for the principal realizations of the basic {i.e., simplest nontrivial 
highest weight) representations of the corresponding affine Lie algebras. 

The generators ip = D(0(x, [w]) of higher symmetry algebras for Liouville-type equa- 
tions are factored by matrix operators □ in total derivatives [8j [29J EE]. For the Euler- 
Lagrange Liouville-type systems, the operators provide Noether symmetries (pc with 
4> = 5H(x, [w])/Sw, see also [531 EUj ES], whence the description of all symmetries fol- 
lows. For this class of hyperbolic equations, which incorporates f[T2"]) . commutative Lie 
subalgebras of higher Noether symmetry algebras yield completely integrable KdV-type 
hierarchies pUl [28], [25]. The generators w,w of conservation laws for £l induce Miura's 
transformations w = w[u] between the KdV- and modified KdV-type hierarchies upon 
w and u, respectively. The operators □ factor higher Poisson structures for the evolu- 
tionary systems and, moreover, prescribe the nonlocalities that arise in these structures. 
For instance, the first integral w 1 for systems (fT2l) with a symmetrizable matrix K yields 
the class ip = (u x + A ■ D x ) (0(x, [w 1 ])) of symmetries for f^oda, here A* = Y^jLi this 
always yields the second Poisson structure for KdV in the upper-left corner (^4fc) 11 of 
the Hamiltonian operators determined by the entire operator □ (see [25] for details). 
Several formulas for the operators □ were known from [131 [56] for 2D Toda chains ( [121 ) 
associated with semi-simple Lie algebras of low ranks, but not in the general case. 

Under assumption that the densities ([5]) are known, we obtain the explicit for- 
mula fllQ5p for these operators □. We derive all the commutation relations for sym^; 
they are encoded by bi-differential operators c£ that act on 0's. We prove that the 

r-tuples (f) obey the variational covector transformation laws (j) i— > <fi = [(^& )*] 1 (0) 
under reparametrizations w = w[w] of ([5]). Therefore the operators □ assigned to 
Liouville-type systems £l are natural examples of well-defined Frobenius structures. 

The second motivation to study linear differential operators subject to ([!}) is much 
more abstract, c.f. [40J. Let us compare geometry of ODE and of PDE from the fol- 
lowing viewpoint. For a linear k-space V, one can study representations q — > Endj c (l/) 
of Lie algebras $j on V and thus endow the linear spaces of endomorphisms with Lie 
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algebra structures. On the other hand, let T be a commutative associative algebra 
with unit, e.g., a k-algebra = C°°(M n ) or an J^.^-algebra F(ir) = J°°(ir)), 
and let d: T — > be a derivation. Consider two left JF-modules P, Q and the space of 
linear differential operators Diff (P — > Q) . A question: Are there any natural algebraic 
structures on this linear space? If P — Q (e.g., P = Q = sym£ of a differential equa- 
tion £, and we deal with recursion operators, see El]), then one has the associative 
composition A o B and the formal commutation A o B — B o A for A, B G P, but what 
else? And what if P ^ Ql In this paper, we give an affirmative answer on the above 
question whenever T = ^"(vr). 

In geometry of PDE [51 00], US], classical constructions such as the differentials of 
Hamiltonians, bi- vectors A J^, etc., appear as zero-order termg^ in Taylor expansions 
of the variational derivatives, variational bi-vectors, or variational Poisson-Nijenhuis 
structures [T7], [22j H2J- Likewise, we discover that the Lie-type brackets, which we 
introduce on linear spaces A C CDiff(f — > x(7r)J of Frobenius operators, encode bi- 
differential Christoffel symbols of flat connections. At the same time, we show by 
counterexample H] on p. [23] that Frobenius operators are not the anchors of Lie algebroids 
over the jet spaces J°°(7r). 

We begin with the classical Hamiltonian formalism for ODE on finite-dimensional 
manifolds M. Bearing it in mind, we pass to the PDE setting and formulate the 
assertions. 

Let V G r(/\ 2 (TM)) be a bi- vector with vanishing Schouten bracket = 0. 

Using the coupling ( , ) on TM x T*M and a nondegenerate Poisson bi- vector V, one 
transfers the Lie algebra structure [, ] on Y{TM) to [, \p on T(T*M) 3 if)i,il>2 and 
obtains the Koszul-Dorfman-Daletsky-Karasev bracket |9J 

[ipu^p2]v = Lr^iCVb) - Lp^(V>i) + d(P(^!,^ 2 )), (13) 

here L is the Lie derivative. By [33], the bracket (TIB"]) is uniquely defined if, first, 

• [ , \p is a derivation of C°°(M) with respect to the C°°(M)-module structure of 
r(T*M), that is, if we have that 

bPuf- = • V>2 + / • bPuHv (14) 

for any / G C°°(M) and Vi,^2 e V(T*M), and second, if 

• [dh, dh 2 ] = -d h 2 } P holds for any hi,h 2 G C°°(M). Then ^(T*M, [ , ] P ) ^ 

(TM, is a Lie algebroid [_59J with the morphism V over the smooth manifold M. 

The de Rham differential ddR on f\'(T*M) is defined in the complex over the Lie 
algebra (T(TM), [, ]) by using Cartan's formula. If the Poisson bi- vector V has the 
inverse symplectic two-form V~ l such that "P _1 [x,y] = [V^^-x, V~ l y\-pi then the dif- 
ferential ddR is correlated with the Koszul-Schouten-Gerstenhaber bracket [,Jj> on 
/\'(T*M) by d dR = [P' 1 , -jv The differential d dR on /\*(T*M) is intertwined pj [37] 

1 To facilitate the exposition, in appendix[21we summarize the notation for Hamiltonian evolutionary 
PDE and its analogy with the structures for ordinary differential equations. However, it is a very 
delicate matter to restore the full parallel between finite-dimensional manifolds M and horizontal 
bundles over J°°(7r) with its own Cartan connection. Essentially, theorems either are converted to 
definitions or become false. This is the case of Lie algebroids, which we discuss in sections EOl and 151 
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with the Poisson differential d-p = {V, •] on f\'(TM) by 



(A 'P)([P~ 1 ,*h) + [PAA p )Wl = > V*e/\ (T*M). 



The trivial infinitesimal deformations 



[x, y) N := [Nx, y) + [x, Ny] -N([x, y\) 



_d_ 

dA a=o 



[^^(11)] 



(15) 



of the standard Lie algebra structure [ , ] on the tangent bundle TM 3 x,y over smooth 
manifolds M were described in [35] using the recursion operators N: T(TM) — ► T(TM). 
If the Nijenhuis torsion 



for an endomorphism N vanishes, then the Lie brackets [, ] N * obtained by iterations of 
the Nijenhuis recursion N are pairwise compatible. 

The Nijenhuis and Poisson structures (N, V), whenever satisfying two compatibility 
conditions p21 [35] , generate infinite hierarchies of pairwise compatible Poisson struc- 
tures N fc o V, k > 0, which means that linear combinations AiN fcl o V + A2N fc2 o V 
remain Poisson for any Ai : A2. This assertion [35], which is valid for finite-dimensional 
manifolds M, admits a straightforward generalization to the infinite-dimensional case 
when the base manifold is the jet space J°°(tt) over a fibre bundle it: E m+n — > M n , 
see P3B 22J, such that the concept of Poisson-Nijenhuis' structures is applicable ver- 
batim to evolutionary PDE [17]. Also, this construction yields the Lie algebroids 



Let P be a Hamiltonian operator for a system of differential equations £ = {w 1 = 0, 
. . ., w r = 0}. The left-hand sides w % belong to some T (it) -module f of sections of an r- 
dimensional fibre bundle over J°°(7t), and the operator P: f — > k(tt) takes sections from 
the dual of f to evolutionary fields. In practice, the equations u\ = fi[u] in determined 
evolutionary systems £ are labelled by the dependent variables u l in the bundle tt, which 
establishes the isomorphism f ~ k(tt). However, the system £ as a geometric object 
in J°°(tt) admits arbitrary reparametrizations w = w[w] for its components. They are 
not anyhow correlated with admissible changes u = u[u] of the dependent variables. 
Summarizing, we see that the Hamiltonian operators for PDE are well-defined under 
unrelated transformations of their domains f and images k(tt). Let us formalize this 
property for a wider class of operators. 

Frobenius linear differential operators in total derivatives, with images closed under 
the commutation, are well defined as follows. Let it and £ be fibre bundles over M n , 
let w be a fibre coordinate in £, and construct the infinite jet bundle ^ : J°°(£) — > 
M n . Consider the jF(£)-module of sections x(£) = r of the induced fibre 

bundle an d denote its A"-dual by x(£). Suppose further that there is a Miura 

substitution J°°(tt) — > r(£) which embeds both .F(£)-modules into r(7r^ (£)). (For 
instance, the substitution determines a system of differential equations £ = {w 1 ^} = 
0, . . . , w r [u] = 0}.) We denote the substitution by the same letter w, because from now 
on we take the restrictions of and x(£) onto its image. 

• Frobenius operators of first kind are A: x(£)| — > x(tt). Under any diffeomor- 
phisms u = u[u] : J°°(it) — > T(tt) and w = w[w] : J°°(£) — > T(£), the operators A of first 



[N,N] fn (x, 2 /) = [Na;,N l /]-N([x,y] N ) 



(16) 



(( 
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kind are transformed according to 



A ^ A = o A o /(™) 



w=w[u\ 
u=u\u] 



(17) 



• Frobenius operators of second kind are linear mappings A: x(£)| — > x(7r). For 
any differential changes of coordinates u = u[u] and w = w[w] in 7r and £, respectively, 
the operators obey 

A^A = tt ] oAo{^y , y (18) 



u=u\u\ 



• Finally (the degenerate case), if no gauge of T(£) is allowed in a given setting, then 
Frobenius operators A: r(7r£j(£)) — ► x(7r) are transformed by A i— ► A = £^ o ^4| u=u [~] 
under « = 

Frobenius operators of first and second kind generalize linear recursion operators 
R G End]R x(7r) and Hamiltonian operators P : k{ji) — > *c(tt), respectively. The domains 
of A may not be composed by sections of the tangent bundle 7r^(7r) to J°°(7r) or, 
respectively, of its dual, which was constructed in jUJ (see also [TjJl 122] )• Thence 
Frobenius operators are described neither by the Poisson bi- vectors nor by the recursion 
(1, l)-tensors. 

For the same reason, the coupling ( , ) on Q x is missing. However, Frobenius 

operators A: f C T(£ n ) — ► x{tx) transfer the Lie algebra structure [, ] for evolution- 
ary vector fields to [,]a on the quotients ^ 1 (^ 7r ) = r(£,r)| /kerA of the respective 
jF(7r)-submodules f of r (£„-). The Koszul bracket [,]a is defined by A(Jfa, fa] a) = 
[Afa,Afa] for any fa, fa £ ^ x (£"-)• Thus the operators determine the Lie algebra 
IR-homomorphisms 

A: (^(UUa) ->(«(*),[,])■ (19) 
For any 7r-horizontal g-module K there is a flat connection 

V A : Der Int (fi 1 (^),0(7r)) - Der( fl (7r), if) (20a) 
that lifts inner derivations of f2 1 (^ 7r ) to K- valued derivations of q(tt) by the formula 

The connection is f2 1 (^ 7r )-linear, V^ x r^ .j A = A(<fi) x VA w.r.t. the Lie multiplications 
(p x ■ = [(p, -] A by any <p e fi 1 ^), here A(<p) x ■ = [A(<f>), ■]. The connection V A is flat 
due to the Jacobi identity. Thence the commutative hierarchies whose flows belong to 
the image of A are the geodesies with respect to V A . 

The 'variational' analogue [, ] A of the Koszul-Dorfman bracket (fT3]) equals 

[fa,fa]A=^Af 1 (fa)-^A^p 2 (fa) + {{fa,fa}}A, fajfaeft 1 ^), (21) 

where the evolutionary derivations are (j28j) . Generally, the bi-differential term {{ , }}^ 
is neither a Lie algebra structure nor a cocycle. The chain rule (I78p for the brackets 
{{ , }}a and {{ , }}aoA determined by two Frobenius operators A and AoA follows from 
the evolutionary summands in (j2ip . 
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An explicit formula for the Zhiber-Sokolov bracket^ {{ , }}p in the Hamiltonian case 
P: x(7r) — > x(7r) is given by (j4"fi]) . We reconstruct {{ , }}□ for a class of Frobenius 
operators^ □: cosym£ K dv — * sym£ pmK dv C sym£ Toda that factor higher symmetries 
of the 2D Toda chains ( Tl~2l) and the Hamiltonian structures for the associated KdV- 
type evolution equations [ID]. Thus Zhiber-Sokolov's pioneering idea in [6D] to study 
operators with property (00) as generalizations of the Hamiltonian operators P: x(tt) — > 
g(7r) was motivated by the examples □: f — > x(7r) of a different geometric origin with 
f % x(7r). The bracket {{ , }}r for Frobenius recursions i? G End^ x(7r) is unknown. 
For Frobenius operators A that are nondegenerate in the sense of we formulate 
an inductive procedure for reconstruction of {{ , }}a- 

The Leibnitz rule (1141) . which is the first axiom for the Lie algebroids, is lost for the 
Koszul brackets (1211) even in the Hamiltonian case f = x(7r). Therefore the under- 
standing of Frobenius structures as the Lie algebroids over the infinite jet spaces, which 
is one of our goals, is essentially nontrivial and requires the construction of a formal 
differential complex over (f2 1 (^ 7r ), [, ]a)- Namely, to each operator (fT9l) we assign the 
Frobenius complex^ 

0(7r ) ^ Hom M (fi 1 (e 7r ),g(7r)) ^4 Hom^/yV^BM) 

iHom R (/\ 3 fi 1 (^),g(7r)) . (22) 

The first inclusion in (f2"2"]) consists of the commutation S/Afo with fixed elements ip G 
g(7r), the second arrow is the composition of the Koszul bracket [, ]a with A, and 
the third arrow calculates the right-hand side of the Jacobi identity. The homological 
vector fields Q 2 = that encode the Frobenius complex (|2"2"1) assign the Lie algebroid 
structures over the jet spaces J°°(7r) to Frobenius operators using the approach of [59J. 

We say that the Frobenius operators are linear compatible if their arbitrary linear 
combinations retain the same property (00). Families of N linear compatible operators 
Ai, . . ., An induce deformations of ( fl9j) such that the Sokolov brackets obey 

JV 

i=l % — 1 

The linear compatibility extends the frames L^IR- A{ to the linear spaces A = ® i M.-Ai. 

The operators A%, . . ., An are called strong compatible if, by ((3D, the commutators 
of evolutionary fields £>A k {-) i n their images belong to the sum of these images. That is, 
for any i,j G 1, . . . , N and p, q G fi 1 (£„.) = f/ Pj^ ker A, we have 

N N 

[Mp), A^q)] = A^Eam^)) - M*Mi)(P)) +J2M r iM<i)) e E im ^- ( 23 ) 

k=l l=\ 



2 The bracket {{ , }}p is not the Adler-Kirillov-Kostant bracket applied in 10J to construction of 
the Drinfel'd-Sokolov systems. 

3 We shall denote by □ the Frobenius operators of second kind that are associated with Liouville-type 
systems of PDE. 

4 The presence of the two evolutionary summands in (|21[) is the reason why we stay in the Horn- 
category and do not pass to the more narrow CDiff-category of the operators in total derivatives. 
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In this way, the Koszul brackets (|2ip in the domains of each Ai merge to the collective 
decomposition of the commutators. For example, the commutation closure @ shows 
that the operators (jSJ) are strong compatible and that 

rg n = D x <g> 1 - 1 <g> D x , rg n = D y ® 1 - 1 <g> D y , 

T° D = D V ®1, T° u = -Kg)D x , r^=-l®D y , T% a = D x ®l, (24) 

where the notation is obvious. Note that q) = rjEj^Qo, q) = T^ n (q,p) = (q,p) = 

on ^Liou for any p(x, [w]) and q(y, [w]). The bi-differential symbols satisfy T^(p, q) = 
— r^(g,p) and other normalizations. They are transformed by a suitable variant of the 

classical formula T = gTg^ 1 + dg g^ 1 for connection 1-forms T and reparametrizations g. 

If the closure condition ([3]) is fulfilled for N Frobenius operators A±, . . ., An, then we 
take [A*, Aj](jp, q) : = L4j(p), ^-(g)] as the definition of a Lie- type bracket (jlj) between 
them. Each operator A £ spans the one- dimensional algebra that is described by T £ u = 
{{ , }}A e , see (!2~T|) . We extend the bracket of operators onto the linear span ®^ =1 K • A t 
if, in addition, the Frobenius operators are linear compatible. The constants are bi- 
differential Christoffel's symbols for the flat affine connection {V^, 1 < t < N} defined 
in ( l20i) . The Magri schemes for completely integrable systems j9l [35] yield examples of 
commutative algebras of Frobenius operators. 

Relaxing the normalization Y k u = {{ , }}a ( ■ 8\ but retaining (jSj), we arrive at a wider 
class of linear differential operators A\ and Lie-type brackets (jlj) on the spaces A = 
• A. Namely, the image of a particular operator may hit the images of other 
operators under the commutation. Now we formulate the assertion. 

Theorem. Let k be a real constant nondegenerate symmetric (mxm)-matrix. Consider 
a hyperbolic Euler-Lagrange system £ L = {5C/5u = 0} which, in a suitable system 
of coordinates, is determined by a Lagrangian C — [L dxdy] with the density L = 
— \{ku x ,u v ) — Ki,(u;x,y). Set m = dL/du y . Suppose further that the system £l is 
Liouville-type and admits r conserved densities w[m] = (w 1 , . . . ,w r ) G ker D y \ . Let 
them be minimal such that / G kerZ)^ implies / = f(x, [w]). Introduce the operator 

□ = (4 m T- ( 25 ) 

Then we claim the following: 

(i) All (up to x «-> y) Noether symmetries (pc of the Lagrangial £ for £l are 

ip c = □ ^ ^ J for any H. 

(ii) All (up to x y) symmetries ip of the system £l are 

up = n(<f){x, H)) for any <f) = (<f>\ . . . , <f) r ) G f = H0\ w=w[m y 

(iii) In the chosen system of coordinates, the image of the operator □ is closed with 
respect to the commutation in the Lie algebra sym^L- 

(iv) Under a diffeomorphism w — w[w], the r-tuples <fr are transformed by 

Therefore, under any reparametrization u — u[u] of the dependent variables u = 
t (u 1 , . . . ,u m ) in equation £l, and under a simultaneous change w = w[w], the 
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operator □ obeys ([15]) . Consequently, the operator □ satisfies (JT]) in any system 
of coordinates, and hence □: MO | m - m [ TO ] — * s Y m ^-L is a Frobenius operator of 
second kind. 

(v) The operator 

V = U* oU^Y oU: k(i)\ I , 

V m / v ' ' I ui=-(ii[m[«J] y ^ ' \w=w[m[u\\ 

is Hamiltonian. 

(vi) The bracket {{ , }}□ on the domain f of the operator □ satisfies the equality 

{{,}}□ = {{, }}v- 

Its right-hand side is calculated explicitly by using formula (H6l) that is valid 
for Hamiltonian operators V. This yields the commutation relations in the Lie 
algebra sym£ L . 

(vii) All coefficients of the operator V and of the bracket {{ , }}□ are differential 
functions of the minimal conserved densities w for Si,. 

2. Hamiltonian operators for PDE 
In this section we recall necessary definitions and introduce standard notation, which 

follows El EH 1221 SQl EU- 

2.1. Symmetries and conservation laws. Let n: E n+m — > M n be a vector bundle 
over an n-dimensional manifold M and J°°(7r) be the infinite jets over this bundle. 
By definition, set 77^ : J°°(7r) — > M n . Put f[u] = f(u,u x , . . . ,u a ), here u denotes the 
components of a vector t (u l , . . ., u m ) and \a\ < oo. 

On J°°(tt), there is the Cartan distribution C of n-dimensional planes that project 
without degeneration onto TM n under tc^^. The distribution is spanned by the total 
derivatives D x i, 1 < i < n, whose action on ^(n) = C°°( J°°(tt)) is defined by restriction 
onto the jets joo(s) of sections s G r(7r), 

3oo (s)(DAf)) = a? (*»(*)(/)), Vs G rW ' V/ G *W (26) 
In coordinates, we have 

„ d \ ,■ d 

D xi = — + > u 3 a , i. -. 27 

The Cartan connection ' ' : A \— > D T i is flat. 

Denote by x(vr) the J r (7r)-module r(vr^(7r)) = T(7r) ® c °°(m») C°°( J°°(tt)). The sec- 
tions ip = (if 1 , . . . , ip m ) G r(7r^ (7r)) of the induced fibre bundle "look like" sections 
s(x) = t (s 1 (x), . . . , s m (x)) G r(7r), but their components (fi l (x, [u]) G J 7 {it) are functions 
on the jet space J°°(tt), see also (1124ft . In coordinates, the evolutionary derivations 
with sections if G x(n) ar^ 

£, = £W)-^. (28) 



5 The notation £ v makes no confusion with £ = {u t = ip(x, [u])}, because almost always we identify 
such evolutionary systems with 7r-vertical derivations E v £ D v ( J°°(tt)). A synonimic notation Q v or 
d v is used for in the literature. 
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Their invariant definition is 18.^, D x i\ = for £ v G D v (J°° (tt)) . We shall use the 
permutability of evolutionary fields (12H1) with total derivatives (12T1) many times. 

The commutators of evolutionary fields (12"S1) endow x(tt) with a Lie algebra structure, 

[£^,£j = £[ VltVa ], where [<^i,^ 2 ] = £^(^2) - £ V2 0i) (29) 

which is defined by the componentwise action. In what follows, we identify evolutionary 
vector fields 8,^ with their generating sections (p. We denote by q(tt) the Lie algebra 
(x(7r), [, ]) with the bracket (1251 . and we call it the standard Lie algebra structure 

on x(it). By definition, a recursion R for g(7r) is a linear differential operator x(7r) — > 
x(7r) in total derivatives. 

The linearizations (the Frechet derivatives), 

— —f ■ Ar, ip G r(7r^(0) for some bundle f over M™, (30) 

are correlated with (128]) by = ^ (</?)• 

Under a change w = of fibre coordinates in ir, the generating sections of (I2"H1) 
are transformed by t— > y5 = £^ \(p). Hence the recursion operators obey 

R^R = tfoRolf (31) 

u=u[u\ 

Let dh = ddR(M«) be the Cartan lifting of the de Rham differential on M n , 

n 

d h = ^dx i <8> D x *. (32) 

i=l 

Let A*(7r) be the JF^) -module of i-th horizontal forms on J°°(7r). Take the cohomology 
A*(7r)/{imd/ l } w.r.t. the horizontal differential d^ and denote it by H l (n). So, elements 
[H(x, [w])dvol(M n )] G H n (it) are equivalent if they differ by exact terms dhf], where 
i] G H n ~ 1 (7i). The equivalence classes 7i G H n (n) of the highest horizontal cohomology 
will be called the Hamiltonians or the Lagrangians, whenever it is appropriate. 

The Cartan differential on J°°(tt) is de = ddR(j°°(7r)) — d^. The pair of differentials 
(dh, dc) generates the variational bi-complex and its C-spectral sequence. Passing to the 
horizontal cohomology at each term of the variational bi-complex, construct the first 
term Ei(n) of the C-spectral sequence. Choose a restriction of dc onto H n (7r), which 
is a corner of the bi-complex, to be the Euler variational derivative. (That is, apply 
the Cartan differential dc and then perform multiple integration by parts, which does 
not alter the equivalence class modulo dh-) Denote by E the Euler variational deriva- 
tive (S/Su 1 , . . . , 5 /5u m ) with respect to u = t (u 1 , . . . , u m ). Recall that the operator of 
linearization for the image of the variational derivative is self-adjoint, 

^(H)=t E{n) , HeH n i7i). (33) 

For any H = [#dvol(M n )] G H n (ir) we have E(H = £* H (1). Consequently, under any 
change of variables u = u[u] : </°°(7r) — > r(7r), a section ip = E(7i) in the image of the 

variational derivative is transformed by ip 1— > iJj = [(4^)*] ^VO L-^j- 

For any fibre bundle £ over M n , the J~(^)- modules r(7r^ (^)) will be called horizontal. 
For any horizontal module f, denote by f the dual T{ti) -module f = Rom^^ (f, A n (7r)) . 
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We study differential equations in jet bundles, considering them together with all 
their differential-algebraic consequences, which we assume to exist. Let 8 C J°°(tt) 
be an equation, and let the differential ideal 1(8) of the equation be generated in 
coordinates by S = {F 1 = 0, . . . , F r = 0}, here F % G In other words, the r-tuple 

t (F 1 , . . . , F r ) is a section of a horizontal bundle f = r(7r^ (^)) for £ with r-dimensional 
fibres over M n . Note that £f- x(n) — > f. 

Remark 1. There are two particular examples of £ that, unfortunately, often lead to 
confusion and tempt us to ignore the general setting: 

• Determined autonomous evolutionary systems 8 = [u\ = (p 1 (x 1 , . . . , x n , [u]) , 
1 < i < m}. Then f ~ x(n). 

• Euler-Lagrange systems 8el = (E(£) = 0, C G H n (ir)\. Then f ~ x(n) = 
Hom^ (7r) (x(7r),A"(7r)). 

In both cases, the number of equations r equals the number of unknowns m, which 
can be used misleadingly to index the equations. We emphasize that, in principle, 
the reparametrizations of F G f that determine the system 8 = {F = 0} and of the 
dependent variables u in J°°(ir) are not correlated at all. 

In this paper, we consider 

• Hamiltonian operators P : k(8) — > x(8) for determined evolutionary systems 8 
(see section 12.21) , 

• Noether non-skew-adjoint operators A: k(8) — > x(8) for the same class of equa- 
tions (see section 12.41) and inverse Noether operators uj : x(8) — > k(8) (see Ex- 
ample [H]) , 

• Frobenius operators □ : k(8\) — > sym£ 2 that factor symmetry flows 82 on Euler- 
Lagrange equations 8el (see section EJ), 

• Frobenius recursion operators R G EndiR x(ix) (see Example [HI), and 

• Frobenius operators A: f C r(7r^ (£)) — > x(n) of general nature (see section [3]). 

Let 8 = {F = 0} be a system of differential equations given by a section F G f. Let 
a class [rj\ G fP" 1 ^ be a conservation law. In other words, the divergence 

d hV = V(F) (34) 

vanishes on the ideal 1(8) of algebraic-differential consequences of 8, and this is realized 
by r operators V = (V 1 , . . . , V r ) in total derivatives. The coefficients of V belong to 
C°°(8) and the operators take values in the module A n (7r) of highest horizontal forms. 

Represent the n-th horizontal form dhf] as V(F) = /l,V(F)), where V coincides 
with V in any local coordinates but takes values in f instead of A" (71-). Here 1 stands 
for the r-tuple (1, . . . , 1) and ( , ) : f x f — » A n (7r) is the coupling. Integrating by parts 
and introducing the adjoint operator V*, we obtain (1, V(F)) = (V*(l), F) mod imd^. 
The section V*(l) G f is called the generating section of a conservation law [77] for the 
equation 8 = {F = 0}. By construction, if) = V*(l) G f has as many components as 
there are equations in the system. Under a reparametrization F = A(F) of equations 
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that constitute the system £, where A is a lineal differential operator in total deriva- 
tives, the section if) G f is transformed by if) i— > if> = A* (if)). The generating sections if) 
of conservation laws are solutions of the equation £ F (if)) = on £\ this follows from the 
formulas = E(d^r/) = 1*^^(1) = t F {if>) + £^(F) = 0, where the second summand is 
set to zero by restriction onto the consequences of £ = {F = 0}. 

We conclude that if) G >c(£) for determined evolutionary systems (this is why the 
notation if) G cosym£ is used), and if) G x(£el) f° r Euler-Lagrange systems £el- 

Lemma 1 ([10]). Let £ = {F = ut — f{x, [u]) = 0} be an evolutionary system. For any 
conservation law [77] = [pdx + - ■ •] such that dh(rj) = V(F) = (V*(l), F), the generating 
section if) = V*(l) of [77] is the 'gradient' if) = E(p) of the conserved density p. 

For the Euler-Lagrange equations £el = {E(£) = 0}, the generating sections if) G 
f (£el) of conservation laws have a geometric nature of symmetries. This is indeed so. 

Theorem 2 (Noether). Let £el = {E(£) = 0} be the Euler-Lagrange equation for a 
Lagrangian L G H n (ir). Then the evolutionary derivation 6^ is a Noether symmetry of 
the Lagrangian, 8><p(C) = 0, if and only if if is the generating section of a conservation 
law [rj\ such that d^r? = on £ EL . 

Lemma 3 ( [40] . see a proof in [28] ). The relation 

E(£„(p)) = £ v (E(p))+^(E(p)) 

holds for any (p G x(tc) and p G A n (7r). 

Consequently, any Noether symmetry (fc G sym c £^ of a Lagrangian is a symmetry 
of the Euler-Lagrange equation £ EL = {F = 0} (and Lemma [H] shows also why the 
converse is not true). For this reason, mappings of generating sections of conservation 
laws to symmetries of differential equations will be called Noether in what follows. 
We have demonstrated that the identity transformation is Noether for Euler-Lagrange 
equations £el = {E(£) = 0}. 

In the sequel, we apply Lemma [3] for other purposes as well; namely, it will allow us 
to induce a bracket ( 155bl) on f starting with a bracket (I55al) of conservation laws for 
evolutionary systems £. 

2.2. Poisson structures. In the preceding section we have recalled that the image of 
Euler derivative is dual to the module x(7r) of evolutionary fields. The standard concept 
of Hamiltonian dynamics for PDE (see [21 1221 HO]) is based on a class of Hamiltonian 
operators that map — > x(tt) and which induce a Poisson structure on H n (ii) due 
to the existence of a coupling ( , ) : k{ir) x x{ti) A n (7r). 

In this section we consider Hamiltonian operators. First we regard them as mappings 
of ^ r (7r)-modules, and then as Lie algebra homomorphisms. Next, we restrict these 
homomorphisms onto evolutionary systems and, finally, to hierarchies of systems. 

Now we recall necessary algebraic constructions. 

6 If the nondegenerate differential substitution A is nonlinear, e.g., F = F + F 2 , then tp becomes a 
nonlinear operator on the module f. This is a difficulty of the theory, see footnote^ on p.[2S] Hence 
only linear coordinate transformations are declared admissible in the module of equations; in the sequel, 
nonlinear changes of fibre coordinares will be allowed for horizontal bundles of other geometric nature. 
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For any vector space V over R, let A G Hom^A^ V, V) and V G Hom^A' V, V). De- 
note by A[V] G Hom R (A fc+ ^ 1 V, 7) the action A[-] : Hom^A* V, 7) -> Hom^A*"*"* -1 7, 
7) of A on V, which is given by the formula 

A[V](ai, . . . ,a fe+ z-i) = ^2 (- 1 )' JA (V(a (T (i), • . . , a a ^), a^j+i), . . . , a CT(fc+z _i)), 



where ai E V and 5^ C S'm denotes the unshuffles. The unshufHes are permutations 
such that er(l) < cr(2) < . . . < cr(/c) and a(A; + l) < . . . < cr(ra) for all cr G 5^; note that 
cr(i) is the index of the object placed onto the z-th position under the permutation, unlike 
in [21]. The Schouten (Richards on -Nijenhuis) bracked [A, V] G Hom^A^ 1 7 V) 
of A and V is P [27] 

[A, V] = A[V] - (-l^-W-^VfA]. (35) 

Thus the commutator [X, Y\ = [X, Y] of two evolutionary vector fields X, Y G q(tt) is 
skew-symmetric. The bracket [A\, A2} £ /\ 3 V of two bi- vectors is symmetric w.r.t. A\ 
and A2, but it is skew-symmetric w.r.t. its arguments ip2, ip3 G V*. 

Definition 1. A linear skew-adjoint [m x m)-matrix operator A: x(tt) — > x(7r) in 
total derivatives is called Hamiltonian if, for any 0,1,0,2 G H u (tt), the bi-linear skew- 
symmetric bracket!! 



{ , } A : #»(tt) x #«(7r) -> ^(tt), {H x ,H 2 }a ■= (E(Wi), A(E(W 2 ) 
satisfies the Jacobi identity 

X){{Wi,W2UWs}a = 0. (36) 

o 

By construction of the Poisson bracket { , }a, its equivalent definitions are 

{O u 2 } A = (Vx, A(^ 2 )> ~ £a(^)(Wi) ~ -£a(^)(W 2 ) (37) 

modulo (~) exact terms, here ^ = E(7ij). The bracket { , }^ endows H n (ir) with a Lie 
algebra structure over M. Two Hamiltonian operators Ai and A2 are called compatible 
if their linear combinations X±Ai + A2A2 are Hamiltonian as well. 

We postpone examples of Hamiltonian operators to section 12.41 where Noether op- 
erators will be considered. In example H2] on p. EH we investigate a wider class of 
Hamiltonian operators defined on a-modules over commutative k-algebras a. 

Lemma 4. Suppose that evolution equations £ = {F = u t — f[u] = 0, F G f} are 
enumerated by the dependent variables u, which means that the isomorphism f ~ x(7r) 
is being used. Then, under a fibre coordinate change u = u[u] in ir that preserves the 
evolutionary form of the equations £, a Hamiltonian operator A is transformed by 

A^A = ^oAo(i^y (38) 

u=u[u\ 

Formula (1381) is also valid for any Noether operator A: x(7r) — > x(7r). 



7 The definition of Schouten bracket is valid over any field k such that char Ik =/= 2. 

8 We respect Dirac's notation. Our choice of the signs in ([57)) is such that we multiply by bra- 
covectors {ip\ from the left and have the ket-vectors \ip) standing on the right. However, sometimes we 
use the reversed coupling ( , ) : f X f — * A n (ir) for convenience. 
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Hamiltonian operators A can be regarded [H [22] as the variational Poisson bi- vectors 
with vanishing Schouten brackets [A, A] = such that the Poisson bracket {Hi, 0.2} A — 
11 A, 7i 2 ], Hi] is a derived bracket in the sense of [32] ■ The compatibility condition for 
the Poisson bi- vectors is \Ai, A 2 ] = 0. 

By definition, put 

WV):=(^W)W ( 39 ) 

for any ip G x(tt), ip G f = P(tt^(£)) for some £, and an operator A G CDiff (f, >c(it)) in 
total derivatives. Note that £a,ip is an operator in total derivatives w.r.t. its argument <p 
and w.r.t. ip (but not w.r.t. the coefficients of A), and hence the adjoint l* A ^ is well 
defined. We emphasize that the notation £a,i/i is not the same as the linearization £am) , 
which was introduced in fl30l). 



Lemma 5 ([10]). Let A be a matrix operator in total derivatives as in ( l39i) . Then one 
has 

for any sections ipi, ip2 of a horizontal module over J°°(ti). 

Proposition 6 (A criterion of {A, A] = 0, |40j). A skew-adjoint operator A G CDiff (x(7r), 
x(7t)) in total derivatives is Hamiltonian if and only if the relation 

£a,^(AM) ~ ^ 2 (AM) = A{l\^i)) (40) 

holds for all ipi, ip 2 G x(7r). The r.h.s. of formula ( 14*01) is skew-symmetric w.r.t. ipi, ip2- 

The proof is informative in itself. It amounts to a straightforward calculation of 
the value of variational Schouten bracket [A, A] on three Hamiltonians Hi, H2, H3 G 
H n (n). Let ipi — E(Hi) be the respective gradients. The Jacobi identity {A, A] (Hi, H2, 
H3) = can be expressed as (b^i, 1P2), IP3) — 0, where b is a differential operator w.r.t. 
each argument and ( , ) : f x f — > A n (7r). Since ^3 G f is arbitrary, we have b^i, 1P2) = 
for all ipi,ip2 £ f- The calculation shows that h(ipi,ip2) is equal to the left-hand side 
minus the right-hand side of 



Lemma 7 ( [401 §7.8]). Consider an operator in total derivatives J G CDiff n\ (x(7r), P) 
which is skew-symmetric w.r.t. its I arguments that belong to k(ir) and which takes val- 
ues in an .F(7r)-module P. If for all Hi, . . ., Hi G H n (it) one has J(E(Hi), E(Hi)) = 
0, then J = 0. 

Lemma [7] implies that the Jacobi identity J(ipi,ip2,' l Pz) = can be verified for ele- 
ments ipi G imE only (in particular, this is done in the proof of Proposition [fj]) • 

Proof of Proposition® Let H a , Hp, and H^ be Hamiltonians. The Jacobi identity is 

{{H a , Hp} A, H-y}A + {{Hj3, H-y}A, H a }A + {{H-y, Ha) A, Hp} A = — 2_j ^Aiip^dAfyajjl/jp)) 

o 

= -^[(ZAfaMMa),^) + (A(E A( ^)(ip a )),^) + (A(ip a ),E A ^)(^))] = 0. (41) 

o 
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Consider the elements of the second sum, 

(A(E. A ^)(ip a )),ipp) = (V>a, A^Afa^a))) = -(A{ijjp)^A{^){^ a )) 
= -(A^J^A^))) = -(i; a (A(^)),A(^)) = (by 

= -<£^(A(^)),i4(^)> = -{A(^),e^A(^)). 
Substituting this back in ( I4T1) . we obtain 

o = - X)<( £ W^))W.tk) + [£ MM^MO) - J2( A ^ a ),U M^)) 

O 

= -((Ea^)^))^),^) - ((^A^ a ){A)){ip ),ipy) - ((Ea^sM))^),^)- 
Now set a = 3, /3 = 2, 7 = 1; thence we have 

= -(&M*)(A))(fo)M ~ {{^Am{A)){^)M - ((£a(w(A))(^),^>- (42) 
Consider the first summand, 

((£a Wi )(^))(^),^2> = ((W^i)))^2> = (A(^),^ 3 (^)) = (by LemmaED 

= -(W(^)).W' ( 43a ) 
Next, the second summand in ( 1421) is equal to 

(43b) 



Now consider the third term in the right-hand side of (j4"2"|) . 

((£a ( «(A))(^i)^3) = ((W^)),^)- (43c) 
Substituting (J4"3"j) in (J4"2"j) . we finally obtain 

{£a,^(A&i))M + (A^MM - ((e A ^(A(^)),^s) = 0, 

whence follows (|40p . The proof is complete. □ 

Actually, Proposition [6] states that images of Hamiltonian operators are closed w.r.t. 
the commutation. This is readily seen from ( I45bl) below, which contains the left-hand 
side of (j4"0|) . Note that, by Lemma [3, the assertion holds for Hamiltonian evolutionary 
vector fields £a(v>) which may not possess Hamiltonians TC such that ip = E(7Y). How- 
ever, the Lie algebra of Hamiltonian evolutionary vector fields £a(e(«)) £ fl( 7r ) that do 
possess the Hamiltonians Ti is correlated by a morphism of Lie algebras with the Lie 
algebra (H n (ir), { , } A ) as follows. 

Proposition 8 (0 §27]). The Euler derivative E and Hamiltonian operators A deter- 
mine the Lie algebra morphisms 

(fl»(7r), { , }a) ^ (*(*). [ , ]a) (0 (tt), [ , ]) , (44) 

such that 

4([V>i, ^2] a) = [A0i, ^1,-02 G x(vr) 
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ancO 

[E(W 1 ),E(W 2 )]^ = E[{Hi,H 2 }a), H l: H 2 G H n (n). 

The correlation between the Poisson bracket { , }a on H n (ir), the Koszul-Dorfman- 
Daletsky-Karasev bracket [, ]a on the quotient 

^(tt) : = x(»/kerA, 
and the standard commutator [, ] of evolutionary fields £a(-), see (l2l?l) . is 

[^i,^ 2 ]a = SagmM - £a(v> 2 )(^i) + {{V , i,<M}a = E((Vi,A(V» 2 ))), (45a) 

where 

A({{^i,^}}a) = (£a^(A))(^ 2 ) - (£a^ 2 (A))(^i). (45b) 

For A Hamiltonian, both ker A C (x(vr), [, ]a) and imA C (fl(vr), [, ]) are ideals in the 
respective Lie algebras. 

Formula ( |40l) provides the expression for the Sokolov bracket {{ , }}a, 

{{AMh = t A ^2), M2 e x(vr), (46) 

which is valid for Hamiltonian operators A if they are nondegeneratd^l in the sense 
of (1551) . see below. In coordinates, we have that for a Hamiltonian operator A = 
HEr^A-H and ^1,^2 ex(vr), 

{{AMYa = £(-1)" (a, o [J2 dM) ■ -q^-\ ) W) • 

cr,a r,/3 CT 

Another way of calculating Dorfman's bracket [,}a for any Hamiltonian operator A is 
based on the equivalence (l3"Tj) and further use of Lemma El see (155bl) in section 12.41 

The property A: x(£) — > x(£) of Hamiltonian operators remains valid for their 
restrictions onto the equations they determine. We have it as follows. 

Lemma 9 (@Dj). Let A be a Hamiltonian operator and consider a Hamiltonian evolu- 
tionary system 

ut = A(E(H)), H G H n {7i). (47) 
Then the operator A takes the generating sections if) = E(p) G of conservation 

laws [77] = [pdx+- ■ ■ } for ( 1471) to generating sections <p G of symmetries E v G sym£ 
for this system. 



Next, recall that the Schouten bracket (1351) of variational bi- vectors satisfies the 
Jacobi identity 

[[A 1; A 2 ],A 3 ] + [[A 2 , A 3 ], AJ + [[A 3 , A], A 2 ] = 0. (48) 

Hence the original Jacobi identity {A, A](?/>i, ip 2 , ^3) = for the arguments of A im- 
plies that Da = {A,-} is a differential, giving rise to the Poisson cohomology H\(n). 



9 Hence the image of the Euler operator E is closed with respect to the commutation. 
io A 

coordinate-free definition of the nondegenerate Hamiltonian operators means that they have 
trivial kernels on the ^"-coverings over the jet bundle J°°(7r), see [55] for the construction of A-coverings 
over PDE. 
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Obviously, the Casimirs c G H n 0n) such that {A,cJ = for a Hamiltonian operator A 
constitute H\0k). 

Theorem 10 (The Magri scheme IB]). Suppose [A U A 2 \ = 0, H G ^(tt) a 
Casimir of A\, and H\ (n) = 0. Then for any k > t/iere is a Hamiltonian Hk G H u (k) 
such that 

[A 2 ,n k -i] = [A u n k ]. (49) 

Put := Ax(E(Wfc)j. TTie Hamiltonians Hi, i > 0, pairwise Poisson commute w.r.t. 
either A\ or A 2 , the densities of Hi are conserved on any equation u% k = (fk, and the 
evolutionary derivations E lfk pairwise commute for all k > 0. 

We emphasize that Theorem [TOl holds for the Hamiltonians which belong to the linear 
subspaces S C H n (ii) spanned by the descendants of the Casimirs Ho G Ok). The 
commutativity of flows <p k and formulas (l3"Tj) imply that the phase volume J f s {ipi, A0ip2)) 
dipidip2 = is conserved on the subspaces S. 

The Magri scheme starts from any Hamiltonians Hk-i, Hk G H u 0k) that satisfy f j49l . 
but we always intend to operate with maximal subspaces of H n (n), and therefore we 
require [A^Ho] = such that the sequences of Hamiltonians can not be extended 
with k < 0. 

Proof of Theorem [TU1 The main homological equality (j49p is established by induction 
on k. Starting with a Casimir Ho, we obtain 

= {A 2 ,0j = [A 2 , [A^-Ho]] = -[A 1} {A 2 ,Hojj mod {A X ,A 2 \ = 

using the Jacobi identity (|48p . The first Poisson cohomology H\ 0k) = is trivial 
by assumption of the theorem, and hence the closed element [y4 2 ,7i ] in the kernel of 
{Ax, •] is exact: [A 2 ,7Yo] = [^i^i] f° r some H\. 
For k > 1, we have 

[A u [A 2 ,H k }j = -{A 2 , [Ai,Hkjl = -{A 2 , [A^H^jj = 

using the Jacobi identity (jUJ) for the Schouten bracket flHol) and by {A 2 , A 2 \ = 0. Hence 
[A 2 ,Wfe] = [v4 1 ,7i fe+1 ] by H\ 0k) = 0, and thus we proceed infinitely □ 

Correlated Magri's schemes for coupled hierarchies are further considered in section^ 
see diagram (j!13p . 

2.3. Lie algebroids. Hamiltonian operators determine Lie algebra homomorphisms ( 14*41) 
that map the Poisson bracket ( 1371) to the standard Lie structure ( 1291) on k{k). In 
the finite-dimensional case, this yields the Lie algebra homomorphisms P: Y(T*M) — > 
T(TM) to sections of the tangent bundles for smooth manifolds M. Recall that a very 
important class of the Lie algebra homomorphisms to T(TM) is then provided by the 
Lie algebroids. 

Definition 2 Q59J). A Lie algebroid over an n-dimensional manifold M n is a vector 
bundle Q — > M n whose sections TQ are equipped with a Lie algebra structure [, ] 
together with a morphism of vector bundles a: Q — > TM, called the anchor, such that 
the Leibnitz rule 

bPl,f-Ha = f- tyuH* + H^)f) ■ ^ (50) 

holds for any ^1,^2 G Ttt and / G L7°°(M"). 
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Equivalently, a Lie algebroid structure on is a homological vector field Q on HQ 
(take the fibres of Q, reverse their parities, and thus obtain the new super-bundle 
HQ over M n ). The homological vector fields, which are differentials on C°°(HQ) = 
T(/\*Q*), equal 

Q = y ia ?w-^ + \y j y ic U^ Q 2 = o, m 

where 

• (x a ) is a system of local coordinates near a point x G M n , 

• (y l ) are local coordinates along the fibres of Q and (y l ) are the respective coor- 
dinates on HQ, and 

• a(ej) = af(x) ■ d/dx a is the image under the anchor a and [e», e^] a = c^(x)ek 
give the structural constants for a local basis (e,) of sections TQ, respectively. 

The proof of equivalence is straightforward, see j5£l [33] for details and other equivalent 
definitions. 



Example 2. Consider the Cartan connection (1271) . and set the anchor a: C —>■ T(TM n ) 
to be a = 7TQO,* o Joo(so)*, which is the restriction joc(so)*(X) of horizontal fields X onto 
the jet joo(so) of a fixed section s £ r(7r) and further projection under tToq-. J°°(ir) — > 
M n . The corresponding homological vector field is the restriction j o(so)*(d/ l ) of the 
horizontal differential (I3"2l onto joo(so)- This is a generalization of a classical example 
Q = TM n , see 



Example 3. Let P: k(ir) — > x(7r) be a Hamiltonian operator and {, }p: H n (7r) x 
H n (7i) — > H n (-K) be the Poisson bracket. Consider the one- dimensional bundle over 
J°°(7r) with the fibre dx 1 A ... A dx n = d vol(M n ) at each point of M n . The equivalence 
classes of sections constitute an jF(7r)-module and are endowed with the Lie algebra 
structure { , }p. Set the anchor a := £p E- By ( J37j) . the Leibnitz rule (J50l) holds: 

{H 1 ,f-H 2 }p = f-{H 1 ,H 2 } P + a(H 1 )(f)-H 2 , e H n (ir), f e ^(tt). (52) 

The sections ?^ G H n (ir) are parameterized by the Hamiltonian densities /i G Let 
/i be the functional coordinates in Ilff n (7r), whose fibre is parity-dual to ft,. Then the 
homological vector field Q that encodes the Lie algebroid (H n (ir), { , }p) — > (x(tt), [ , ]) 
is equal to 

f r 1 ^ ■ 

Q = / ft ■ £p(E(Mvoi(M»)])) + ^ft 2 fti {[ftidvol(M n )], [ft 2 dvol(M")]} p ■ AF(tt). 

J^(Tr) L ^ aft J 

Proposition [H] states that a, which is defined by a Hamiltonian operator P, is a 
morphism of Lie algebras. This property was not included into the definition of Lie 
algebroids because, even in a general situation not related to Hamiltonian operators, 
the property is a consequence of the Leibnitz rule (1501) and of the Jacobi identity for 



Proposition 11 (|35j). The anchors a map the brackets [, ] between sections of bun- 
dles Q over finite-dimensional manifolds M n to the Lie bracket [ , ] on the tangent 
bundle to M n . 
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The converse is not true. Namely, there are morfisms of Lie algebra structures in 
the modules of sections of fibre bundles over the base manifold that do not respect the 
Leibnitz rule (*5"U1) . Moreover, and this is our claim of crucial importance, the converse 
of Proposition [TT] is not valid for Hamiltonian operators P in total derivatives and the 
morphisms a = £pq: ^(tt) — > q(tt), where Dorfman's bracket [, ]p on f2 1 (7r) is ( 145 aft . 
Indeed, there is no Leibnitz rule over J^(vr) for [, }p. This is the obstruction to an 
extension of Dorfman's bracket [, \p of 1-forms ipi-,ip2 G T(T*M n ), 

[i>l,i>2]v = Lp^i^2 - Lp^-01 + d d R(Af") (v(i>iM), Per(/\ 2 TM"), (53) 
to Lie algebroid structures over the base manifolds J°°(tt). 

Counterexample 4. Consider the first Hamiltonian structure P\ = D x for the KdV 
equation ( fTTl) . and consider two conserved densities p\ — w and p2 = |w 2 . Obviously, 
{pi dx, p2 dx}p 1 = 0, and hence Dorfman's bracket [VtjV^Jfi of the gradients ipi = 
E w (pj) also vanishes. Now multiply p\ by any nonlinear f(x) G C°°(M) and, applying 
Euler derivative E^, get / • ipi. First let us commute / • ipi with ip 2 by [ , ]p 1 and obtain 
the equivalence class in the quotient x(7r)/kerPi, 

[f(x)-l,w] Dx = f'(x) + const. (54) 
On the other hand, the Leibnitz rule ("5*0"]) prescribes that this is equal to 

f(x) • [1, w] Dx - 8. Dx{w) (/(a;)) -1 = 0-0 = 0. 

This zero value can be achieved at any point Xq G M by choosing const := —f'{xo) 
in (I54p . but not on the entire M at once. 

Analogously, the Leibnitz rule does not hold for multiplication by f(x, [u]) G J 7 {ft), 
and we do not repeat the reasonings for the sake of brevity only. 

This counterexample manifests the fundamental difference between jet bundles J°°(7r) 
over M n and smooth manifolds themselves, which can be regarded as the jet bundles 
over the point {x }, n = 0. 

We conclude that, depending on a problem, we have the choice between, first, as- 
suming that the Lie algebra homomorphism A: fi 1 (^ 7r ) — > is enough for further 
reasonings, and, second, postulating the Leibnitz rule (15*01 . The former case (14*41) is 
realized by Frobenius operators. In this sense, they can be regarded as the anchors in 
the analogues of Lie algebroids over infinite jet bundles J°°(7r). 

However, if the Leibnitz rule is still needed, then we propose to resolve the difficulty as 
follows. First we assign formal differential complexes ( J72|) to Frobenius operators, and 
the representations of the complexes through homological vector fields will determine 
the Lie algebroids (see (jBTl) in the second part of definition's]). This scheme can be 
further applied in the Batalin-Vilkovisky (BV) formalism to construction and analysis 
of the quantum Poisson manifolds [TJ [3] . 

2.4. The bracket of conservation laws. Let us convert Lemma [9] to a definition. 
Suppose £ = {F — | F: J°°(n) — > f} is a determined evolutionary system, hence 
f(£) ~ and the coupling ( , ) : f x x — > A n (7r) is well defined. In this section we 

consider the class of Noether (or 'pre-Hamiltonian' [22, 24J) operators A: cosym£ — > 
sym£ that map generating section of conservation laws for £ to symmetries but may 
not be skew-adjoint or, even if it is so, may not define a bracket that satisfies the Jacobi 
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identity flHUl) . In all cases, if a Noether operator for an autonomous evolution equation 
£ = {F = u t — f(x, [u]) = 0} does not depend on the time t, then A satisfies the 
operator equation £ F o A = A o £ F on £. 

Using Noether operators, one can induce skew-symmetric brackets between conser- 
vation laws r] G H n ^ 1 {£) and, second, between their generating sections f(£). 

Let 7] = pdx + • • • G A n_1 (£^) be a conserved current for an evolutionary system £. 
Then, by Lemma (U its generating section ip v = E(p) G f(£) is the Euler derivative of 
the conserved density. For any ip G sym£, the current £^(77) is obviously conserved 
on £ and its generating section E(£ (/ ,(p)) equals E^ip) + £t{ip) by Lemma El 

Given an operator A: cosym£ — > sym£ and any two conserved currents 771,772 £ 
A"- 1 ^), setE! 

((, )) A : A n -\£) x A n -\£) - A"- 1 ^), (55a) 

((Vi,V2))a = - £a^ 2 )(^i) + K^i, M^lh)) - I (^1 ), $2), 

and (see (JSSD, which is [35, Eq. (3.2')]) 

lj A :}(£)x}(£)^}(£). (55b) 

UiMa = £a Wi) (^) - £^ 2 )(V0 + + !E«tfcL,i4(V*)) - Wi),V>2», 

By (BID , we have {T^T^Ia = — £a(e(Hi))(7^2) for Hamiltonian operators A, and there- 
fore ((7Yi, H2)) a = — {7^1,7^2}^- Hence we conclude that formula (155bj) gives an alter- 
native way to calculate the bracket {{ , }}a, see f)45ap and (|4"5"j) . Indeed, this is achieved 
by omitting the first two summands in (I55bl) and setting to zero the application of 
evolutionary derivatives to ipi and ^2 in what remains. The equivalence between (|46|) 
and ( I55bl) implies nontrivial identities even for the simplest scalar Hamiltonian oper- 
ator D x . Formula ( H6l) is more elegant because it contains only the action by total 
derivatives onto ipi,^ G f(£) (note that the partial derivatives djdu % a are applied in 
it only to the coefficients of A). At the same time, the last four summands in ( 155bl) 
incorporate an application of the partial derivatives to ip\ and ip2- However, the skew- 
symmetry {{ipi, ^2}} a = ~{{4'2,i , i}}A of the expression (pTOl) is not obvious. Again, we 
recall that the notation £ A ,^, see is not the linearization £a{iP) given in (13TJ]) . 

Remark 2. Only skew-adjoint operators T*M n — > TM n , which can be represented 
by bi- vectors r(/\ 2 TM n ), were considered in [35l §3.2]. This was why the antisym- 
metrization was not performed in (1531 . which is the realization of (155bl) in the case of 
finite-dimensional manifolds M n . 

Due to the antisymmetrization, brackets (1551) are invariant under adding self-adjoint 
operators to A. This shows that the brackets (( , ))a and [ , }a satisfy the Jacobi identity 
if and only if A = P + S, where P is a Hamiltonian operator and S* = S. 

Example 5. Consider the dispersionless 3-component Boussinesq-type system [18] 

u t = ww x + v x , v t = —uw x — 3u x w, w t = u x , (56) 

which is equivalent to the Benney-Lax equation. Note that (1561) is scaling-invariant; 
e.g., the homogeneity weights are \u\ = 3/2, \v\ = 2, |w| = 1, \D X \ = 1, and \D t \ = 3/2. 

llr The arising algebra (H 11-1 (£) , [ , ] a) of conservation laws converts the Lie algebra of Hamiltonians 
with a Poisson bracket to an algebra with the bracket specified by a Noether operator A : cosym £ — > 
sym£ . This algebra could be called the "current algebra" if the term were not already in use. 
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In [21], two compatible Hamiltonian operators Aq and A 2 , and a self-adjoint Noether 
operator A\ : cosym £ — > sym £ were found for ([56]) by a calculation on the £*-covering 
over £ and using the scaling invariance of the system. The first Hamiltonian structure 
for (|56|) is given by the operator 



D 



X 




A = ( -AwD x - 2w x D x I . (57a) 

Sokolov's bracket {{ , }}^ o in the inverse image of A is obtained using formula ( H6l) . 
and we have that 

{{ , Hi = « , }}% = 0, {{p, q}}\ = 2 (p v q v x ~ PW) • (57b) 

The self-adjoint Noether operator A\ is 

ww x + v x —3u x w — uw x u x \ 

A\ = J —3u x w — uw x —3w 2 w x — Av x w — uu x v x . (57c) 
u x v x w x J 

By Remark [21 the brackets (( , ))a ± on H x {£) and [, on k{£) vanish identically. The 
second Hamiltonian operator for ( 156]) is the linear combination 

A 2 = A 1 + A 2 , (57d) 

where, following [21] , the operator A 2 : cosym £ — > sym £ is equal to 

((2w 2 + Av) D x + ww x + v x — lluw D x — (2u x w + 8uw x ) 3u D x \ 
— lluw D x — 3u x w — uw x h\ D x + h Av D x , (57e) 

3u D x + u x Av D x + 2v x 2w D X J 

here we put h = — {2uu x + 8vw x + Av x w + 6w 2 w x ) and hi = — (3u 2 + 16vw + 6w 3 ). The 
right-hand side of system (1561) belongs to the image of each of the three operators Aq, 
Ai, and A 2 : 

'u t \ (\w 2 + i\ fl\ /f 

v t = i u = A x = A 2 

W V ™ / V / V°, 

Again, by Remark [21 the bracket [, \a 2 coincides with Dorfman's bracket [, ]^ for 
the second Hamiltnian operator, and both of them satisfy the Jacobi identity. Thus 
Sokolov's bracket {{ , J}^ can be caclulated not by ( fl6l) . but using ( I55bl) and omitting 
all terms with evolutionary derivations. 

In what follows, the Noether operators Ao, Ai, and A 2 : cosym. £ — > sym£ for ( J56l) 
will be met again. We claim that the images of Ai and A 2 are closed w.r.t. the commu- 
tation ( 1291) . The images of arbitrary linear combinations AoA) + A1A1 + A2A2 are also 
closed under commutation, and hence the three Frobenius operators are linear com- 
patible. Moreover, the image of a pre-symplectic operator A^ 1 is closed w.r.t. either 
[,]_!„ or [, ]_£ on x(£), and therefore both recursion operators i?o = A) ^i" 1 an d 
R 2 = A 2 o A^ 1 : sym£ — > sym£ are Frobenius as well. 
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In the previous subsection we addressed the construction of Lie algebroids over J°°(7r) 
and noticed a very important property of (|55p : For operators A = P + S in total 
derivatives, the Leibnitz rule is valid for the bracket (j55aj) but does not hold for (155bj) . 

Not every skew-adjoint Noether operator is Hamiltonian since, in general, the brack- 
ets ( 1371) and ( 155bl) may not satisfy the Jacobi identity. 

Lemma 12 ([23]). Let £ = {u t = f(x, [u])} be an evolutionary system of differential 
order ord£/ > 1, and let its symbol be nondegenerate on an open dense subset of £. 
Then any skew-adjoint Noether operator A : cosym £ — > sym £ is Hamiltonian. 

Unfortunately, the non-vanishing of the determinant of the symbol (which is the 
higher-order matrix Yl\a\=max(^f)v) depends on the system of local coordinates. 

Counterexample 6. Consider the Kaup-Boussinesq system 

u t = uu x + v x , v t = (uv) x + e 2 u xxx , (58) 

whose symbol is degenerate and which extends the assertion of Lemma [12], see [23] . 
System ( 1581) is transformed by w = v — eu x to the second-order Kaup-Broer system 

u t = uu x + w x + eu XX} w t = (uw) x - ew xx . 

After the transformation, the determinant of the symbol does not vanish anywhere (of 
course, we assume e ^ 0). At the same time, it is clear that the third-order linearization 
of the right-hand side of the Kaup-Boussinesq system (15*81 is nondegenerate in the 
following sense. 

Definition 3. We say that a matrix operator A in total derivatives is nondegenerate if 

P|kerA a = {0}, where A = ^ A a ■ D a . (59) 

a cr 

We shall consider nondegenerate Noether and recursion operators that are transformed 
using (|3"T1) and (IHHj) . respectively, under nondegenerate reparametrizations u — u[u) of 
fibre coordinates. This makes the nondegeneracy (15*9"]) well defined. 

Example 7. The Hamiltonian operator ) is degenerate. The operator □ = u x + 
\D X , which was introduced in ([H]) and which will be made well defined by (fijoT) below, 
is nondegenerate, although its kernel is spanned by exp(— 2u). 

A conjecture that the assumptions of Lemma [12] can be weakened is already present 
in [23] , but its proof [15] , which is quoted there, hinted no way how it can be generalized. 

Conjecture 13. The assertion of Lemma [T2l holds for skew-adjoint Noether operators 
A: cosym £ — > sym£ on evolutionary systems £ = {u t = /} with ordif > 1, if the 
linearizations if of the right-hand sides are nondegenerate, see ("5"9"]) . on open dense 
subsets of £. 

3. FROBENIUS OPERATORS 

In this section we introduce the well-defined notion of Frobenius operators, reveal- 
ing the nature of their domains ^ 1 (c; 7r ) C r(7r^ D (^)), establish their properties, and 
give examples of Frobenius operators. In particular, we recognize Frobenius recursions 
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x(7r) — > x(7r) as higher differential order solutions of the classical Yang-Baxter equa- 
tion for the Lie algebra g(vr); here we construct examples of the projector solutions 
which are related to Liouville-type systems (e.g., to the 2D Toda chains). We estab- 
lish the chain rule for factorizations of Frobenius operators and correlate the Sokolov 
brackets for Frobenius recursions that produce sequences of Lie algebra structures. We 
show that Frobenius operators A induce flat connections in triples (fi 1 (£„-), g(vr), A) 
of Lie algebras and their morphisms. Finally, we assign formal differential complexes 
over fi 1 ^) to Frobenius operators. 

The Frobenius theorem [8], which deals with involutive distributions over finite-di- 
mensional manifolds, gave the proper name to the class of operators defined here. 

3.1. Main definitions. By Proposition El the Poisson bracket { , }a in the Lie al- 
gebra of Hamiltonians H n (ir) is correlated with the commutator [, ] in the Lie alge- 
bra q(tt) of vector fields £a( )- The skew-adjoint operator A in the Poisson bracket 
{H\,H2}a = (E(7ii), A(E(H.2))) coincides with the operator A whose image is a sub- 
algebra of evolutionary fields £a(e(W))- Thus we have 

[£a(E(Wi)), £a(E(W 2 ))] = £a(E({Wi,« 2 }a))- 

This is readily seen from the Jacobi identity ( 1481) for the Schouten bracket ( 1351) . Indeed, 
the commutator of two Hamiltonian vector fields equals 

[[A,Hi], [A, H 2 \] = [A, [H 2 , [A WJH + [W2, [[A,Ki],A]], 

whence the derived bracket [7i 2 , [A Will = {^1,^2} a appears in the right-hand side 
of (IBUjl and d\(Hi) = is contained in the second summand. 

In this section we reverse the status of the two Lie algebra structures [ , ] . A and { , }a, 
giving the priority to the commutators of vector fields and thus considering involutive 
distributions of evolutionary vector fields 6^ G g(7r) whose generators if = A(-) G x(tt) 
belong to images of matrix operators A in total derivatives. 

Definition 4. The definition of a Frobenius operator consists of two parts; all notation 
is correlated with the previous material. 

§1 (see [60J). Prohibit any changes of coordinates in all modules of sections. Con- 
sider a matrix linear differential operator A: f — ► x(tv) in total derivatives, where 
f Q r(7r^ (^)) is a certain sub module that will be specified in §2 of this definition. The 
operator A is Frobenius if, first, its image in the ^ r (7r)-module x(n) of generators of 
evolutionary vector fields is closed w.r.t. the commutation: 

[imA, imA] C imA. (61) 

The commutator [, ]|. . induces^! a skew-symmetric Koszul bracket [, ]a in the quo- 
tient f2 1 (^ 7r ) = f/kerA of the domains of A: 

[A^), AM] = A([^,^}a), Vi,^ e n 1 ^). (62a) 

12 Two sets of summands appear in the bracket of evolutionary vector fields A(ipi), A(ip2) that belong 
to the image of a Frobenius operator A: 

[Afa), Afo)] = A(Z A{ ^ 2 ) - E A(lh) {1>i)) + (E MM (A){fM - E A (^)(A)^i)). 

In the first summand we have used the permutability of evolutionary derivations and operators in total 
derivatives. The second summand hits the image of A by construction. 
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The Koszul bracket, which is defined up to ker A, equals 

[ipi,ip2U = £a(^)(^2) - £a(^ 2 )(^i) + {{^1^2}} a- (62b) 

It contains two standard summands and the Sokolov bracket {{ , }}a- 

§2. Let id be a fibre coordinate in the bundle £. Allow coordinate reparametriza- 
tions, both in £ and tt. Suppose that there is a differential substitution J°°(7r) — > T(£) 
that yields the embedding : f <^-> r(7r^ c (£)) of an jF(£)-module f. From now on, re- 
strict f onto the image of the substitution. Thence we denote the substitution by the 
same letter w[u] and retain the notation f for the restriction of the ^-"(£)-module. Put 
fl 1 (^) = f| M /kerA 

/ 1 \ w=W[u\' 

We postulate that, under diffeomorphisms u = u[u}: J°°{ti) — > T(tt) and w = 
w[w] : </°°(0 — > r(£), the transformation rules for the operator A are uniquely defined 
by the fibre bundles 7r and £. This implies that f is one of the following^ 
f = x(£) , Frobenius operators of first kind are 

' \^ / 1 W=W[U\ r 

A: x(OL-x(tt). (63) 

Under any diffeomorphisms u = u[u] : </°°(7r) — > T(tt) and w = w[w] : </°°(0 — > 
T(£), these operators are transformed according to 



A ^ A = t 1 ^ o A o 



MO I w- w [u] ' Frobenius operators of second kind are linear mappings 



(64) 



A: x(OL^x(tt). (65) 

For any differential changes of coordinates u = u[u] and w = w[w], these oper- 
ators obey 

A ~JU#Wo (#>)*,. (66) 



u=u\u\ 



Sections ip that constitute the domains ^ x (£^-) of Frobenius operators (I6"31 of first kind 
are transformed as vectors = if) • d/dw + • • • ; in the case (1531) . the sections ^ satisfy 
the rules valid for variational covectors if> = S(-)/Sw. 

If Frobenius operator A is a recursion x(7r) — > x(7r), then it; = id, and the trans- 
formations of the domain and image of A are uniquely correlated. If A is a Noether 
operator x(7r) — ► x(7r), then 10 determines the system £ = {to [it] = 0} of differential 
equations, and, by default, we use the misleading isomorphism f ~ x(7r) for evolution 
equations (see Remark [1] on p. [T5l) . Otherwise, there may be no constraint between 



13 The axiom we accept means that, in this paper, we consider only the tangent and cotangent 
bundles to (infinite-dimensional) smooth manifolds. Therefore we do not attempt to define Frobenius 
operators on .F(£)-modules other than listed. Second, we do not consider operators ® i=1 >f(£) C5> 
1 r 1 — > xM whenever k + £ > 1. This would generalize 52 of this definition to the same 

— i \^/\ w — w[u\ v ' 

extent as the Schouten bracket [ , ] incorporates the commutator [ , ] and the Gcrstenhaber bracket 
[ , \p contains the Koszul bracket [ , ]-p. However, the full generality must not be discarded (e.g., when 
the isomorphism f ~ >c(tt) is not used for Noether operators A, see Remark [TJ. Finally, note that, 
when allowing nonlinear differential rcperametrizations F — F[F] £ f of equations, we arrive at the 
theory of nonlinear Frobenius Noether operators. 
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transformations of the domains and images. In particular, the gauge group of f can be 
trivial, meaning that transformations of w in £ are still prohibited, as in §1: 

f = r(7r^ (£)): These Frobenius operators A: r(7r^(£)) — > x(n) are transformed 

by A i— > A = o A\ r -, under u = u\u\. 

J u \U=U[u\ L J 

Note that no substitution w[u] is needed in the third (degenerate) case. 

Remark 3 (On the "Frobenius theorem"). Frobenius operators A: f — > x(7r) specify 
involutive distributions of vertical symmetries £>a(-) £ D v (J°°(7r)) of the Cartan distri- 
bution C = (d/dx l ) C D(J°°(7r)), see ( 1271) . which is itself Frobenius. 

Assume that for such an involutive distribution (£a(-)) there is an integral manifold 
(typically, it would be infinite-dimensional), and suppose further that it is a differential 
equation S. (Note that by an equation we mean the infinite prolongation, which does 
not always exist.) The definition of Frobenius operators implies that £ admits infinitely 
many symmetries ip = A(<f) which contain free functional parameters (f) 6 f. This prop- 
erty is close but not equivalent to Definition [H] of Liouville-type differential equations 
(see p. W7\ and Remark [T71 that follows). 

Lemma 14. The bracket [, ]a is IR-bilinear, skew-symmetric, and transfers the Jacobi 
identitj0 from the Lie algebra g(7r) of evolutionary vector fields to fi 1 (£ 7r ), see flBT|) . 
The kernel of A is an ideal in f . The image of a Frobenius operator may not be an ideal 
in the space of evolutionary derivations although it is a Lie subalgebra by definition. 

Remark 4. The bi-differential skew-symmetric bracket {{ , }}a £ CDiff (/\ 2 fi 1 (£ 7r ), fi 1 ^)) 
does not generally satisfy the Jacobi identity. Indeed, for the Koszul bracket [, }a we 
have 

= 5^[[^i> V'aU V'a]^ = 5Z[ £a Wi)(^2) - £a(^ 2 )(^i) + {{tpi,tp2}}A,^3] A 

o o 

+ E a( {{ ^ } u) - £ ^)({{^i^ 2 }}a) + {{{{AMU, ^}}a}. (67) 

The underlined summand with a derivation of the coefficients of {{ipi, A}} a, which 
belong to ^(vr), may not vanish on the ^-"(£)-module f, when the action of £a(-) on t° the 
basis of sections ipi is set to zero (see Remark [201 on p. [621) . Note that the Jacobi identity 
for [, ]a then amounts to the last line of ( |671) . which equals ^j ) [{{'0i,'V ; 2}}yi,'?/ ; 3] /l = 
and which is a half of the equation for infinitesimal deformations of [, ]a- 

Remark 5. In this paper, we study only those operators whose images in x(tt) are 
closed w.r.t. the standard bracket ( |29l) . Suppose, however, that an endomorphism 
R G End]R x(7r) is Frobenius and hence induces a new Lie algebra structure [ , ]r on the 
quotient x(n) / ker R of its domain. Contrary to Proposition [HI the images of standard 
Hamiltonian operators are generally not closed w.r.t. the Koszul brackets [, ]r induced 
by Frobenius recursions R. 



In what follows, the nondegeneracy condition ([59]) specifies the case when the Jacobiator 
Jtyi, i>2, V's) f° r [i }a vanishes identically, not being a nonzero element of ker A as in [50] , 
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Although one could repeat the whole construction for operators — > x(7r) whose 
images are closed w.r.t. [, } R , this does not produce a new formalism. Indeed, the 
operators A, whose images in are closed w.r.t. [, ]r, generate standard Frobenius 
operators R o A, and vice versa. At the same time, the new formalism appears for 
operators A whose images in x(tt) are closed w.r.t. the r-brackets ( 1761) given by r- 
matrices for the Lie algebra fl(7r). 

Remark 6. Specifying the Koszul bracket [, ]a on fi 1 ^) by a Frobenius operator A, 
we define Frobenius generalizations uj: x(7r) — > Q 1 ^) of symplectic structures >c(tt) — > 
x(7r). To this end, we require that the images of uj are closed w.r.t. [, ]a- Since all 
modules are already known, the adaptation of §2 of Definition H] for uj is obvious. 

3.1.1. Flat connection. Let A: f — ► >c(tt) be a Frobenius operator. It provides the 
homomorphism of Lie algebras 

A: (^(UUa) ->(*(*),[,]). (68) 

Let .fT be a g(7r)-module; typically, K = J-(ir) or any other horizontal JF(7r) -module 
(e.g., q(jt) itself). By the homomorphism A, the 0(7r)-module K is an Jl 1 (^-module as 
well. The Jacobi identity implies that the adjoint action by an element of Lie algebra 
fi 1 (^ 7r ) is a derivation. We further impose the condition of semi-simplicity: we claim 
that 

bl>i, -)a = [fa, -)a ipi = ifa, 

for equivalence classes ipi,^2 £ ^ 1 (Ctt) — f/kerA 
Now we define a connection V A in the triple (1681) . 

V A : Der Int (fi 1 (^),g(7r)) ^ Der (q(ti) , K) . 

The construction is analogous to the definition of connections in triples (A, B, i), where 
A is a k-algebra, B is an Aalgebra, and %: A — > i3 is an algebra homomorphism [38J. 
In its turn, this generalizes the connections in fibre bundles E n+m — > M n , when .4 = 
C°°(M n ), = C°°(£ n+m ), and z: A B. We set 

V A : Ao[^.]^[i(^),-]. (69) 

The above definition is f2 1 (^ vr )-linear. Indeed, for a derivation A = [ip, -}a we have that 

V? xA = ^(/)xVl, /Gft 1 ^), AeDerw^y, (70) 

where the multiplication x by / and by its image under A is the adjoint action. Note 
that the right-hand side of the analogue of (170!) in a classical definition of C°°(M n )- 
linearity of connections in fibre bundles 7r does contain the image i(f) of the identical 
embedding A B and not f E A itself. 

Proposition 15. Connection ( 1691) is flat: 

(V£ o - o V£ - V£ dA ) fa) = 0, Vp.gGfi 1 ^), V?e (7r). (71) 

Proof. The Jacobi identity for the bracket (I29p of generating sections of evolutionary 
vector fields, 

[A(p), [A(q),<p]] + [A(q), [<p,A(p)]\ + [<p, [A(p),A(q)]] = 
is the flatness condition (ITTj) . □ 
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Corollary 16. If the flows of a commutative hierarchy 21 belong to the image of a 
Frobenius operator A, then the hierarchy is a geodesic w.r.t. connection (1ST?]) . 

Proof. Indeed, for any curve ip(r): R -> fi 1 ^) located in the inverse image of 21 un- 
der A, the covariant derivative V ^ T ^A{ip' {r)) of the velocity ip'(r) vanishes along the 
curve. □ 

Remark 7. The Frobenius operators (1681) induce a morphism H 1 (^ 7r ) — > /\*g(vr) to the 
Schouten algebra of evolutionary polyvector fields, which is endowed with the Schouten 
bracket ( J35l) . Flat connections ("69l) in triples ( |68l) are naturally extended to connections 
in (O 1 (£„.), A* fK 71 ")' which are flat in the graded sense. 

Remark 8. The notion of connections for the triples (which generalize A = C°°(M n ), 
B = C co (E n+m ) 1 and 1: A B) was proposed by Krasil'shchik in [3B]. It leads to the 
Cartan cohomology, see [lOj Ch. 5], which allows to interpret symmetries and recursion 
operators as equivalence classes. Therefore we expect to encounter some cohomology 
for the connection ( |69i) in the triples (l6~8"l) . 

Of course, connection (l6~9l) in the triples ( 1681) is not the Cartan connection ( 1271) 
on J°°(7r). This is readily seen from the fact that everything at hand is it- vertical and 
is projected to zero vector fields on TM n under ir^*. (Everything would be projected to 
the point Xq if the jet bundle were a finite-dimensional manifold N m for it: N m — > {xq}, 
see (|54l) .) Therefore, instead of a dual description of the n-dimensional distribution 
C C D(J°°(7r)) by Cartan 1-forms, we face the problem of dual representation of the 
involutive distributions (£a(-))- 

3.1.2. Frobenius complex. Now, using Cartan's formula, we construct a differential com- 
plex on the Chevalley cohomology HomR(/\ fc fi 1 (£„-), q(tt)) with values in an fi 1 ^)- 
module q(tt). This is the Frobenius complex assigned to an operator A: $V-(£, n ) — > fl(7r). 

For any k > and u>k E HoniR(/\ fe fi 1 (£?,-), q(tt)) , define the differential d: ujk 1— > du^ 
by setting 

d^fe^o,---,^) = ^(-l)*V$.u;(^o,---,?i,---,^Jfe) 

i 

+ ^(-l) i+ ^ 1 ^fc([^,^-]A,^0, ...,$i,...,ll>j,...,ll>k)- 

i<j 

Hence we obtain the complex 

0(7r ) ^ Hom M (fi 1 (e 7r ),g(7r)) Hom M (/\ 2 fi 1 (e 7r ),s(7r)) 

i>Hom R (/\ 3 fi 1 (e 7r ),0(7r))^---. (72) 

The first inclusion in (1721) consists of the commutations VA</?o with fixed elements 
V^o G 0(vr), the second arrow is the composition of the Koszul bracket [ , ]a with A, and 
the third arrow calculates the right-hand side of the Jacobi identity fIBTj) . The Poisson 
complex is a special case of (172*1) . extending it from the left with the ^ r (7r)-module of 
n-th horizontal cohomology H n (ir) for J°°(7r). Again, we assume the semi-simplicity: 
no elements (p G g(7r) commute with the entire Lie algebra q(tt). 
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The home-logical vector field representations (e.g., see (|32p ) for the differential d 
allow to encode the Lie algebroids r(£„.) — > k{ti) over infinite jet spaces J°°(ir) in terms 

of flED). 

3.2. Examples of Frobenius operators. Throughout this paper, we consider only 
the commutative case when the base manifolds M n are not Z,2-graded and all coordi- 
nates on the jet bundles J°°(7r) are permutable. Also, we confine ourselves to the local 
setting and consider Frobenius matrix differential operators which are polynomial in 
total derivatives. However, the search for nonlocal Frobenius (super-)operators can be 
performed using standard techniques, see Remark [201 on p. ______ 

Let us have some examples of Frobenius operators of second kind and of Sokolov's 
brackets on their domains. We start with Hamiltonian operators, which are the most 
well studied examples of Frobenius structures. 

Example 8. Every Hamiltonian operator is Frobenius. Indeed, the criterion in Propo- 
sition [6] gives formula ( H6T) for the bracket {{ , }}a on the domains k{n) of nondegenerate 
Hamiltonian operators A, which obey the law ( 15B1) under a change of coordinates. This 
demonstrates that the formalism of Frobenius operators is a true generalization of the 
Hamiltonian approach to nonlinear evolutionary PDE. The use of Frobenius operators 
can be a helpful intermediate step in the search for and classification of the Hamiltonian 
structures. An advantage of this approach is that it is easier to solve first equation (1U21 
w.r.t. operators A and filter out skew-adjoint solutions rather than to solve the Ja- 
cobi identity in the form of fHUj) . Thus, by Proposition [HI a nondegenerate skew-adjoint 
Frobenius operator A is Hamiltonian iff the Sokolov bracket {{ , }}a is equal to the r.h.s. 
of ( fl6l) up to ker A. Also, we note that we do not require a restriction of the operators 
A onto differential equations £ such that A: cosym£ — > sym£, unlike in [22] . 

Remark 9. A recent version [____l[__5] of the Darboux theorem for (1 + 1)D evolutionary 
systems implies that Hamiltonian operators A for non-exceptional systems can be trans- 
formed to const ■ D x . The bracket {{ , }}a is then zero, which is readily seen from f|46|) . 
Hence the actual statement of the Darboux theorem for PDE is that the bracket {{ , }}a 
can be trivialized for such Hamiltonian operators. 

Remark 10. One bracket {{ , }}a can correspond to several operators A that satisfy (loTj) . 
For example, the second structure A^ dY = — ^Dl+2wD x +w x for the KdV equation (fTTT) 
determines the bracket {{p, gll^Kdv = pq x — p x q, which is also induced by the opera- 

_____ 2 /Q\ 

tors (1751) . see Example fTUl below, and by A\ = D x o( y D x + u), see Example fT51 Actually, 
this Wronskian-based bracket (c.f. [27]) is scattered through the text, see Example _____ 
on p. __6j Hence there are fewer brackets {{ , }}a than there are Frobenius operators A. 

Formula (H6l) does not remain valid for arbitrary Frobenius operators, which are 
generally non-skew-adjoint. In appendix _B] we describe an inductive procedure that 
reconstructs the bracket {{•,•}} a for nondegenerate operators A, see (loT?]) . From now 
on, we consider only nondegenerate Frobenius operators. 

Example 9. Noether operators (loT|) for the 3-component dispersionless Boussinesq- 
type system ( !56l) are Frobenius. For the Frobenius operators Ai,A2~. cosym£ — > sym£, 
the transformation formulas ( 1581) classify them to the second kind; the substitutions 
J°°(tt) — > are the identity mappings. The same is obvious for Hamiltonian op- 
erators A and A 2 for this system. However, unlike [, }a 1 =0 and [, ]a 2 — [, li 2 ' 
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see (]55bp and Remark [21 the components of Sokolov's brackets for A\ and A 2 are ob- 
tained using the inductive algorithm, and the result is (I128p on p. Bi-differential 
representations for components of these brackets are formulated in appendix O 

In section we shall consider a class of Frobenius operators □ : cosym E\ — > sym £ 2 of 
second kind that map (co-)tangent bundles for two hierachies of evolutionary systems 
related by the substitution w: E 2 — > 

Now we list several operators which are known in a fixed system of local coordinates 
and whose images are then closed under commutation. The transformation rules for 
their domains are unknown, although it may occur that these reparametrizations are 
uniquely determined by a change of coordinates in the images. 

Example 10 (KdV scaling weights). Let us fix the weights \u\ = 2, \D X \ = 1 that 
originate from the scaling invariance of the KdV equation ut = —\u xxx + 3uu x , see 
also (II ip ; we have that \D t \ = 3. Using the method of undetermined coefficients, 
we performed the search for scalar Frobenius operators that are homogeneous w.r.t. 
the weights not greater than 7. We obtained two compatible Hamiltonian operators 
Af dV = D x and Af dy = -\D\ + 2uD x + u x , the generalizations D 2 x n+1 of D x , and the 
operator 

uuxxx + 3uu xx D x + 3uu x D 2 X + u 2 D x . 
Also, there are four non- skew-adjoint operators that satisfy (j6ip . 

Af =u 3 -u 2 x , Af = 2u 2 x - uu xx - 2uu x D x + u 2 D 2 , 

A { p =u 2 x D x - 2uu xx D x - Auu x D 2 x - Au 2 D\, {{p, q}} A m = u 2 ■ (pq x - p x q); 

Ag ] = -2u x u xx - u 2 x D x , {{p, q}} A (7) = 8u xx ■ (pq x - p x q) + 2u x ■ (pq xx - p xx q). 

Finally, we have found the operators that contain arbitrahry functions: f(u)D x and 
f(u)u 2 with vanishing brackets, and also we have got 

A 3 = f(u)u x , {{p, q}} Ai = f(u) (p x q - pq x ) ; 

A 4 = f(u) (u D x - u x ) , {{p, q}} Ai = f(u) (pq x - p x q) . (73) 

Formula ( |46l) is not valid for any of these non-skew-adjoint operators. 

Example 11. There is a three-parametric family of scalar first-order Frobenius oper- 
ators A = a{u,u x ) D x + b{u,u x ) with brackets {{p, q}} A = c(u, u x ) ■ (pq x — p x q). They 
are given by 

a = y(u), b = (3{u) -u x + a ■ y(u), c=-(3(u), (74) 
where a is a constant and the functions (3, 7 are arbitrary. 

3.3. Frobenius recursion operators. Within this section, we set £ := it and w — u. 
Thus we consider Frobenius recursion operators R: x(-k) — > x(7r) whose images are 
closed w.r.t. (|29l) . The standard Lie algebra structure [, ] is transferred to [,]r on the 
quotients 1 (7r 7r ) = / ker R by the classical Yang-Baxter equation for g(ir), 

[R<pi,R(p 2 ] = R([(fi, <p 2 )r) , <fi,<f2 e x(vr)/keri?. (75) 
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For R id, the new Lie bracket [,}r is different from the original commutation of 
evolutionary vector fields in the image. 

Remark 11. Frobenius recursion operators provide higher differential order solutions of 
the classical Yang-Baxter equation for the algebra g(7r). The parallel with the zero- 
order theory of r-matrices is as follows. Let g be a Lie algebra with the Lie bracket [ , ] . 
The classical r-matrix [U [2JJ E2] is a linear map r : g — * g that endows g/ ker r with the 
second Lie product [ , ] r ; for any a, b G g set 

[a, b] r = [ra, b] + [a, rb}. (76) 

A sufficient condition for an operator r of differential order zero to be an r-matrix is 
that r satisfies the Yang-Baxter equation YB(a): 

[ra, rb] — r([a, o] r ) — a[a, b], a,6Gg, a = or 1. (1751 ) 

Now let a recursion differential operator R G EndjR x(ir) be Frobenius, then the Koszul 
bracket ( 162bl) satisfies the Yang-Baxter equation ( 1751) with a = 0. 

We conclude that, for Lie algebra q(it) of evolutionary vector fields on jet spaces J°°(7r), 
the classical Yang-Baxter equation ( 1751) admits solutions [, }r) of form (162bl) other 
than the standard r-brackets ( 1761) . The two notions of r-matrices and Frobenius re- 
cursions R G EndiR x{n) as solutions of equation (1751) are very close. However, the 
distinction between them is expressed by the respective brackets [ , ] r and [ , ]r on their 
domain the Koszul bracket (162bj) is calculated and the r-bracket [, ] r is postu- 

lated. This explains why §1 of Definition H] incorporates the form of the Koszul bracket. 

Example 12 (Projector solutions). Instead of the entire jet space J°°(7r), consider a 
Darboux-integrable Liouville-type hyperbolic system [56, 60J 

£l = {< y = f(x, y; [«]), ...,<,= f m (x, y; [«])} c 

see Definition[8]on p.HTJand a comment on it in Remark [TBI Essentially, £l — {J°°(tTx)® 
J°°(7ry) j ~, plus the rules for calculating u x .., X y...y\ , where ir x , % y : R m xR — > M. are trivial 
bundles over R with the base coordinates x and y, respectively, and ~ glues together 
their fibres with coordinates (it 1 , . . . , u m ). 

In this notation, the symmetry algebra sym^L is decomposed along the x- and y- 
characteristics such that any ip G sym£ L is of the form (11001) . see p. SHI 

where sections G x(^) | and G x(^) | belong to certain restrictions of the modules 
and the operators □ and □ are Frobenius. This is discussed in section [5] in detail. 

Without loss of generality (we consider the x-components of all structures, but one 
could replace x <-> y in what follows), set 

R - £ nw+Dw ^ e sym^ L . (77) 

The Liouville-type systems are so special that R is a recursion on sym^L, although 
keri? is very big under this projection. The image of R is closed w.r.t. the restriction 
of bracket ( !29l) onto £l- Thus we obtain an example of a Frobenius recursion operator. 

Surprisingly, we can apply here the algebraic construction of Poisson brackets, which 
was introduced in [UJ, to such projector recursions R. We recall that a homomorphism 
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I! G Hom a (V, End V) is Hamiltonian for a commutative k-algebra o and an a-module V 
if the Poisson bracket {vi,v 2 }u '■= (n(t>i)) (t^) of vi,v 2 G V satisfies the commutation 
closure condition U({vi,v 2 }n) = [n(t>i), n(u 2 )l . The latter reduces to (!60l) if we choose 
k = E, a = F(ir), V = H n (-7r), and n = £a(b(-))- 

In our case, we have a = M. and V = sym^L- Define II G Hohik( sym^L, EndjR sym^) 

by 

U((f) = &d R{(p) oR, ip G sym£ L , 
whence we obtain {ipi,ip 2 }u = [R(fi), -^(^2)] and 

[[RM,RM],R(tp 3 )] = [Riip^R^RfalRfa)])] - [R(<p 2 ), R{[RM, R(f 3 



for any (pi,(p2,(p3 G sym^L- Clearly, the above equality holds by virtue of the Jacobi 
identity, since [im R, im R] C im R and R 2 = R. 

This example describes a rare situation when a Poisson structure is introduced 
on sym^L and not on H n (ir). The price we pay for this is that the underlying al- 
gebra a = R and an M-module structure of sym^L are very poor in comparison with 
J-(tt) for H n (n), respectively. 

Suppose a Frobenius operator is divisible by another Frobenius operator; then their 
Koszul brackets are correlated as follows. 

Proposition 17 (The chain rule, c.f. [SI]). Suppose that the Frobenius operator A is 
nondegenerate, see (1591) . and the image of A' = A o A is closed w.r.t. the commutation. 
Then the brackets {{ , }} A and {{ , }} A , are related by the formula 

A({{£i,£ 2 }} AoA ) = £A ( A(a))(A)(6) - £a(a (6 ))(A)(£i) + {{A(£i), Afa)}} A (78) 
for any sections £i,£2 £ f Q r(7r^(£)) of a horizontal fibre bundle in the domain of A. 

This assertion provides a considerable simplification of the search for new Frobe- 
nius operators using known ones and allows to reconstruct Sokolov's brackets in the 
inverse images of Frobenius recursions R G End]R x(7r), factored by Frobenius operators 
A: fi 1 (^ 7r ) — > x(7r). A weaker version of Proposition [T71 was formulated in [54| (where 
the operators were considered in a fixed system of local coordinates). This property is 
valid in the Hamiltonian case although is not well known. 

Proof. Suppose £1, £2 £ f and ipi = A(£j), (fi = A(ipi) for i — 1, 2. We have 

[^2] = (A o A)(£^(£ 2 ) - £ V2 (£x) + {{£i,£ 2 }Ua)- (79a) 
On the other hand, we recall that ipi — A(£,) and deduce 

[<p u <p 2 ] = A(E vl (i> 2 ) - £^1) + {tyiMU) (79b) 
= (AoA)(£^ 1 (£2)-£ V2 (£ 1 ))+A(£ Vi (A)(£ 2 )-£ ¥ , 2 (A)(£ 1 ) + {{^,^}}a). 
Now subtract (179aH from ( 179bl) and, using the nondegeneracy ( 1591 and hence omitting 



the operator A, we obtain the assertion. □ 

Example 13. Historically, the class of scalar operators, first regarded in view of §1 in 
Definition IU was studied in [60] using local coordinates. In [51] , Gelfand's symbolic 
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method [131 was applied in that setting, and it was argued that the class is infinite and 
contains the (presumably for A\ etc., well-defined Frobenius) operators 

A? = D X1 A? = D x o (D x + it) 5 4 3) = A? o (D x + it) ; 

Af = A? o (D x + u), Af = Af o (D x + 2it); 

here the superscripts denote the differential order and the subscripts enumerate opera- 
tors of equal order. Further, there are two operators 

A (n) = A (n-1) q + A (n) = q ^ + ^ _ ^ 

for any odd n > 5, and there are four operators 

4 n) = A f- 1] o p s + it) , A™ = AS n - 1} o (D x + 2u), 

AS n) = 4 n_1) °(^ + (n-3)«), 4 n) =4 M °(^ + («-2)?t) (80) 

for any even n > 6. These operators are homogeneous w.r.t. the weights \u\ = \D X \ = 1, 
hence the weights coincide with the differential orders. 

None of these operators is Hamiltonian. The brackets {{ , are reconstructed 

from the vanishing bracket for A\ by using the chain rule (see Proposition fTTt) ; the 
first four of them are 

{fe <?}} A ( 2 > = Pxq ~ pq x , 

{{Pi <?}} A (3) = 2u (PxQ - PQx) + PxxQ ~ PQxx, 

{{P> <?}} A (4) = 3(u 2 + u x )(p x q - pq x ) + 3u(p xx q - pq xx ) + p xxx q - pq xxx , 

q}} A W = 6 (« 2 + 2u x )(p x q - pq x ) + 6u(p xx q - pq xx ) + 2(p xxx q - pq xxx ) + p x q xx - p xx q x 



It remains unknown whether operators ( 1801) . which are factored to products of primi- 
tive first order operators with integer coefficients, exhaust all Frobenius operators with 
differential polynomial coefficients and homogeneous w.r.t. the weights \u\ = \D X \ = 1. 
The transformation rules for these operators are not known and the substitutions w for 
them, if they are not of the third (degenerate) kind, have not been constructed. 

Consider an operator R = A 2 o to factored by a Hamiltonian operator A 2 : k(ir) — > 
x(7r) and a pre-symplectic structure u: x{ji) — > x(7r) whose image is closed w.r.t. Dorf- 
man's bracket [, ]^. The recursion operator R is Frobenius, but the bracket {{ , }}r, 
which can be reconstructued using (1781) . is generally nonlocal. The exception is given 
by operators u of zero differential order if the determinants of their matrices do not 
vanish identically. 

Example 14. Again, let us consider Hamiltonian and Noether operators (1571) for 
the dispersionless 3-component Boussinesq-type system (|5*U|) . which is a true source 
of illustrative examples for the exposition! 15 ! The image of the Hamiltonian opera- 
tor Aq\ cosym£ — > sym£ is closed w.r.t. commutation ( 1291) . Note that the determinant 



15 The operator was labelled 'first' and A\ was 'second' with the lexicographical order w -< u -< v 
of their arguments ip = * [ip w , ip u , ip v ) & and images <p — * (ip w , tp u , ip v ) 6 k(tt) in [21] . 
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of the zero-order operator (I57cl) does not vanish identically, hence lo = A± : sym£ — > 
cosym£ is an isomorphism of modules for an open dense subset of solutions s G r(7r) 
of (I5TJ]) . Therefore A^ 1 is a Frobenius structure w.r.t. [, ]^ . By this argument, define 

R = A o A^ 1 : sym£ — > sym£. (81) 

This recursion operator for ([5"E]) is Frobenius. Indeed, its image is contained in the image 
of the Hamiltonian operator A , whose Dorfman's bracket [, ]^ is pushed forward by 
the zero-order operator Ax to [, ] Ro on sym£. 

The bracket {{ , can be calculated explicitly, although the resulting formulas are 
relatively long. For that, we use a Maple program for the Jets environment [17]; 
the listing is contained in appendix 03 Each component of Sokolov's bracket {{ , }} Ro 
for (1811) is a fraction of differential polynomials and contains about 15, 000 summands 
in the numerator. The denominators are the cubes of the determinant of the matrix A\. 
Perhaps, these fractions are reducible. 

In the same way, we construct the Frobenius recursion R 2 = A 2 o A^ 1 for floTJl) using 
its second Hamiltonian structure f!57dj) and the bracket^ {{ , }}s = {{ , }}a ± + {{ , }}a 2) 
see (11281) . The bracket {{ , }}r 2 is computed using a slight modification of the program 
that is contained in appendix [Dl The modification amounts to a substitution of A 2 and 
{{ > })i 2 for and {{ , }} if) , respectively. 

Finally, we specify a condition upon the brackets {{ , }} R e for sequences of Lie alge- 
bra structures [, ] R e induced on x(tt) by Frobenius iterations R £ of a recursion R G 
End]fjx(7r). The condition flH2"j) is a consequence of Proposition [T71 

Proposition 18. Let R: x(ir) — > k{ti) be a Frobenius operator and assume that 
its powers R 2 ,. . . ,R k are Frobenius for some k > 1. Then the brackets {{ , }} R e for 
each i G [1, . . . , k) satisfy the relation 

1-2 

+ {{R e -%,R e -%}} R , (82) 

where £i, £ 2 belong to the domain of R e . 

Under assumption that R, . . . ,R k are Frobenius for k > 1, is there any condition for 
R +1 , . . . , R k+e to be Frobenius as well? A well-known condition for Nijenhuis opera- 
tors [HI [35] originates from the standard cohomology theory for Lie algebras [HI [12], E] . 
At the same time, the nontrivial finite deformations [, ] R k of the standard Lie alge- 
bra structure [, ] on q(tt) are not the trivial infinitesimal deformations [, }^e, see ( Tl5|) . 
which are obtained using powers of Nijenhuis recursion operators N G End^ x(tt) with 
vanishing Nijenhuis torsion [N, N] fn = 0. 

4. Compatibility of Frobenius operators 

We introduce two types of compatibility for Frobenius operators Ax, . . ., A N : fi 1 (^ 7r ) — 
g(7r), where now fi 1 (^ 7r ) = f/ f)f =1 ker A{. The linear compatibility means that linear 

In the next section we show why this compatibility occurs. 



e R ^ l) (R)(m 2 )-8. Ri(! : 2) (R)(R^x) 
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combinations of operators remain Frobenius and hence the 'individual' Koszul brack- 
ets are correlated. The strong compatibility of N Frobenius operators means that the 
sum of their images is an involutive distribution in the Lie algebra q(tt) of evolutionary 
vector fields. We endow the spaces of both linear and strong compatible Frobenius 
operators with a bi-linear bracket that satisfies the Jacobi identity, and then we relate 
the Lie-type algebras of operators to an affine geometry with bi-differential Christoffel 
symbols. Let us consider these notions in more detail. 

4.1. The linear compatibility. First recall that any linear combination of two com- 
patible Hamiltonian operators is Hamiltonian by definition and hence Frobenius. Ex- 
ample [10] shows that a Hamiltonian operator can be decomposed as a sum of operators 
that satisfy (EE]), e.g., Af dY = -\ ■ D% + 2 ■ uD x + u x is a linear combination of ei = D x , 
e2 = uD x , and = u x . The decomposition may not be unique due to the existence of 
several linear dependent Frobenius operators that appear in the splitting; indeed, the 
operator A^ dY can be also obtained using v4 4 = uD x 

Definition 5. Frobenius operators A±, . . ., An are linear compatible if their arbitrary 
linear combinations A x = YliLi are Frobenius for any A £ M. N . The operators are 
linear compatible at a point X G M N if A^ o is Frobenius for a fixed linear combination. 

Example 15. There are two classes of pairwise linear compatible Frobenius scalar 
operators (j74p of first order. The first type of pairs is 72 (u) = const • ji(u) with any 
cki, ot2 G M and arbitrary functions Pi(u), ^{u). The second class is given by letting 
cti = a.2 G R, while the functions Pi, Pi, 71, and 72 remain arbitrary. 

Proposition 19. The Sokolov bracket induced on the domain of a linear combination 
A^ = J2iLi of linear compatible Frobenius operators is 

N 

= £v{{,}w 

Proof. This is readily seen by inspecting the coefficients of Af in the quadratic polyno- 
mials in Aj that appear in both sides of the equality [A^(p), A^(g)] = A x {\p, q]A x ) upon 
the Koszul bracket, here p, q G Q 1 ^). 

Consider the commutator XiA^p), ^ . XjAj(q)\ . On one hand, it is equal to 

= J2*MMpIMq)} +J2 X *M e >A i ( P )(q) - £A i(q) (p) + {{p,q}} At ). (83) 

i^j i 

On the other hand, the linear compatibility of A4 implies 

= (£a xW (?)) - A % (8. Ax{q) (p)) + A x ({{p, q }} Ax ) . 

The entire commutator is quadratic homogeneous in A, whence the bracket {{ , }}a x is 

linear in A. From (I83p we see that the individual brackets {{ , }} A . are contained in it. 
Therefore, 

{{p, q}}A X = £ V {{P, q}}A e + £ X * ' 7<(P> ?)> 
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where 7^: fi 1 (^ 7r ) x fi 1 (^ 7r ) — > fi 1 (^ 7r ). We claim that all summands •), which do not 
depend on A at all, vanish. Indeed, assume the converse. Let there be £ G [1, . . . , N] such 
that jeip, q) 7^ 0; without loss of generality, suppose t = 1. Then set A = (1, 0, . . . , 0), 
whence 

\Ai{p)^ X 3 A A<l) 

i j 

+ (A1A1) (£(AiAi)(p)(g) - ^(X 1 A 1 )(g)(p) + \i{{p,q}} Ai y 

Consequently, 7^0°, <?) G ker/bj for all p and q. Since each Ai is nondegenerate by 
assumption, we have that 7^ = for all £, which concludes the proof. □ 

Corollary 20 (Infinitesimal deformations of Frobenius operators) . Two Frobenius op- 
erators are linear compatible iff for any p, q G fi 1 (^ 7r ) one has 

[B(p),A{q)] + [A{q),B(p)] = A([p, q] B ) + B([p, q) A ) , 

which is equivalent to the relation 

E A(p) (B)(q) + 8, B(p) (A)(q) - E Aiq) (B)(p) ~ ZB {q) {A){p) = A({{p,q}} B ) + B ({{p, q}} A ) . 

Example 16. The three Noether operators Aq, Ai, and A2, see f )571) on p. [25], are linear 
compatible Frobenius structures for system ( 1561) . Any linear combination \ A + \iAi + 
X2A2 is Frobenius again, and therefore Sokolov's brackets for operators ()57cH57eT) are 
correlated by 

= {{,}}* + {{,}}*, 

which we claimed in Example [TH 

Remark 12. The operators remain Frobenius when multiplied by a constant, there- 
fore pass to the projective frame A G WP N of iV G N Frobenius operators. Then in 
CDiff (fi 1 (^ vr ), x(7r)) there is a basis of Frobenius operators which either are isolated 
points or which generate Frobenius cells with a nontrivial topology of attaching the 
simplexes together. 

An illustration is given by Example [TOl For e*i = D x , e 2 = uD x , and £3 = u x , the 
curve A\^ AV = (A : 2 : 1) is Hamiltonian and the ray A4 = f(u) -(0:1: —1) is Frobenius. 

4.2. The strong compatibility. We impose an additional specification on the struc- 
ture of the commutators of evolutionary vector fields whose generating sections belong 
to images of several Frobenius operators. 

Definition 6. Frobenius operators Ai, . . ., A^\ fi 1 (^ 7r ) — > x(7r) are strong compatible 
if the commutators of evolutionary fields in the images of any two of them belong to 
the sum of the images of all the N operators such that, for any p,q G fi 1 (^ 7r ) and 
l<i,j<N, 

N N 

[A t {p),A 3 {q)} = A 3 (Z Ai[p) {q)) - A(8, Aj{q) (p)) +Y,A k {r* J (p,q)) G ^imA, (84) 

fc=i i=\ 



(X 1 A 1 )(p),(X 1 A 1 )(q) = (X 1 A 1 )(X lll (p,q)) 
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Example 17. Consider the Liouville equation U xy = exp(2U). Let q x be the same 
projection onto the ^-characteristics as in (1771) . and similarly for q and the coordinate y. 
Frobenius operators □ = (U x + \DA o q x and □ = (U y + ^D y ) ° S x , c.f. ([HD, are strong 
compatible. The bi-differential coefficients are given in ( 124"1) on p. [121 

Definition (continued). The common domain f2 1 (^ 7r ) of the operators Aj is an F^)- 
submodule. Therefore, in view of the functional arbitrariness of sections p, q G il 1 (£„-), 
we say that the involutive distribution of evolutionary vector fields in the images of 
linear independent strong compatible operators Ai, . . ., An has reduced dimension (the 
rank) N. 

Example 18 (The Magri schemes). Completely integrable bi-Hamiltonian hierarchies 
determine commutative distributions of reduced dimension N = 2. All commutators 
vanish for the restrictions A\, A 2 of Hamiltonian operators onto the hierarchies; here 
l<i,j,k< 2. 

In other words, the strong compatibility of Hamiltonian operators A^ : Q 1 ^) — > x(-7r) 
is achieved on linear subspaces of Q 1 (7r). This is valid for the linear span of the gradients 
ipi = E(7ij) of the Hamiltonians Tii G H n (7r) which descend from the Casimirs Tio G 
Hq 1 ^) in the Magri scheme, see Theorem [TOl Indeed, one has im^ C imAi whenever 
both Hamiltonian operators are restricted onto the descendants of the Casimirs for A±, 
and hence the commutators (however, which vanisb[3 by the same theorem) belong to 
the image of A\. Thus the iteration (1491) of the Magri scheme corresponds to involutive 
distributions of reduced dimension two. 

Let us formulate the properties of the bi-differential symbols T\- G CDiff ((^) 2 fi 1 (£„.), 
^ (£*•))> which are determined by strong compatible Frobenius operators A±, . . ., An- 
Note that the symbols T^- depend on a point 9 G J°°(7r). 

Property 1. By construction, for any number iV of strong compatible Frobenius op- 
erators, we have 

r| = # • {{ , }} Al , for each i, 1 < i < N. (85) 
Hence a Frobenius operator yields the involutive distribution of reduced dimension one. 

Property 2. The symbols are not uniquely defined. Indeed, they are gauged by 
the conditions 

N 

E^^W^-^ll^ + rJfe?)) =°> P,qeQ\^); (86) 

k=l 

again, we assume the semi-simplicity of all the images in q(tt): [Ai(ip), q(tt)\ = implies 
ip G ker A%. 



For example, the first and second Hamiltonian structures for the KdV equation (TIT]) are not 
strong compatible unless restricted onto some subspaces of their arguments. On the linear subspace 
of descendants of the Casimir [w dx] , we have im A2 C im A\ and, since the image of the Hamiltonian 
operator A\ = D x is closed, we have [im A\, im A%] C imii. We emphasize that we do not exploit the 
commutativity of the flows. 
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Property 3. If, additionally, two strong compatible Frobenius operators Aj and A, are 
linear compatible, then their Sokolov's brackets are 

{{P, v}}Ai = Iy(p, q) + ^(p, q) and {{p, q}} Aj = T^p, q) + T^p, q) 

for any p, q G ft 1 ^)- 

Proof. For brevity, denote A = A iy B = Aj and consider the Frobenius linear combina- 
tion p,A + z/B. By Proposition [T9l we have 

(fiA + vB) ({{p, q}}^ A +vB) = 

= fJ. 2 A({{p, q}} A ) + p,v ■ A({{p, q}} B ) + p,v ■ B({{p, q}} A ) + u 2 B{{{ Pl q]} A ). 
On the other hand, 

(vA + uB)(p),( f iA + vB)(q) 

= ^ 2 [A(p),A(q)] +tiis[A(p),B(q)] - fiu[A(q), B(p)] + v 2 [B(p), B(q)]. 
Taking into account f[8~4T) and equating the coefficients of fiu, we obtain 
A({[p,q}}B)+B({{p,q}U)=A(ri B (p^ 

Using the obvious formulas T AB (q,p) = — T BA (p, q) and T AB (q,p) = — T BA (p, q), see ( 1SS1) 
below, we isolate the arguments of the operators and obtain the assertion. □ 

Definition 7. From now on, we consider Frobenius operators that are both linear 
and strong compatible. Such operators will be called totally compatible. By definition, 
totally compatible operators span linear spaces 

N 

A = @A k -R (87) 

k=l 

of Frobenius operators. 

Example 19. The operators □ and □, which were introduced in Example [I7J are 
totally compatible. This construction admits a straightforward generalization (11051) for 
other Euler-Lagrange Liouville-type systems, see Theorem [2HJ 

Remark 13. Within the Hamiltonian approach [22], it is very productive to think that 
the cosymmetries ip G x(7r) are odd. Indeed, in this particular situation the homomor- 
phisms ip G Hom^^) (x(ti), A n (7r)) are identified with Cartan 1-forms times the volume 
form dvol(M") for the base of the jet bundle. 

We preserve this understanding for domains fi 1 (^ 7r ) of Frobenius operators (1651) of 
second kind. The new Z-grading must be never mixe with any (e.g., Z 2 -) gradings 
of any variables on J°°(it) and with the arising gradings of the sections F G f) and 
ip G f) (here f = f)). Hence if tt is a super-bundle with Grassmann- valued sections, 
then Frobenius operators A are bi-graded [39]. Their proper Z-grading is \A\% = —1 
because the images in q(ti) have degree zero; the Z 2 -degree of the operators A can 
be arbitrary. In particular, Hamiltonian operators P produce Cartan's 0-forms by 
£ P( .): ^(vr) -> 0(tt). 



'We owe this remark to Yu. I. Manin (private communication). 
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In what follows, we assume for simplicity that all coordinates on J°°(7r) are per- 
mutable, whence fi 1 (^ 7r ) is even w.r.t. the Z 2 -grading. 

Property 4. For any i, j, k e [1, . . . , N] and for arguments p,q G ^ 1 (^7r) of Z-degree 1 
for strong compatible Frobenius operators of second kind, we have 

r* (p,g) = -r*(g,p) = (-l)W--lffl« . r* (g,p) (88) 

due to the skew-symmetry of the commutators ( |29i) . Hence the symbols T^- are sym- 
metric w.r.t. the Z-grading in the case (I65I) . 

Now consider the space of flat connections V Ak defined by 0691) for each Frobenius 
operator A^. Let us reveal the standard behaviour of the Christoffel symbols T^-. 



Property 5 (Transformations of Iy). Let w = w[w) be a nondegenerate change of 



fibre coordinates in the bundle £. Recall that the sections p, q 6 f are reparametrized 
by p i— > p = E(p) and g t— > g = 5(g), where H = £^ for Frobenius operators ( J63l) of first 
kind and H = 1 f° r ^ ne operators of second kind. Consequently, Frobenius 

operators A 1: . . ., A N : fi 1 (^ 7r ) — > q(tt) with a common domain fi 1 (^ 7r ) = f/f^ker-Aj 
are transformed by A i— »• As = A; o r .,. Then the bi-differential symbols 



r§ G CDiffp 1 ^) x fi 1 ^) - fi 1 ^)) obey 

r£(p,?) -r|( P ~g~) = (soi^Xh-^h- 1 ^^^ 

(89) 

This is a direct analogue of the standard rules r = 2 T 2 1 + dS ■ H 1 for the connection 
1-forms T under reparametrizations 5. Formula (18"9"1) represents this rule in the case of 
connections over the infinite jet bundles. 

Proof. Denote A = Ai and B = Aj; without loss of generality, assume i = 1 and j = 2. 
Let us calculate the commutators of vector fields in the images of A and B using two 
systems of coordinates in the domain. Then we equate the commutators straighfor- 
wardly, because the fibre coordinate in the images of the operators is not touched at 
all. So, we have, originally, 

TV 

[A(p),B(q)] = B(E A(p) (q)) - A(£ B{q) (p)) +A(ri B ( P , q )) + B(r%(p,q)) +J2M rk AB(p 

k=3 

On the other hand, we substitute p = E(p) and q = E(q) in B(q)] , whence, by 

the Leibnitz rule, we obtain 

[A(p),B(q)] = B(E, Am (E)(q)) + (BoE)(E Am (q))-A(8, m (E)(p)) + (ioH)(£ 5(a) (p)) 

JV 

+ {A o H- 1 ) (r^(Sp, Eq)) + (B o S" 1 ) (r|^(Hp, Eq)) ^) ^ 
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Therefore, 

rluM = o r^)(s P ,sg) - (s- 1 o S fl(g) (H))(p), 

rf B (p,g) = (S- 1 or^)(Sp ) Sg) + (E- l o8. A{p) (E))(q), 
T k AB (p, q) = (E- 1 o ry (Hp, Sg) for k > 3. 

Acting by S on these equalities and expressing p = H _1 p, g = we conclude the 

proof. □ 

Corollary 21. The bi-differential symbols 1^ constitute symmetric flat connections 

V(A) = J2k ^k^ Ak in the graded triples (Q 1 ^), q(it), A) determined by Frobenius 
operators of second kind. 

Remark 14. A straightforward calculation shows that two operators (|74|) are strong 
compatible if and only if they are proportional. Next, consider operators (|57p for the 
Boussinesq-type system (15 6p. If these linear compatible Frobenius operators are not 
restricted to a subspace of ^(S), and thus the arguments of Ak are generic, then the 
commutators [im Ai, im Aj] , < i < j < 2, are not decomposed using ( 184")) w.r.t. the 
images of any two operators Ai and Aj. We argue that these negative examples have a 
very deep motivation, which has been indicated in Example [T8j 

Indeed, by definition, the commutator f l84l) always contains the standard first two 
terms in the r.h.s. and takes values in the entire sum of images ^ fc imAfc for generic 
p, q G Q 1 (£„-) . Hence the commutation relations, which determine the Lie-algebraic type 
of the involutive distribution in the sum of images of Ai, . . ., A N , depend on a linear 
subspace S C fi 1 ^) that contains p and q. We may choose it ourselves in such a way 
that some of the operators become restricted on their kernels (up to the condition (186jl ). 

Therefore we expect that the strong compatibility and the decompositions w.r.t. the 
three linear independent operators (ISTI) are restored for restrictions of the Frobenius 
operators onto the spans of Hamiltonians for the dispersionless 3-component Boussinesq 
hierarchy. This will be discussed elsewhere. 

4.3. Algebras of Frobenius operators. Linear spaces of recursion operators R e 
EndiR x(7r) are equipped with the associative composition o, and therefore the recursions 
constitute monoids. Their unit is the identity mapping, and there appear relations 
between the operators or between their restrictions onto differential equations. For 
instance, the structural relations for for recursion operators of the Krichever-Novikov 
equations are described by hyperelliptic curves, see [7j. Taking compositions of two 
recursions Ri and Rj, we also obtain their formal commutators (I3"5"j) by setting 

[Ri, Rj] =Ri<3 Rj - Rj o R^ (90) 

Nontrivial examples of relations for the algebra structures (1901) are known, e.g., one 
has [i?i,i? 2 ] = R\ for the dispersionless 3-component Boussinesq system fl5*B1) . see [24J. 
In this way, the Richardson-Nijenhuis bracket (I3"5"j) endows the linear spaces of recur- 
sions with a graded Lie algebra structure (39J, HO]. If the Nijenhuis torsion [N, N] fn 
vanishes for an operator N: x{n) — > x(tt) and the Frolicher-Nijenhuis bracket [, ] fn , 
then it produces trivial infinitesimal deformations (fToT) of the standard bracket [ , ] on 
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another linear space, that is, on x(7r). By their turn, Frobenius recursion operators (1751) 
determine nontrivial finite deformations [, of the Lie algebra structure of q(tt). 

We propose a reverse scheme that starts with the standard structure (T29"]) in the 
Lie algebra g(7r) and then endows linear spaces of Frobenius operators with a Lie-type 
bracket. 

First, let Frobenius operators A\, . . . , An : f — > g(ir) be totally compatible such that 
each point of the linear space ( 1H71) is a Frobenius operator. Consider the commutation 
closure relations (184"1) . which are specified for vector fields with generating sections 
that belong to the images of A4. These relations express the decomposition of the 
commutators in the left-hand side w.r.t. the images again. Note that formula (1841) is 
linear w.r.t. each operator from A. 

We suggest to take the decomposition (1H41) and define the commutation rules 

[A l ,A J }(p,q)=[A l (p),A J (q)] (91) 

on A, here 1 < i, j < N and p,q G fi 1 (^ 7r ) = f/ f\kerAj. The structural constants 
of the algebra A are encoded by the bi-differential symbols r^-, whose properties were 
described in section 14.21 By construction, the Jacobi identity holds for the bracket (1911) 
which, under % «-> j, p «-> q, is skew-symmetric for Frobenius operators (lolll) of first 
kind and which is symmetric (in the graded sense \p\% = \q\% = 1, see Remark [T3l) for 
Frobenius operators ( 1651) of second kind. 

Consequently, each Frobenius operator Ai G A spans the one-dimensional algebra, 
and the Koszul bracket [ , ]a, measures its noncommutativity. The assumption that the 

operators are linear compatible implies that each line (A), where (A) = f IR'X^i C 
A, is a one- dimensional subalgebra of A. 

By this argument, we obtain the third operation on the space ( 1H7|) of totally compat- 
ible Frobenius recursion operators x{ir) — > x(7r), in addition to the composition o and 
the commutation (190]) . At the same time, if the Frobenius operators are not recursions, 
then the compositions of the operators Ai : f — > k{tt) are not defined, and the algebra 
structure (19 ip on the linear space A is the only one remaining. 

Example 20. Proposition [8] proves the existence of one-dimensional algebras of Hamil- 
tonian Frobenius operators of second kind. Theorem [10] yields the two-dimensional 
algebras A, see Example fT8l Restricting the Poisson pencils (Ai,^) onto the hierar- 
chies, we obtain analogues of the solvable two-dimensional Lie algebra with a relation 
[ai, a 2 ] = a\. Moreover, the Magri scheme shows that these algebras A are commutative. 

Second, let us consider a wider class of operators A»: f — ► k(tt) that may not be 
Frobenius. Namely, let all the operators be transformed by either (IUH|) or (|65p . which 
makes them well defined, and suppose that the commutation closure 

im Ai ' J2 im M - J2 k im Ak ( 92 ) 

is valid for the images of the whole iV-tuple (Ai, . . . , An). This incorporates the previous 
case of totally compatible Frobenius operators that satisfy (I6TI) . Now we introduce a 
Lie-type structure on the linear subspace A C CDiff (f — > *c(ir)). 
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Let [Ai, Aj] be the same bracket 
other notation, by ( I92l we assume thai 



[Ai, Aj] 



of the operators A\ 



N 



A k o c 



k 

iji 



Am G A at hand. In 



(93) 



k=l 



where c^- : f x f — > f are the bi-differential structural constants of the new algebra A. 
Obviously, the Christoffel symbols r*-, which were introduced in (JM1) . are encoded by 
cfj up to the gauge (1861) . The structural constants are symmetric or skew-symmetric 
w.r.t. the lower indexes simultaneouly with T^. Admitting a slight abuse of language, we 
continue calling these algebras A with relations f[9"3l by algebras of Frobenius operators, 
although the image of an operator Ai may be not closed under the commutation for 
some £ e [1,...,N]. 

For instance, we have c^. = mod ker A& for the algebras A of rank two that 
are generated by restrictions of Poisson pencils (A\,A2) onto the Magri schemes, see 
Example [181 Other examples of such algebras A will be given in the next section, where 
we assign Frobenius operators (I105P of second kind to Liouville-type integrable systems 
(in particular, to the 2D Toda chains associated with complex semi-simple Lie algebras, 
see ([U). 

Thus we arrive at the structural theory problem for the operator algebras A with 
generators Ai : f — > and relations (l93l . 

Conjecture 22. There are horizontal jF(^)-modules f c — > r(7r^(£)) and linear differ- 
ential operators A\, . . ., A^: f — > 53(71") such that, by a choice of appropriate subspaces 
S C f in their domain, one recovers a bi-differential extension of the classical struc- 
tural theory of (e.g., nilpotent, semi-simple, or solvable) Lie algebras of reduced dimen- 
sion N for the linear spaces A, which are generated by the restrictions of the operators 
A 1: . . . , A N onto S, and for the algebra structures (I9"B1 on them. 

Remark 15. Recursion operators for differential equations S = {F = 0} can be under- 
stood as Backlund autotransformations between symmetries, which are solutions of the 
linearized systems £f{<p) = 0, see [5], [39]. Hence the structures ( 1931) of the algebras A 
of Frobenius recursion operators for 8 are inherited by these classes of Backlund auto- 
transformations. This determines the brackets between non-abelian coverings over the 
equations £. 



5. Factorizations of symmetries of Liouville-type systems 

In this section we describe an infinite class of Frobenius operators □ and calculate the 
brackets {{ , }}□ induced by them. These operators appear in the description of symmet- 
ries of the hyperbolic Liouville-type Euler-Lagrange nonlinear systems [HI EU EHl ED] . 

To start with, we extend the collection of known Frobenius operators with the one 
that factors point symmetries of the non-evolutionary (2 + 1) -dimensional A^-Toda 
equation. 

19 Thus we postulate that the bi-differential dependence on the arguments in is entirely absorbed 
by the structural constants cf . 
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Example 21 Q28J). Generators of the point symmetry algebra for the 'heavenly' Toda 
equation u xy = exp(—u zz ) have the form cp x = CP(0(:e)) or ip y = where 
0,0 G C°°(M) and 

LT = u x + \z 2 D x , O y = u y + \z 2 D y = (x^ y) (U x ) . (94) 

Clearly, the commutator of any two point symmetries of x- or y-type is a point symmetry 
again such that the action of the operators □ on the spaces of the free functional 
parameters 0, is given by the Wronskian, {{0i, 02 j}^ = 0i • (02) x - (0i) x • 02- 

Operators (19~4"1) factor the right-hand sides of the symmetry flows u t = 0(0) on the 
heavenly equation. This means that a bigger differential equation is the representing 
object for a transformation group generated by the flows. This approach to constructing 
Frobenius operators is very productive. First, let us recall a procedure that assigns 
hyperbolic Euler-Lagrange systems to hierarchies which are Hamiltonian w.r.t. A\ = 
const --Dx, see [25l [28] and also Remark[9]on p.[32j The method is based on the canonical 
coordinate-momenta formalism [31] for PDE and on a representation of the KdV-type 
hierarchies as commutative Lie subalgebras of Noether symmetries for the hyperbolic 
systems. 

Example 22. The dispersionless 3-component Boussinesq-type system, see p. [2H 

u t = ww x + v x , v t = -uw x - 3u x w, w t = u x , (ESD 

is not written here in the form of a conserved current only due to an unfortunate choice 
of local coordinates for it and for structures (1571) in [2^. Indeed, let us return to the 
notation of jH] and let the new dependent variables m = (u, v, w) be 

u = v + w 2 , v = w, w = u, (95) 

which are densities of the Casimirs w.r.t. the first Hamiltonian structure (157aj) for (I56p . 
Thence we obtain the system 

1 0\ (8/8u\ 

1 O A, S/6v (Hi), (96) 
1/ \5/8wJ 

where TCi = [( uw — ^v 2 w) dx] . Likewise, the dispersionless Benney equation [HI 
Eq. (19)] acquires the same form ( 1961) in these coordinates. 

Now we introduce the conjugate variables U = (U, V, W) such that 

(«")'=-a80^ 

in other words, the adjoint linearization (l^Y of the canonical momenta m w.r.t. 
the canonical coordinates U is proportional to the first Hamiltonian structure for ( 1961) . 
Thus we set 

u = V x , v = U x , w = W x . (97) 
The potential variables satisfy the system 

u t = w x , v t = -u x w x , w t = -\u 2 x + v x . pi) 
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Using the notation (j9~5|) and (j9~7|) together, we cast f[9"BT) and (1961 j) to the canonical form 
TT SHi TA 6H1 TJT SHx SHi 6H1 8H X 

Ut = -fo> Vt = -fo> Wt = l^ ] ^ = -lcT' ^"W 

(98) 

where H x = [(V X W X - \V 2 X W X ) dx] . 

Next, let us find the hyperbolic Euler-Lagrange equation 

£el = {Eu(£) = | C= [Ldxdy],L = -i(m,U 2/ )-H L (U)} 

such that the bi-Hamiltonian hierarchy of commuting flows for ( |98i) is composed by 
symmetries of £el- This is done straightforwardly. We note that system (1981) is scaling- 
invariant w.r.t. the homogeneity weights \u\ = 2, \v\ = 1, |w| = 3/2; |Z7| = 0, |V| = 1, 
|W| = 1/2; and \D X \ = 1, |D t | = 3/2. The symmetries^ of (EHD with time weights 
3/2 = \D t \, 2, and 3 fix the Hamiltonian Hl(U) for S^l uniquely. The only system that 
is ambient w.r.t. the whole hierarchy of (I98p is the wave equation 

£ ={U xy = O, V xy = 0, W xy = 0}. 

The wave equation £ is exactly solvable, because the conditions 

U x , V x , W x G kerD y \ £z , U y , V y , W y € ker^l^ (99) 

can be integrated immediately, and we obtain XJ = f(x) + g(y). At the same time, by 
the argument illustrated in this example, we have naturally arrived to the definition of 
a class of nonlinear hyperbolic systems such that conditions ( 1991 become nontrivial. 

Definition 8 Q60J). A Liouville-type system £ L is a system {u xy = F(u,u x ,u y ; x,y)} 
of hyperbolic equations which possesses the integrals w±, . . ., w r ; Wi, . . ., w f G C°°(£l) 
such that the relations D y \ £ Jwi) = and D x \ £ (iBj) = hold by virtue (=) of and 

such that all conservation laws for £l are of the form [f(x, [w]) dx\ © [g(y, [w]) dy\. 

Example 23. The m-component 2D Toda chains (fP2"l) associated with semi-simple 
complex Lie algebras [12] constitute an important class of Liouville-type systems, here 
u = (u l , . . . , u m ). By [56J, they possess the complete sets of 2m integrals W\, . . . , w m \ 
Wx, . . . , w m and are integrable iff K is the Cartan matrix. This class is covered by the 
ansatz in Proposition [23 below. 

Remark 16. Consider Liouville-type systems that possess complete sets of 2m integrals: 
r = m and f = m. By slightly narrowing the class of such equations, let us consider 
the systems £l whose general solutions are parameterized by arbitrary functions / 1 (x), 
. . ., f m (x) and 5 fl (?/), . . . , g m (y)- These equations are represented as the diagrams [26J 

(f) m J°°(ir x ) © J°°(<k v ) ^S L ^ fff J°°(ir x ) © J°°(^), 

where tt x and n y are the trivial fibre bundles R x 1 R such that /, w and g, w 
determine their sections, respectively; the first arrow is given by the formulas for exact 
solutions and the second arrow is determined by the integrals. 

20 For \s\ = 2, the symmetry of jM) is U s = V x , V s = |i7| - \Wl, W s = -U X W X ; the flow with 
|t| = 3 is given by U T = 2U X V X + W%, V r = \U* + - 2U X W%, W T = 2V X W X - 2U^W X . There is no 
symmetry of (|96l|) with time weight 5/2. 
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Remark 17. Consider now a class of hyperbolic quasilinear systems 8^ whose sym- 
metries if = □(</>) are determined by Frobenius linear differential operators □, c.f. 
Remark [3] on p. [2H1 Definition [H] of the Liouville-type systems, which postulates the 
existence of the integrals, is not equivalent to the definition of Soo with the functional 
freedom in symmetries. However, the two notions are very close. In Corollary EU see 
below, we prove that symmetries tp of Euler-Lagrange Liouville-type systems £l are 
factored by differential operators □, □ that originate from the integrals w, w for 8^: 
we hav^j] 

ip = D(0(x, M)), Cp = Q(0(y, [«>])), (100) 

where components of the sections <j) and <j) are arbitrary smooth functions. We show that 
the operators □, □ are Frobenius, which yields decomposable operator algebras A = 
□ • R © □ • R of rank two. 

Let us remark on the history of the problem. Scalar Liouville-type equations were 
studied in and their symmetries have been further analyzed in [8J. Operators □ that 
factor Noether symmetries of Euler-Lagrange Liouville-type systems u xy = F(u;x,y) 
were constructed in [25]. The general y) of the Euler-Lag- 

range Liouville-type systems (see (11041) ) was considered in [58]; however, no method for 
reconstructing the brackets {{ , }}□ is described there. 

The problem of construction of operators □ that assign (possibly, not all) symme- 
tries of non-Euler-Lagrange Liouville-type systems to their integrals w, w is much less 
transparent. A considerable progress has been achieved here in a recent paper [57], see 
section 15. 2\ where it is shown that the existence of differential operators □ is based on 
the existence of w, and respectively for □ and w. We analyse mainly the former case 
of Euler-Lagrange systems. 

Example 24. Consider the parametric extension of the scalar Liouville equation ([6]), 

u xy = exp(2w) • y/l + 4e 2 u 2 x . (101) 

This equation is ambient w.r.t. the hierarchy of Gardner's deformation of the potential 
modified KdV equation, see [26]. The contraction U — U(e, [u(e)]) from (11 Oil) to the 
non-extended equation U xy = exp(2U) is U = u + \ arcsinh(2£:u x ); it determines the 
third order integral for (1 10 II) using the integral ([7J) at e — 0, see Example [TJ However, 
the regularized minimal integral of second order for (1 10 II) is 

w = (1 - Ae 2 u 2 x ) /2e 2 + u xx /^/l + Ae 2 u 2 x ; (102) 

such that all other x- integrals for (II 011) are differential functions of (11021) . The second 
integral for (1 1011) is w = u yy — u 2 — e 2 - exp(4w) 6 ker D x \ g . The operators □ = u y + \D y 
and 

12e A u u 2 

□ = |(1 + Ae 2 u 2 x - 2e 2 u xx ) ■ D x + u x + Ae 2 ul - 2e\ xxx + - - J 2 J (103) 

assign symmetries (HOOp of (110 II) to its integrals. We emphasize that operators in 
the family f 1 1 3 j) assign higher symmetries ip = D(0(x)) of (1101)) to functions on the 

21 It is believed that Liouville-type systems do not possess any symmetries other than (|100[) . unlike, 
e.g., the KdV equation (fTTjl that admits the scaling and the Galilean symmetries in addition to its 
commutative hierarchy 21. 
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base of the jet bundle whenever e 7^ 0, while the operator □ always determines point 
symmetries (p = Q(0(?/)). 

Both operators □ and □ satisfy flBTl) . The bracket {{£>,<?}}□ = p y q — pq y for □ is 
familiar; the bracket induced in the inverse image of □ is calculated in appendix [Bl 
The surprisingly high differential orders of {{ , }}□ with respect to its arguments and 
coefficients is motivated by the presence of higher order derivatives of u in (11031) . 

Remark 18. Deformations of the algebras of Frobenius operators appear by virtue of 
the Gardner deformations for differential equations, see Example EH and [26]. Let S(fi) 
and £{y) be the extensions of an equation £(0) and let m £ : £{e) — > £(0) be the Miura 
contraction at e G I C R. Symmetries (p^ G sym£(yu) and ip v G sym£(z/) induce the 
symmetries of £(0) through the contraction, although the induced flows can be formal 
sums of infinitely many generators of higher orders. However, consider the commutator 
of the (formal) symmetries. For all e(/j,, v) such that the commutator is lifted to a 
true symmetry of £(e), we define the product o: (/^, v) i— > \i o v — e{\i,v). Thus we 
obtain the multiplication M o N = {e} C 1 on the sets M = {//}, N = {v} C 
IR. For the Liouville-type Gardner's extension (11011) . we further obtain the product 
o: (D(yu), D(z/)) t— > n(fioh') of the Frobenius operators (1 1 03H . 

5.1. Frobenius operators for Lagrangian systems. For Euler-Lagrange Liouville- 
type systems £el = {F = E(£) = 0}, the existence of factorizations (llOOp for at least a 
part of symmetries is rigorous and can be readily seen as follows. For integrals w such 

that D y (w) = V(F) and for any I(x, [w]), the generating section ipj — V* o [t^ )* ° 

(i^Y (1) for a conservation law [I dx] solves the equations = on £el- The 

Helmholtz condition ( 1331) for £e(c) implies that the vector 

¥#] = [V* o (4 U) )*] N)) G ker£ E(£) |^ L (104) 

is a symmetry of £el for any = (£y' )*(1) = E w (/dx). A standard reasoning (see 
Lemma [7] or Corollary [2H below) shows that formula (11041) yields symmetries of the sys- 
tem £ EL even if sections do not belong to the image of the variational derivative E w . 

In this section, we construct a class of Frobenius operators □ associated with sym- 
metries of Euler-Lagrange Liouville-type systems. We express the operators □ and 
the brackets {{ , }}□ in terms of minimal integrals w. This is done by the following 
argument. 

Proposition 23 Q25J). Let k be a nondegenerate symmetric constant real (m x m)- 
matrix. Suppose that C = [L dxdy] with the density L = — | J2i j K ij u% x u y ~~ Hl(u; x, y) 
is the Lagrangian of a Liouville-type equation £ L = {E(£) = 0}. Let m = dL/du y be 
Dirac's momenta [31] and w(m) = (w 1 , . . . ,w r ) be the minimal set of integrals for £^ 
that belong to the kernel of D y \ £ ■ Then the adjoint linearization 

□=(4 m) r (los) 

of the integrals w.r.t. the momenta factors all Noether symmetries <pc of £^, which are 
given by 

<p c = n(5H/8w) (106) 
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for any 7i = [H(x, [w]) dx] . 

Outline of the proof. The assertion follows from 

• the structure ip = — D~ x (E u (7i)) of the generating sections of conservation 
laws for hyperbolic systems u xy = k." 1 (SH^/Su) resolved w.r.t. the second-order 
derivatives, 

• the correlation 

between the variational derivatives w.r.t. u and w, and 

• the correlation if) = mpc between generating sections of conservation laws and 
Noether symmetries cp£ of systems u xy = k' 1 (5 Hi,/ '5u) , see Theorem El 

□ 

Corollary 24. Under the assumptions and notation of Proposition [23j the section 

ip = n((f)(x, [w])) e x(7r) 

is a symmetry of the Liouville-type equation £l for any = t (<f> 1 , . . . , r ) G f. 

Proof. The proof is standard and analogous to the one for Lemma [7J with the only 
alteration in the jet space at hand. Consider the jet bundle J°°(£) over the fibre bundle 

(:K r xK^K (107) 

with the base l9i and the fibres W with coordinates w 1 , . . . , w r . By Proposition 1231 
the statement of the theorem is valid for any <f> G f = in the image of the variational 
derivative E^, where w = w[m\. Obviously, the image contains all r-tuples whose 
components <p l (x) G C°°(R) are functions on the base of the new jet bundle. The 
substitution w = w[m[u]\ : f — > converts sections G f to elements of r(7r^ (^)). 
Now recall that □ is an operator in total derivatives ( |27|) . whose action on r-tuples of 
functions on the jet space J°°(7r) is defined by ( 1261) through their restriction onto the 
jets of sections w = <p(x), whence the assertion follows. □ 

Theorem 25. The operators (11051) are Frobenius of second kind if the integrals w are 
minimal. 

This assertion is continued in Proposition 1271 on p. [5*31 

Proof of Theorem [251 The commutator of two Noether symmetries (p' c , (p£ is a Noether 
symmetry </?£, and hence the conservation law corresponds to it. The geometry of 
the Euler-Lagrange Liouville-type equations £l — {^ _1 E U (£) = 0} is such that the 
conservation law is represented by an integral, D y (H) = on £ L . By assumption, the 

integrals w that specify the symmetries ip' c = □(0'[u>]) and ip" c = □(</>"[w]) are minimal, 
meaning that any integral is a differential function of them, hence H = H(x, [w]). 
Let the gauge of the minimal integrals be fixed. Then the factorization (I106P for the 
new symmetry ipc follows from Proposition [23 Under differential reparametrizations 
w = w[w] of the integrals, the sections = 5H/5w are transformed by = (^' > )*(0), 
thence □ is a well defined Frobenius operatos of second kind (|65p . □ 
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Remark 19. Frobenius operators □ generate symmetries of Euler-Lagrange Liouville- 
type systems using arbitrary r-tuples of integrals. Therefore let us consider columns 
of these (m x r)-matrix operators separately: 

□ = ((□!),■■■,(□,)). (108) 

Generally, the image of a A;-th column is not closed under the commutation. This is 
the case when the k-th components <p k of sections G f are coupled not only in the k-th. 
component of the bilinear bracket {{ , }}□. Likewise, the commutators of symmetries 
in the images of any two columns Dj are decomposed with respect to the images of 
other operators Di, . . ., D r as well. 

However, the decomposition (11081) does not produce operator algebras A of rank r. 
Indeed, under a reparametrization w = w[w] in the domain of a Frobenius operator □ of 

second kind, which yields [Do )*] (<j>) = Q(0), the columns LJ fc are not transformed 
individually. Therefore they are not an r-tuple of objects, but they do constitute a single 
well-defined operator (11081) . By construction, it generates the algebra A of rank 1. 

Example 25. Consider the Euler-Lagrange 2D Toda system associated with the simple 
Lie algebra sl 3 (C), see H21ES], 

£ T oda = {^ y = exp(2[/-V), V xy = exp(-U + 2V), #=(_?"!)}. (109) 

The minimal integrals for (11091) are [^31 156] 

w 1 = U xx + V xx - U 2 X + U X V X - V 2 , (110a) 
w 2 = U xxx - 2U X U XX + U X V XX + U 2 X V X - U X V 2 . (1 10b) 

All symmetries (up to x <-+ y) of (I109p are of the form = \3((j)(x, [w 1 ], [w 2 ])), where 
<f = t (4 >1 , <p 2 ) G and the (2 x 2)-matrix operator in total derivatives i@ 

n= fU x + D x 2D\ + 3U X D X + U 2 + 2U X V X - 2V 2 - U xx + 2V XX \ , . 

\V X + D X D l-2U XX + V xx + 2U 2 X -2U X V X -V 2 )' {LLL) 

By Theorem [23| the image of operator (II lip is closed w.r.t. the commutation. The image 
of the first column, of first order, is itself closed under commutation (see, e.g., [25]), but 
the image of the second column of □ = (di, D 2 ) is not. Indeed, the bracket {{ , }}□ in 
the domain of the entire operator □ is equal to 



= PW ~ P^l + 3p 2 ql x - 3p l xx q 



,2 



+ 6p 2 qLx - q pIxxQ 2 + Qujl ■ (phi - pW) , ( 112 ) 
{{p, q}}h = pW - p l( il + 2 pW - 2 p 2( iI + 3 pIx<i 2 - 3 P 2 <i 2 xx- 



Here we box the individual bracket {{ , for the (2 x l)-matrix operator Dl, and 
we underline the couplings of components in the domain of O^; under commutation, 



22 For convenience, we have multiplied the second column of operator (jllip by —3, which of course 
does not preclude it from assigning symmetries of p09[) to arbitrary arguments 4> 2 . However, this 
correction, which is introduced by hand, is the only reason why the operator (Ay) in Example 1281 on 
p.EUis not Hamiltonian; see [HI EI] for the bi-Hamiltonian pair of the Boussinesq hierarchy associated 
with the semi-simple Lie algebra 513(C). 
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they hit both images of Di and \Z\ 2 . The remaining terms specify the bi-differential 
structural constants, which we denote by T^- with an obvious meaning. We have 



r{ 2 = -3Dl 



i, = 3-ii 



-1®D X + 2D X 



r 2 

1 21 



Dt 01 — 2 • 1 £x) D q 



In Proposition [23 we explain why the coefficient &w l in T\ 2 is expressed by a function 
of integrals for the Liouville-type system 



Examples of Frobenius operators associated with the Dynkin diagrams via the 2D Toda 
chains (fT2"j) are discussed in [51] . 

The integrals w[m] of Euler-Lagrange Liouville-type systems Si determine the Miura 
substitutions from commutative mKdV-type Hamiltonian hierarchies 03 of Noether 
symmetries on Si to completely integrable KdV-type hierarchies 21 of higher symme- 
tries of the multi-component wave equations S = {s xy = 0}, see below. The potential 
modified KdV equation ([Toll , which is transformed to the KdV equation ( TTTT) by 0, 
gives a natural example. This is discussed in detail in [25] . 

The hierarchies share the Hamiltonians 7ii[m] = W[w[tti]l through the Miura substi- 
tution. The Hamiltonian structures for the Magri schemes of 21 and 03 are correlated 
by the diagram 

Ho Hi H 2 



> $2 ( > 02 (hierarchy 21) 

□ (113) 
(hierarchy 23). 



H Hi 

Frobenius operators □ map cosymmetries fa for the hierarchy 21 to symmetries (fi+i 
that belong to the modified hierarchy 23. The first Hamiltonian structure B\ = (^m^)* 
for 03 originates from the differential constraint m = dL/du y upon the coordinates u 
and the momenta m for £ L . 




Lemma 26. Introduce the linear differential operator 

A k = □* o B x o □, 



(114) 



which is completely determined by the Euler-Lagrange Liouville-type system Si,. The 
operator ( 1114ft is Hamiltonian and determines a higher (that is, k = k(\3,m) > 2) 
Poisson structure for the KdV-type hierarchy 21. 

Proof. By construction, the bracket (|3T|) given by satisfies the equality 

{WiH,W 2 H} ifc = {HAw[m]],H 2 [w[m}}} ^. 

Therefore the left-hand side of the above equality is bi-linear, skew-symmetric, and 
satisfies the Jacobi identity. Fourth, it measures the velocity of the integrals w along 
a Noether symmetry of Si. Since evolutionary derivations are permutable with the 
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total derivative D y , the velocity {^1,^2} A k ^ es m ker | £ ^ and hence its density is a 
differential function of the minimal integrals w. □ 

Proposition 27. The bracket {{ , }}□ in the inverse image of operator (11051) is equal 
to the bracket {{ , }}^ k , which is induced by the Hamiltonian operator Ak and which is 
calculated by formula f |46l) . 



{{0', 0"}}n = {{0', 0"}} A k = V), 0', 0" e cosymSl C sym£ . (115) 

The coefficients of the Hamiltonian operator Ak and of the bilinear terms in the bracket 
{{ ? }}□ ar e differential functions of the integrals w. 

The multi-component wave equation £ = {s xy = 0}, whose symmetries contain the 
hierarchy 21 and such that A\ = potentiates the image of the Miura substitu- 

tion, is not a priori unique. Again, the constraint between the coordinates s and the 
momenta w for 21 determines the first Hamiltonian operator A\ for 21, but the con- 
straint appears apparently from nowherd^l; the shift of the field or the frozen point 
argument [2] are customary procedures here. Our paradoxal conclusion is that the 
first structure Ai = A^ 1 for 21 is chosen such that Ai factors the higher Hamiltonian 
structure for 23. We thus have 

B k , = DoA 1 o □*, k' = k'(D, (4*))*) > 2, 

which specifies the required nonlocalities. 

This means that Frobenius operators and the factorizations they provide are helpful 
in the bi-Hamiltonianity tests for integrable systems. It is likely that one can reveal 
a similar origin of the nonlocal first Hamiltonian structures for the Drinfeld-Sokolov 
hierarchies [10J associated with the Kac-Moody algebras, whose Cartan matrices are 
degenerate. 

Reciprocally, formula (I115p calculates the bi-differential brackets {{ , }}□ induced by 
Frobenius operators fll05[) . Setting to zero all but one components of sections 0'- lj, <$>" ■ 
lj G f = x(£) in domains of operators □ (respectively, except z-th and j-th components, 
1 <hj < r), we obtain the coefficients in the bracket {{ , }}□. Proposition [2T1 shows 
that the functional class for these coefficients is very narrow. 

Example 26 (The modified Kaup-Boussinesq equation) . Consider an Euler-Lagrange 
extension of the scalar Liouville equation [2S], 

A xy = -\Aexv(-\B), Bxy = \ex V (-\B). (116) 

Denote the momenta by 

a = \B X and b = \A X . (117) 
The minimal integrals of system (11161) are 

Wi = —\a 2 — a x , w 2 = ab + 2b x , 

such that D y {wi) = on (I116p . i = 1,2. Hence the operator 

□ = ( et b LY = i ~^ B * + Dx 1 „ ^\ „ \ (118) 



' m ' W2 J \ ±B X -2D X 



2 



2:! 



This difficulty of the theory was pointed out to us by B. A. Dubrovin (private communication). 
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factors (Noether) symmetries of (11161) . The bracket {{ , }}□ induced in the inverse image 
of the Frobenius operator □ is 

{{^x}}n-5-^ r x X 2 -^ 2 xl )' ( 9) 

where ip = t (ip 1 , ip 2 ) and % — *(x\ X 2 )- Split the operator □ = (di, D 2 ) to two columns. 
The bracket (11191) determines the structural constants T^- for the commutation of images 
of D 1 , D 2 . For example, setting the components ip 2 -,x 2 an d i n (HI 91) to zero, 

respectively, we isolate the brackets 

{{ , }}□! =0, {{ , }}n 2 = \(D X ® 1 - 1 ® D x ) 

for Di and D 2 . Thus we have obtained an extension of the Wronskian-based bracket 
for the second Hamiltonian structure of KdV, see Remark [TUJ At the same time, 
bracket (11191) contains all other structural constants for the rank one algebra generated 
by □ = (Dx, D 2 ) , e.g., T\ 2 = \ (D x <g> 1 - 1 <g> D x ) , and T{ 2 = 0, whence Y\ x and !%_ are 
deduced by the skew-symmetry. 
Consider a symmetry of ( 11 161) . 

At = \A X A XX + | (jA 2 ~ l) B x , B t = —2A XXX + ^A X B 2 — ^A X B XX . (120) 

Using the constraint ( 11 17[) between the coordinates A, B and the momenta a, b, we 
cast (I120p to the canonical form 

sn sn n m . m /101 , 

A t = — , a* = ; B t = — , b t = , (121) 

* 5a ' * <L4 5b' SB' K J 

where W = [(±A*B% + \A X A XX B X + \A 2 XX - \B 2 X ) dx] . 

Let us choose an equivalent pair of integrals u = w 2 , v = W\ + \w\. It is remarkable 
that the evolution of u and v along ( 11201) . which equals (see p. [26]) 



u t = uu x + v x , v t = (uv) x + u xxx , 

is the Kaup-Boussinesq system, and (11201) is actually the potential twice-modified 
Kaup-Boussinesq equation. The right hand side of the integrable system (11201) be- 
longs to the image of the adjoint linearization □ = (^["'^l)*. The Frobenius operator 

□ factors the third Hamiltonian structure Af B — □* o (jfy^ff)* ° Q for (1581) ; we have 
k = 3 and 

aKB = f uD x + \u x Dl + (W + v)D x + \{u 2 + 2v) x 

3 ' \D't + {\u 2 + v) D x + \v x \(2u D\ + 3u x D 2 X + (3u xx + 2uv)D x + u xxx + (uv) x ) 

By Proposition [27J the bracket {{ , }}q is equal to {{ , }}^kb, which is given by for- 
mula (HoT) . We obtain 

= {{lx}}^ = (f ^<!> - ' !) , (122) 

where Vi = -§ d%+vd x ) and V 2 = -\{d x uD x )- Here we use an alternative 
notation (11221) for the components of the bracket {{ , }}q that acts by the differential 
operators Vi and V 2 on one factor in each coupling. 



INVOLUTIVE DISTRIBUTIONS OF OPERATOR- VALUED EVOLUTIONARY DERIVATIONS 55 



Finally, the operator A\ = ( £ D X ) is the first Hamiltonian structure for ( 1581) ; its 
inverse A\ = A^ 1 factors the second Hamiltonian structure B 2 = □ o A\ o □ for (11201) . 

We recall that the operators ( 1105!) determine recursions for the symmetry alge- 
bras sym £el of Liouville-type systems. 

Proposition 28 (see [28, 29J). Let £el be an Euler-Lagrange Liouville-type system 
that meets the assumptions of Proposition [23j w[m] be the integrals, and □ be Frobenius 
operators (11051) . Let the r-tuples <fi = t ((j> 1 , ■ ■ ■ , 4> r ) G f = M£) ^ e the same as in 
Corollary [2H see ( 1107P on p. [501 Then the m-tuples 

u $ = □(dc?(ac, Htn]])J 

of Cartan's 1-forms are generating sections of local recursion operators R?: ip i— > £^ i c^^- 

for <p G sym£ EL . 

Conjecture 29. Under the above assumptions and notation, the recursion operators 
constitute infinite-dimensional Frobenius-Lie algebras by 

<?ef 

— * — * I 

where 0', 0" G f = and p',p" G sym^L are any symmetries of the Euler- 

Lagrange Liouville-type systems Sel- 

5.2. Non-Lagrangian Liouville-type systems. Let = {F = 0} be a Liouville- 
type system; now it may not be Euler-Lagrange. Let w G kerD^J^. be a section 
composed by minimal integrals for £l, thence D y (w) = V(F), see ( 1341) on p. [151 By 
construction of Liouville-type systems, hyperbolic equations in £l are independent from 
each other. Therefore £l is both normal and ^-normal p3, HQ], meaning that A(F) = 
or Ao£ F = on 4 implies A = 0, respectively. By this argument, a section ip G x(ir) is 
a symmetry of a Liouville-type system £l if and only if the evolutionary vector field 6^, 
preserves the integrals, 

D v (E v (w)) = £ V (V)(F) + V(£„(F)) = V(If(^)) on £ L . 
Consider the operator equation 

D y o £ w = V o i F on £ L . 

If, hypothetically, a differential operator □: f — > x(7r) such that £ w o □: kerD y | f — ► 
kerD^I^. were constructed, which thus resembles the right inverse of £ w , then it would 
assign symmetries of a Liouville-type system £^ to arbitrary r-tuples of the integrals, 
see ffTOOl . 

In a recent publication [57J, the following fundamental result has been obtained: 

Suppose £l = {F = 0} is an m- component Liouville-type system that 
admits complete sets w G kerD y \ £ ^ and w G kerD^I^ of m integrals 
which are independent on a dense open subset of £l- Let an integer 
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k 3> 1 be sufficiently large, obeying an estimate in Lemma 4 of loc. cit. 
Consider the operator equation for □: f — > x(7r), 

£ w o □ = l mxm • D k x mod CDiff <fc (ker A,| £l -> kerD^J, (123) 

whose right-hand side is not fixed a priori. Consider the approximation 

fc— ordj; (w)— 1 

□ « a i- Dj-^W-* + 0(1) 

i=0 

of a formal inverse to £w\ and modify its columns following the construc- 
tive procedure of Lemma 4 in loc. cit. This yields the linear differential 
operator □ in total derivatives that solves equation (11231) exactly. 

The profound result quoted above establishes the truncation from below for the sequ- 
ence of coefficients of □ that satisfy some recurrence relations; most remarkably, this is 
a consequence of the presence of a complete set w £ ker D x \ of ^/-integrals for £ L . 

Whenever a local operator in total derivatives □ = (D 1 , . . . , □ r ) is known, the dif- 
ferential order of its z-th column LT can be decreased if the corresponding integral w 1 
is not minimal. Indeed, the reduction w % = D x (w l ) and the relation £^i = D x o £ w i 
means that the order k of the right-hand side of the equation D x o £ wi o = 
mod CDiff (ker D y \ £ — ► kerD^I^. ) has been increased artificially. 

Therefore the most general form of the determining equation (11231) is 

£ w o □ £ CDiff (ker A, | ^ -f ker D y | g J , (Tl2"3l ) 

with an a priori unfixed operator J^^ 1 ^^^, [w]) • D x ~ z in the right-hand side, here 
Qi E ker-D^I^, are (r x r)-matrices. Likewise, the operator □ itself is defined up to 
right multiplication by CDiff-operators from the same ring. As soon as the minimal 
integrals w are differentiated a suitable number of times such that their differential 
orders coincide in the end, and under assumption that the determinant det Qk does not 
vanish almost everywhere on £l, the topmost coefficient in the r.h.s. of (112311) can be set 
equal to unit, as in (11231) . Indeed, this is achieved by a solution □ := □ o Q^ 1 of (1 1231 p . 



The coefficients (Oj)^ of the product l& o □ are generally not equal to the Kronecker 
symbols 5^ u ■ 6*. Thus the operator □ o D~ h is not the right inverse of £ w . We conclude 
that equation (I123ip describes the problem of "right-inverting" operators in a ring with 
nontrivial zero divisors. A constructive algorithm for solving it has been given in [57] . 
However, in most cases the image of the resulting operator □ does not span the entire 
x-component of the Lie algebra sym^L, arid such operators are generally not Frobenius. 

Example 27. Consider the Liouville equation f^iou = { u xy = exp(2w)}, see its 
minimal integral w = u xx — u x , and the operator □ = u x + \D X that assigns symmetries 
to integrals. A straightforward calculation yields that 



'-2w 



o □ = Dl + AwD x + 2w x = l- Dl mod CDiff (ker D V L -> ker D V L ). 

■ ' ^-Liou I t-Liou ' 
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Example 28. Similarly, consider the 2D Toda system (I109p . its minimal integrals (IllOp . 
and the operator □ = (d 1 , D 2 ), which is given in fillip . Then we have 

ti wl o D 1 = Dl + w x D x + \wl 

£ w oD 2 = D x + w x Dl + [2wl - 3w 2 } ■ D x + [w l xx - 2w 2 x ] , 

£ w2 oD l = D X + w 1 D 2 x + 3w 2 D x + w 2 x , 
£i w2 o a 2 = Dl + 2w l Dl + 3w l x D 2 x + 

After the correction w 1 = w x , D 1 = D 1 o D x of differential order for the first inte- 
gral (j 11 (Jap , we obtain the matrix Q 5 = ( 2 |) with unit determinant at the top of the 
operator in the right-hand side of (I123ip . This yields a solution Dofl^ 1 of equation (j!23j) . 

We refer to footnote [221 on p. EH and to [3T] for further comments on this example, 
which is related to the Boussinesq hierarchy. Actually, the operator £ w o\3 is Hamiltonian 
whenever the second column, D 2 in □, is taken in its original scaling (here it has been 
divided by — §)• 

Lemma 30. If a solution □ of (I123ip for a Liouville-type system £l is Frobenius, then 
all coefficients of the bracket {{ , }}□ belong to keiD y \ £ . 

Proof. By assumption, we have that (D y o i$ o □) ([0', </>"]□) = for </>', 4>"(x, [w]) G f. 
This equals 

= (D y o eJ£ o □ ) (£□(»>) - Znw>0) + iW, ^»d) = (4 w) °n)(A/({{^On)), 

because the underlined composition satisfies (j 12310 . Clearly, D y ((f>') and D y (<p") van- 
ish on £l for arbitrary </>', <p" G f. For the same reason, not only the whole bracket 
{{(j)', </>"}}□, but each particular coefficient standing at the bilinear terms lies in ker D y \ £ . 

□ 



D x + 



w , 



|wL + w 1 wl 
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Conclusion 

We introduced the coordinate-independent definition of Frobenius linear differential 
operators A: Q 1 ^) — > q(tt) in total derivatives, whose images constitute Lie subalge- 
bras in the algebra g(7r) of evolutionary vector fields. We demonstrated that, in presence 
of Miura-type substitutions w: J°°(tc) — > there are two main classes of such oper- 
ators. They generalize recursion and Hamiltonian operators, whenever sections of fibre 
bundles in their domains fi 1 (^ 7r ) obey vector and covector transformation laws under a 
change of coordinates, respectively. We defined the notions of linear and strong com- 
patibile Frobenius operators and gave examples of compatible Frobenius structures of 
each kind, indicating both the matrix operators A and the types of transformations 
ip i— ► ip for pre- images ip of vector fields £a(^) • 

We solved a long-standing problem in geometry of integrable systems: We gave a 
complete description of higher symmetry algebras for hyperbolic Euler-Lagrange equa- 
tions £l of Liouville-type (in particular, for all 2D Toda chains associated with semi- 
simple complex Lie algebras of rank r). To this end, we obtained an explicit for- 
mula f ll05p for Frobenius differential operators □ of second kind that assign (Noether) 
symmetries <p = □(</>) to arbitrary r-tuples (p(x, [w]) (respectively, to <p = 5H/5w for 
any TC[x, [w])). Second, we calculated the commutation relations in sym£ L , and we 
proved that the arising coefficients of the brackets {{ , }}□ are differential functions of 
the substitutions w. 

We suggested a new algebraic operation on the linear spaces A = ®^ =1 K • A { of 
operators Ac. fi 1 (^ 7r ) — > whenever the spaces ^Vim/lj are closed under commu- 
tation. In particular, we considered the spaces of both linear and strong compatible 
Hamiltonian operators for the Magri schemes. Independently, we associated another 
class of such spaces to Liouville-type systems. Using the Lie algebra structure of vector 
fields 6^ with ip 6 imA^, we defined a Lie-type bracket L4 i; A,] = J2k=i^-k ° c^- of 
operators A^, Aj e A. This operation on linear spaces of Frobenius operators extends 
the associative composition and the Schouten bracket of recursions q(ti) — > q(ti). 

We have analysed the algebraic and geometric properties of Frobenius structures. 
We described flat connections in the triples (f2 1 (^ 7r ), g(7r), A) such that bi-differential 
Christoffel symbols are encoded by the structural constants c^- of the operator alge- 
bras A. We derived a condition when Frobenius recursions, which are solutions of the 
classical Yang-Baxter equation for the Lie algebra g(vr), generate chains of non- 
trivial deformations for the standard structure on 5(71). By our Counterexample HI 
the concept of Lie algebroids over the infinite jet spaces does not repeat the case of 
finite-dimensional base manifolds, when the anchors are determined by the Lie alge- 
bra homomorphisms A. Therefore we assigned formal differential complexes over Lie 
algebras n 1 (^ 7r ) to Frobenius operators A, such that the representations of the differ- 
entials in the complexes through homological vector fields restores one of equivalent 
constructions of Lie algebroids. 

Our approach not only incorporates the standard Hamiltonian formalism for PDE, 
but also generates new nonlinear models. A rigorous parallel between the bi-differential 
and classical affine geometry allows to interpret known equations of string theory as the 
master equations in the setting over the jet spaces. The suggested algebraic technique 
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can be especially efficient whenever off-shell considerations are in order, without any 
restriction onto differential equations. With all this, we intend to contribute to a pro- 
found relation between integrable systems and Lie algebras, which is well acknowledged 
in mathematical and theoretical physics [II [9j EH EUJ HH1 [52] . 

Discussion. Open problems concerning with Frobenius structures can be assembled in 
two categories: questions about geometry of involutive distributions generated by evo- 
lutionary fields in the images of Frobenius operators, and application of the formalism 
to integrable systems. 

Geometric aspects. Answering a motivating question of M. Kontsevich, who asked us 
whether Frobenius operators are the anchors in Lie algebroids over infinite jet spaces, we 
were stuck at Counterexample H] for A^ dV . We proposed a way out, which is based on the 
formal complex (j72p and the representation of its differential via the homological vector 
fields (in the spirit of the Koszul-Tate resolution for Lagrangian systems). At the same 
time, the difficulty means that a correct variational analogue of the Leibnitz rule ( f50l) 
has not been recovered yet. Variational analogues on J°°(7r) are known from [221 HQ] 
for vector fields, forms, commutators and Schouten brackets, for symplectic structures, 
etc. (see Table[Hon p. ED) but, to the best of our knowledge, the bare Leibnitz rule itself 
has never been under debate. Hence a profound analysis of Koszul-Schouten structures 
on Q k (M n ), which was performed in [33], remains without an analogue on 

^(k) = (g)V(^), 

leaving our ad hoc construction of Lie algebroids over J°°(7r) so implicit. 

We have already indicated that standard classification problems arise for the alge- 
bras A of Frobenius operators with the brackets [A^, Aj] = A^ocfj. For Hamiltonian 
operators P, we know two ways, ( H6|) and ( 155bj) . of calculating Dorfman's brackets [ , }p 
on fl 1 (7r). Is there a similar construction that yields the bi-differential Christoffel sym- 
bols rf-: fi 1 (7r) x Q 1 (tt) — > Q 1 ^) for totally compatible Hamiltonian structures? More 
globally, is there any homological obstruction for the existence of compatible Hamil- 
tonian or recursion operators that constitute non-commutative algebras? We expect to 
use here some techniques which will be more refine than the ones used in the proof of 
Lemma H21 We notice that the assertion of [23] was a reformulation of a much stronger 
result in [T3] that states non-existence of nontrivial variational ^-vectors and /c-forms 
at all k > 3 under given assumptions. However, this is not our case of the Christoffel 
symbols rf-, which are not skew-adjoint in each argument and which are not tensors, 
obeying the transformation rules for connection 1-forms. Let us indicate that, in these 
terms, the approach of geodesic motion is reproduced for the problems of optimization 
and control in the theory of partial differential equations. 

The use of coordinate-independent constructions of geometry of integrable systems [51 
[TTj [221 ED], see section [21 allowed us to introduce the well-defined notion of Frobenius 
operators. This would have been impossible within a 'simple' approach of local coor- 
dinates. We emphasize that now the definition itself manifests the two main classes of 
the operators, whose domains are the analogues of tangent and cotangent bundles to 
smooth manifolds. The general case of Frobenius operators and their Lie-type algebras 
will be the object of another paper. 
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Field -theoretic viewpoint. The definition of Frobenius operators does not depend on the 
number n of independent variables. In this paper, we passed from examples of Frobenius 
operators for ID evolutionary hierarchies (the KdV-type equations) to factorizations of 
symmetry generators for 2D Toda chains, which interpolate between KP and mKP 
(see [2n] )• We expect that Frobenius operators can serve a key to construction of 
integrable systems in multidimension. By this argument, we arrive at the problem of 
finding Frobenius structures for the flows of the KP hierarchy. 

Another recent construction may have applications in this theory, as soon as examples 
of Frobenius operators are known (e.g., for (1 + Tridimensional hydrodynamic chains). 
Namely, consider Frobenius complexes (1721) with the differentials d, and study their 
strong homotopy deformations, which involve Schlessinger-Stasheff 's algebras with Lie- 
type brackets of many arguments (see [H El EZj and references therein) . We know that 
in our situation, which is based on evolutionary fields on the jet spaces, the Wronskian 
determinants play a central role. Hence their homology-preserving generalizations for 
functions of many variables, which were defined in [27], will be useful. 

Returning one step back from the homotopy Lie algebras to homological vector fields, 
we recall their immanent presence in the Batalin-Vilkovisky formalism, see pQ. A 
description of its ingredients over the infinite jet spaces was our second motivation. 

In this jet setting, Frobenius recursion operators are solutions of the classical Yang- 
Baxter equation ([2]) for the Lie algebra g(7r) of evolutionary vector fields. It is still 
not clear how far one can pattern upon the standard approach to the Yang-Baxter 
equations, what systems £ C J°°(7r) can be obtained, are there any reductions for 
them, and what is the physical significance of the models whose state functions are 
transformed by derivations that belong to g(7r). 

Finally, we note once again that we have always remained in the non-graded setting 
and that the differential operators were local. We are advised to apply the formalism of 
Frobenius structures towards the study of noncommutative differential equations, see 
Remark [T31 on p. EETl and, first, reveal these structures for integrable evolution equations 
on associative algebras [50] • In addition, the use of difference Frobenius operators 
can be a fruitful intermediate idea for discretization of integrable systems with free 
functional parameters in their symmetries (e.g., Liouville-type difference systems). The 
restriction (CQ) produces narrow classes of difference operators such that the discrete 
systems are integral objects for the symmetry algebras. 

We have described a generator of operator algebras A with the commutation rela- 
tions (jlj). In section El we managed to calculate the bi-differential structural constants 
cfj by using the fact that the integrals w for the Liouville-type systems £l induce 
Miura's transformations w = w[u}: 03 — > 21 between bi-Hamiltonian hierarchies 21 
and *B C sym^L- By [60], many substitutions are obtained in this way. Similarly, 
the class of structures (J4j) given by Dynkin diagrams via the 2D Toda chains may be 
standard in the following sense: Is it true that for all spaces A of Frobenius operators 
Ai : fi 1 (^ 7r ) — > 0(71") of second kind there are Hamiltonian operators P : x(£) | — > x(£) 
such that {{, }} Ai = {{, }} P ? 

Final remark (M. Kontsevich, private communication). Algebras with bi-differential 
structural constants appear naturally in the BRST-formalism. Hence the algebras A of 
Frobenius operators may have applications in the string theory. 
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Appendix A. Analogy between Hamiltonian ODE and PDE 



In Table [T] we track the geometric correspondence between Hamiltonian ODEs and 
PDEs; we adapt the analogy to our needs, and therefore it is forced to remain incom- 
plete. The distinction between the coordinates and momenta in the PDE framework, 
which is implemented in section |5] to (symmetries of) Euler-Lagrange systems, is ad- 
dressed in [34, 25J. The concept of A-coverings over PDE, which is convenient in 
practical calculations that arise here, is developed in [22], see Remark [201 on p. [621 

Table 1 . Hamiltonian ODE and PDE. 



No internal structure of a time 
point t G R; 

Symplectic manifold M 2n 3 

(p, q); 

Components X % of vector fields 
X G T(TM 2n ) on M 2n ; 
Vector fields X G T(TM 2n ); 
Integral trajectories (p(t), 
q(t)) C M 2n of the field X; 
Lie algebra (TM 2n , [, ]) of vec- 
tor fields on M 2n ; 
The de Rham differential d; 

The de Rham complex; 

The space of Hamiltonians TC G 

C°°(M 2n ); 

The cotangent bundle T*M 2n 
and its sections if>: TM 2n — > 
C°°(M 2n ); 

The differential dTi of a Hamil- 
tonian H; 

Symplectic 2-form u G 

n 2 (M 2n ), Poisson bi- vector 

V G r(/\ 2 TM 2n ); 

Hamiltonian vector field such 

that X n _iw = dH; 

The Poisson bracket {Hi, 7^ 2 } = 



Smooth base manifold M n and a fibre bundle 
tt: E n+m -> M n . 

Infinite jet space J°°(7r) M n . 

Sections if = (^[u], . . ., v 9 " 1 ^]) ^ x(n) = 
F( 7T oo( 7T )) °f the induced fibre bundle. 
Evolutionary vector fields G D t, ( J°°(7r)). 
Solutions of the autonomous equation u t = (p. 

Lie algebra q(tt) of evolutionary derivations. 

The de Rham differential d = d^ + dc split to 

horizontal and vertical parts w.r.t. tt. 

The variational bi-complex and the C-spectral 

sequence. 

The highest horizontal cohomology H u (tt) 3 7i. 

The dual module f = Homjr( 7r ) (f(7r),A n (7r)) for 
a horizontal module of sections of a bundle over 

J°°(7T). 

Euler's operator E restriction of dc to 

H n {Tr), the 'gradient' E{H) G k. 

Hamiltonian operator A G CDiff(f, x) in total 

derivatives; f ~ for evolution equations. 



Sections tp = A(i/)), if) G x(7r). 

The Poisson bracket {7ii, 
A(E(H 2 ))) = £ A(E(W2)) (H 1 ). 



Appendix B. Reconstruction of the brackets {{•, - }} A 

In this appendix we describe an inductive procedure that assigns the bracket {{-, - }}a 
to a nondegenerate Frobenius operator A, see ( )62l) . The bracket {{ , }}a may be not 
contained in our knowledge that A is Frobenius if, e.g., the operator determines the 
factorization of symmetries of a Liouville-type system and has minimal differential 



62 



A. V. KISELEV AND J. W. VAN DE LEUR 



order. This is precisely the case of operator fl 1 U3j) . which is used in Example 1301 as an 
illustration. 



Remark 20. The algorithm we suggest is based on the use of the A-covcrings [22\ over 
the jet spaces J°°(7r). In our case, it amounts to 'forgetting' the jF(7r)-(sub)module 
structure of fi 1 (^ 7r ) and treating it as a jet (super-) bundle over J°°(7r), see [22J for 
details. Hence, instead of calculating D x {ip[u]) for ip £ fi 1 (£„.), one introduces the 
variable ip x and so on0 setting the derivatives in the Koszul bracket to zero, 

see ( 1621) . Consequently, only the derivatives of coefficients of the operator A contribute 
to the left-hand side of (I62al) while its right-hand side becomes A({{ipi,ip2}}A)- This 
is the structure of equations (I125H1261) below: they do not contain any evolutionary 
summands involving domains of A in the right-hand sides. 

Also, the nature of the assumption f)59p becomes clear. Indeed, this construction tells 
us that there are no notrivial kernels for restrictions of the nondegenerate operators onto 
the new jet bundles with sections ipi, ip2- Note that the same operators A may have 
kernels which are spanned by certain sections ip[u] £ r(7r^(£)) of the induced fibre 
bundles, see Example Hon p. [261 

For brevity, let us technically assume that A = \\ J2k^k ^x\\ * s a matrix opera- 
tor in D x , where A % 1 £ T{tt). We use the notation lj for the basic sections ip = 
*(0i, . . . , 0i_i, l u i+ i, . . . , r ) of ^(Ctt), whence 

r 

ij = J2^-h, ^£^(vr). (124) 

i=l 

Suppose further that the operator A is nondegenerate, see p. 

f|kerA fc = {0} for A = ^A k -D k x . (JMD 

k k 

Hence we encounter no difficulties when resolving inhomogeneous equations A^ipx, 
ip 2 }} A ) = Lp w.r.t. the brackets {{ , }} A £ CDiff(A 2 fi 1 (&.)), which are bilinear in 

Let the (yet unknown) bracket be 



r 



i,j,k=l 

where £ ^{n) and the condition c*j? h = — c£ k follows from the skew-symmetry of 
the bracket {{ , }}a- The coefficients c°° fc can be nontrivial if the dimension r of the 
fibres of £ is r > 1, see Example 1291 below. 

The base of the algorithm is given by the Jacobi bracket of the sections lj, lj: 

[A(l i ),A(l J )]=A(^ k -l k ). (125) 

The choice 1 < i < j < r yields r(r — l)/2 compatible systems of m equations. The 
components of sections in domains of A in the right-hand side of (j!25j) are enumerated 



24 Note that an additional relation ip y = can be introduced for the new variables ip that imitate 
the integrals. These integrals are the sections in domains of Frobenius operators (see ()100p and p07[) ) 
for Liouville-type systems. 
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by k. Since A is nondegenerate, the equations are solvable. Actually, these systems are 
overdetermined whenever the differential order of A is positive and hence the left-hand 
sides of (11251) and (I126p . see below, contain higher order derivatives of u that are not 
present among the arguments of c^[tt]. 

The inductive step is made by using the sections x a ■ lj and x 13 ■ lj. We obviously 
have 

[A(*^),A(xM,)]-a( Yl (* a ) ( ^V) W c£f ■!*) =A{a\j3\dff k -l h ), (126) 

0<a'+f3'<a+f3 

whence the coefficients c*j k on the diagonal a + (3 — const are obtained one by one. 
Having passed through the diagonal < a + (3 = const, with a > [3 or a < (3 in view 
of the relation c£ k = —c*- k , we check the condition 

[a{^a{ x )]=a[ <$-v&4) ( 127 ) 

a+/3<const 

that terminates the algorithm when holds. The differential order of the bracket {{ , }}a 
with respect to its arguments is estimated by calculating the Lie bracket [A(i/;), A(x)] 
and taking into account the Leibnitz rule in the right-hand side of (I62al) . We remark that 
the representation of jet coordinates u a using powers x a of base variables is standard 
in geometry of differential equations [5]. 

Example 29. The operator A\ for the dispersionless 3-component Boussinesq sys- 
tem (1561) has order zero and its matrix (I57cl) is nondegenerate almost everywhere, hence 
the condition (J59l) is valid. We reconstruct the bracket {{ , }}a 1 for this operator per- 
forming two steps of the above algorithm. 

The first step involves six combinations of ipx = 1, and ip2 = lj with 1 < i < j < 3. 
The second step repeats the first, but now tp2 is multiplied by x, and we have nine 
combinations ipi = lj and ip2 = x ■ lj, 1 < i, j < 3. Then the terminal check (11271) is 
fulfilled. This proves that the components of the bracket {{ , j-}^ are 

{{p, q}} u Al = v w x q u - v u q w x + Mp u qI - pI<?) + Mp v qI - pW) 

+ p u x q w - p w q u x + 2w x (p u q v - p v q u ) + u(p v q v x - p v x q v ), 

{{p, q}Y Al = p u x q u - p u q u x + Mp v q v x - pW) + p v x q w - p w ql + p w x q v - p v q w x , 

{{P, q}} W Al = u{p u q v x - p v x q u ) + 3w 2 (p v q: - plq") + 2u x (fq u - p u q v ) 

+ w(p u x q u - p u q u x ) + u{p v q u x - p u x q v ) + p w x q w - p w q w x . (128a) 

The first order operator (I57ep is also nondegenerate, and we obtain 

3w(p u q v x - p v x q u ) + 1w{p"q u x - p u x q v ) + 2w x {p u q v - p v q u ) 

+ p w x q u -p u q w x +Mp v q v x -p v x q v ), 
Mp v ql - p v x q v ) + p u x q u - p u q u x + 2(^<f - «), 
Hp v q v x - pW) + Mp v q u x - p u x q v ) + ^\p v ql - p v x q v ) 

+ w{p u x q u - p u q u x ) + 2u x (p v q u - p u q v ) + u{p u q v x - p v x q u ). (128b) 



{{p, q}} u M = 

{{p, q}} v M = 
{{p, q}} w A2 = 
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Sokolov's bracket for Hamiltonian operator (157dj) is equal to the sum of brackets (I128p . 



because the operators are linear compatible. The result agrees with 

Example 30. Following the above algorithm and using the package [47J, we obtain the 
bracket on the domain of the Frobenius operator fllU3j) : for any p, q G Q 1 (S[), we have 

{{p, <?}}□ = £ 2 ■ (Pxxq x - Pxqxx) - 2e 2 • (p xxx q - pq X xx) 

- l2e A ■ (8e 2 u 3 x u xx - 4e 2 u 2 x u xxx + 4e 2 u x u 2 xx + 2u x u xx - u xxx ) 

x [l + 8e 2 u 2 x + 16e 4 u x - 2e 2 u xx - 8e A u 2 x u xx ] 1 • (p xx q - pq xx ) 
+ (1 + 288e 4 i4 - 288e\ 2 x ii xx + 28e 2 u 2 x - 16e 2 u xx - 288e 6 u x u xx u 



X XX ^xxx 

- 96e 6 ul x + 3072e 10 u^° + 24e 6 u 2 xxx + 24e 4 u ix + 1408e 6 w^ + 3328e 8 u 8 x 

- 768e 10 U4 X u xx u x — 384e 8 u ix u 2 x u xx — 2304e 8 u x u xx u xxx + 384e 8 u 2 xx u x u xxx 

- 4608e l0 u 5 x u xx u xxx + l6e A u 2 xx - 5632e 8 u 6 x u xx - I920e 6 u xx u x + 3328e 8 u 4 x ii 2 xx 
+ 512e 6 u 2 xx u 2 x + 384s 1 V^Lz ~ 960e 10 u^ - 18 '»,//,,,,,. - 3072e 10 u 7 x u xxx 
+ 3072e w u 3 xx u x - 2304e 8 u 5 x u xxx - 576e 6 u 3 x u xxx + 288e 6 u 4x u 2 x + 384e 8 u 2 x u 3 xx 



+ 6U4e 10 u 2 xx u 6 x - 6U4e 10 u xx u 8 x + I152e 8 u 4x u x + 1536e 10 u 4x u 6 x + 192e 8 u 2 xxx u 2 x 
+ 240e 8 u xx + 1536e w u 2 xx u 3 x u xxx - 48e & u Ax u xx ) 

x [1 + 96e 4 u x + 256e 6 ^ + 2hQe 8 u 8 x + 4e 4 u 2 xx - 48e 4 u 2 x u xx + 32e 6 u 2 xx u 2 x 

- 4e 2 u xx - 25Qe 8 u 6 x u xx + 64e 8 u x u 2 xx - 192e 6 u xx ii x + I6e 2 u 2 x ] 1 • (p x q - pq x ). 

The two underlined units correspond to the bracket p x q — pq x on the domain of the 
operator □ = U x + \D X that factors symmetries of the Liouville equation U xy = exp(2W) 
at e = 0. In agreement with Lemma 1301 the non-constant coefficients of bilinear terms 
p xx q - pqxx and p x q - pq x belong to ker D y \ 



Remark 21. The classification problem for Frobenius operators A and the task of re- 
construction of the associated brackets {{ , }}a can be performed using any software 
capable for calculation of the commutators, e.g., [32], H7j designed for symmetry anal- 
ysis of evolutionary (super-)PDE. The implementation of technique of the A-coverings 
[22J is extremely productive here; see [32] for numerous examples. In this paper, we con- 
sidered not only the £- and £*-coverings, which correspond to Frobenius recursions (1751) 
for symmetries of PDE and to Noether operators for determined evolutionary systems, 
respectively. 

Appendix C. Bi-differential representations of the brackets {{ , }} A 

The components of Sokolov's brackets {{ , }}a € CDiff (/\ 2 fi 1 (£„.), r2 1 (^ 7r )) are matrix 
bi-differential operators in total derivatives w.r.t. the components of sections that belong 
to f2 1 (^ 7r ). We illustrate this for the brackets (11281) . which were calculated in Example [23 
on p. ESlfor Frobenius operators f)57cp and f)57el) . 

The notation means that the differential operators standing in the first and second 
tensor factors act, respectively, on the first and second arguments of the coupling (^i | 

{{,}}^2>. 
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The u-, v-, and w-components of the bracket {{ , on the domain k(i\) of A\ are 
the bi-differential operators 

3w ■ (1 ® D x - D x ® 1) + 2w x ■ 1 ® 1 D^ ® 1 - 1 ® D^ 

3u; • (1 ® D^ - D* ® 1) - 2^ ■ 1 ® 1 it ■ (1 ® D x - D x ® 1) 

Da; (g) 1 - 1 ® Da: 

'Da; (g) 1 - 1 ® Da; 

4k; • (1 <g> Da; - Da; ® 1) Da; (g) 1 - 1 <g> Da; 

Da; <g 1 - 1 ® Da; 

W ■ (D x (g) 1 - 1 ® Da;) U • (1 (g Da; - Da; ® 1) - 2u x • 1 ® 1 

U ■ (1 (g Da; - Da; (g 1) + 2u x ■ 1 <g 1 3u> 2 • (1 <g> Da; - Da; (g 1) 

Da; <g 1 - 1 ®D X/ 

The components of the bracket {{ , }}a 2 associated with the operator A 2 are 

3k; • 1 (g Da; - 2w D x <g> 1 + 2^ • 1 ® 1 -1 ® D x 

2w ■ 1 ® Da; - 3u> D x ® 1 - 2Wz • 1 ® 1 2n • (1 ® Da; - Da; ® 1) 

Da: ® 1 

'Da; ® 1 - 1 ® D x . \ 

iw ■ (1 ® D x - D x ® 1) -2 • 1 ® D x , 

2Da;®l / 

K) • (Da: ® 1 - 1 ® D x ) K • 1 ® Da; - 8li Da; ® 1 - 2lia: • 1 ® 1 0\ 

2lia; • 1 ® 1 + 8U ■ 1 ® Da: - U D x ® 1 (8l> + 6k; 2 ) • (1 ® Da; - Da; ® 1) . 

0/ 

Both bi-differential matrix representations are skew-symmetric in their arguments. 



Appendix D. Calculation of {{ , }} R for Frobenius recursions R 

The following program for Jets environment |47J under MAPLE calculates Sokolov's 
bracket {{ , }}r for Frobenius recursion R Q = A o Aj 1 , see (!HT|) . which is factored by 
the first Hamiltonian operator (I57ap and the inverse of Noether operator f)57cl) for the 
dispersionless 3-component Boussinesq-type system (!oT)l) . 

> read ' Jets . s ' ; 

> coordinates ( [x] , [u , v , w , pi , p2 , p3 , ql , q2 , q3] , 5) : 

The dependent coordinates p 4 , q J denote the components of the fibre coordinates £1,^2 
in the new jet space, which is the ^-covering over J°°(7r), see Remark [20] on p. [62] 
and [22J. The notation is in parallel with formula (1781) . 

> A : =Matrix (3,3): ia : =Matrix (3,3): M : =Matrix (3 , 3) : N : =Matrix (3,3): 

> A : =<<w*w_x+v_x , -3*w*u_x-u*w_x , u_x> I 

<-3*w*u_x-u*w_x , -3*w"2*w_x-4*w*v_x-u*u_x , v_x> I <u_x , v_x , w_x>> ; 

We have assigned A = Ai, see (I57cj) . 

> with (Linear Algebra) : 

> ia:=MatrixInverse(A) : 

This yields ia = A^ 1 : sym£^ — * cosym£, which we also denote by u. Now put m = 
A^(^) andn = A? fa). 
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> ml := simplify(ia[l,l]*pl+ia[l,2]*p2+ia[l,3]*p3) : 

> m2 := simplify(ia[2,l]*pl+ia[2,2]*p2+ia[2,3]*p3) : 

> m3 := simplify(ia[3,l]*pl+ia[3,2]*p2+ia[3,3]*p3) : 
> 

> nl := simplify(ia[l,l]*ql+ia[l,2]*q2+ia[l,3]*q3) : 

> n2 := simplify(ia[2,l]*ql+ia[2,2]*q2+ia[2,3]*q3) : 

> n3 := simplify(ia[3,l]*ql+ia[3,2]*q2+ia[3,3]*q3) : 

We have computed ipi = k>(£i); recall that ip G cosym£. 

Next, we apply the Hamiltonian operator Aq and obtain symmetries of (1561) . k = 
A (m) = MAlHti)) and 1 = i (n) = A (A^^ 2 )). 

> kl := simplify(evalTD(TD(ml,x))) : 

> k2 := simplify(evalTD(-4*w*TD(m2,x)-2*w_x*m2+TD(m3,x))) : 

> k3 := simplify(evalTD(TD(m2,x))) : 
> 

> ll:=simplify(evalTD(TD(nl,x))) : 

> 12:=simplify(evalTD(-4*w*TD(n2,x)-2*w_x*n2+TD(n3,x))) : 

> 13:=simplify(evalTD(TD(n2,x))) : 

Now we act by evolutionary derivations on the coefficients of the operator u = A^ , 
see (1751) . We set M = £roKi)(w); note that the matrix M is symmetric. 

> M[l,l] :=simplify(evalTD(kl*pd(ia[l,l] ,u)+TD(kl ,x)*pd(ia[l , 1] ,u_x)+ 
k2*pd(ia[l,l] ,v)+TD(k2,x)*pd(ia[l,l] ,v_x)+k3*pd(ia[l , 1] ,w)+ 
TD(k3,x)*pd(ia[l,l] ,w_x))) : 

> M[l,2] :=simplify(evalTD(kl*pd(ia[l,2] ,u)+TD(kl,x)*pd(ia[l,2] ,u_x)+ 
k2*pd(ia[l,2] ,v)+TD(k2,x)*pd(ia[l ,2] ,v_x)+k3*pd(ia[l ,2] ,w) + 
TD(k3,x)*pd(ia[l,2] ,w_x))) : 

> M[l,3] :=simplify(evalTD(kl*pd(ia[l,3] ,u)+TD(kl ,x)*pd(ia[l ,3] ,u_x)+ 
k2*pd(ia[l,3] ,v)+TD(k2,x)*pd(ia[l ,3] ,v_x)+k3*pd(ia[l ,3] ,w)+ 
TD(k3,x)*pd(ia[l,3] ,w_x))) : 

> M[2,2] :=simplify(evalTD(kl*pd(ia[2,2] ,u)+TD(kl,x)*pd(ia[2,2] ,u_x)+ 
k2*pd(ia[2,2] ,v)+TD(k2,x)*pd(ia[2,2] ,v_x)+k3*pd(ia[2,2] ,w)+ 
TD(k3,x)*pd(ia[2,2] ,w_x))) : 

> M[2,3] :=simplify(evalTD(kl*pd(ia[2,3] ,u)+TD(kl,x)*pd(ia[2,3] ,u_x)+ 
k2*pd(ia[2,3] , v)+TD(k2,x) *pd(ia[2,3] ,v_x)+k3*pd(ia[2 ,3] ,w)+ 
TD(k3,x)*pd(ia[2,3] ,w_x))) : 

> M[3,3] :=simplify(evalTD(kl*pd(ia[3,3] ,u)+TD(kl,x)*pd(ia[3,3] ,u_x)+ 
k2*pd(ia[3,3] , v)+TD(k2,x) *pd(ia[3,3] ,v_x)+k3*pd(ia[3,3] ,w)+ 
TD(k3,x)*pd(ia[3,3] ,w_x))) : 

In the same way, we define the symmetric matrix N = £^ (£ 2 )(u;). 

> N[l,l] :=simplify(evalTD(ll*pd(ia[l,l] ,u)+TD(ll,x)*pd(ia[l,l] ,u_x)+ 
12*pd(ia[l,l] ,v)+TD(12,x)*pd(ia[l,l] , v_x)+13*pd(ia[l , 1] ,w)+ 
TD(13,x)*pd(ia[l,l] ,w_x))) : 

> N[l,2] :=simplify(evalTD(ll*pd(ia[l,2] ,u)+TD(ll ,x)*pd(ia[l ,2] ,u_x)+ 
12*pd(ia[l,2] ,v)+TD(12,x)*pd(ia[l,2] ,v_x)+13*pd(ia[l ,2] ,w)+ 
TD(13,x)*pd(ia[l,2] ,w_x))) : 
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> N[l,3] :=simplify(evalTD(ll*pd(ia[l,3] ,u)+TD(ll,x)*pd(ia[l,3] ,u_x)+ 
12*pd(ia[l,3] ,v)+TD(12,x)*pd(ia[l,3] , v_x)+13*pd(ia[l ,3] ,w)+ 
TD(13,x)*pd(ia[l,3] ,w_x))) : 

> N[2,2] :=simplify(evalTD(ll*pd(ia[2,2] ,u)+TD(ll ,x)*pd(ia[2,2] ,u_x)+ 
12*pd(ia[2,2] ,v)+TD(12,x)*pd(ia[2,2] ,v_x)+13*pd(ia[2 ,2] ,w)+ 
TD(13,x)*pd(ia[2,2] ,w_x))) : 

> N[2,3] :=simplify(evalTD(ll*pd(ia[2,3] ,u)+TD(ll,x)*pd(ia[2,3] ,u_x)+ 
12*pd(ia[2,3] ,v)+TD(12,x)*pd(ia[2,3] ,v_x)+13*pd(ia[2,3] ,w)+ 
TD(13,x)*pd(ia[2,3] ,w_x))) : 

> N[3,3] :=simplify(evalTD(ll*pd(ia[3,3] ,u)+TD(ll ,x)*pd(ia[3,3] ,u_x)+ 
12*pd(ia[3,3] ,v)+TD(12,x)*pd(ia[3,3] ,v_x)+13*pd(ia[3,3] ,w) + 
TD(13,x)*pd(ia[3,3] ,w_x))) : 

We act by the operators M, N on £ 2 and £i, respectively, and calculate the difference 

e = £ ife(a)H(6) - £/?o(6)( w )(^i)- 

> el : =simplif y (M [1 , 1] *ql+M [1 , 2] *q2+M [1 , 3] *q3-N [1,1] *pl-N [1 , 2] * p 2-N [1,3] *p3) 

> e2 : =simplif y (M [1 , 2] *ql+M [2 , 2] *q2+M [2 , 3] *q3-N [1,2] *pl-N [2 , 2] *p2-N [2 , 3] *p3) 

> e3 : =simplif y (M [1 , 3] *ql+M [2 , 3] *q2+M [3 , 3] *q3-N [1,3] *pl-N [2 , 3] *p2-N [3 , 3] *p3) 

Next, we recall that Sokolov's bracket for A is equal to ( 157b jl . and substitute ipi = 
uj(£i) in it. Thus we put s = e + {{m, n}}^ Q . 

> si :=simplify(el) : 

> s2:=simplify(e2) : 

> s3:=simplify(e3+evalTD(2*(m2*TD(n2,x)-TD(m2,x)*n2))) : 

We can now check that formula (157b jl is correct. 

> J:=Jacobi([kl,k2,k3,0,0,0,0,0,0] , [11,12,13,0,0,0,0,0,0] ) : 

> Dl:=simplify(evalTD(TD(sl,x))) : 

> D2:=simplify(evalTD(-4*w*TD(s2,x)-2*w_x*s2+TD(s3,x))) : 

> D3:=simplify(evalTD(TD(s2,x))) : 

> simplify(evalTD(J[l]-Dl)) ; 



> simplify(evalTD(J[2]-D2)) ; 



> simplify(evalTD(J[3]-D3)) ; 



Finally, we act onto the right-hand side of (1751) by the operator u^ 1 = A x , and thus 
we obtain the components Z = A\{s) of Sokolov's bracket {{ , }}r . 

> Zl :=simplif y ( (w*w_x+v_x) *sl+(-3*w*u_x-u*w_x) *s2+u_x*s3) ; 

> Z2 : =simplif y ( (-3*w*u_x-u*w_x) *sl+ (-3*w*w*w_x-4*w*v_x-u*u_x) *s2+v_x*s3) ; 

> Z3:=simplify(u_x*sl+v_x*s2+w_x*s3) ; 

The result is somewhat a surprise: the output contains more than 15,000 lines. Be 
that as it may, the operator R = A o oj is the first known example of a well-defined 
Frobenius recursion for an integrable system. 
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